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TIME VARIATION OF SPATIAL AND ENERGY DISTRIBUTION 
OF NEUTRONS FROM A PULSED SOURCE 


I. G. Dyad'kin and P. Batalina 


Translated from Atomnaya Energiya Vol, 10, No, 1, 
pp. 5-12, January, 1961 
Original article submitted February 29, 1960 


The present paper is concerned with the time dependence of the spatial and energy distribution of 
neutrons from a pulsed source. Time-dependent problems of this type are encountered, for example, 
in geophysical methods of prospecting for oil. 

It is shown that the time distribution of neutrons of a given energy can be described by the 
Poisson probability distribution in practically the entire time interval. A space -energy-time rela- 
tion is derived from which it follows that the space -energy and energy-time distributions are inde - 
pendent in a certain range of distances and times, Outside this range there is a correlation between 
these relationships, and the correlation can be taken into account. An expression is obtained for 
the mean square slowing-down length as a function of time. All these relationships are given in 
the form of simple formulas having a clear physical interpretation. 


Formulation of the problem 
The time-dependent transport equation for the slowing-down of neutrons is of the form 


IQ grad +1) G(1, 0, th du’ \ f(QQ’, u—u’)G(A, r, Q’u', ty, t)= 
v oat 


(1) 
= 6(Q — Q))8(r —r,) (u—u,) (t— = 5°(1-— 0) 6(¢—¢,), 


G (1, 0, ty t)=0, 


where J is the mean free path, Q is the unit momentum vector;u = In(€9/e) = 21n(vo9/v) is the logarithmic neutron 
energy (€, v is the neutron energy and velocity, € 9, vo are the initial values of the energy and velocity); 0 represents 
r, Q,u; lm, QY, wy, etc.; and G(1, 0, ty, t) = G(r, Qyuy, r,Q, u, ty, t) is the Green function for the neutron colli- 
sion density. Furthermore, 


M 


where J yy is the mean free path relativeto nuclei of mass M. The operators in Eq, (1) act on r, Q, t only. 


The scattering is assumed to be symmetric in the center of mass system and J is taken to be independent of 
energy. 


An exact solution of Eq. (1) was given in [1] in the form of multiple integrals, In the present paper an approxi- 
mate solution of Eq, (1) is derived, using the solution obtained in [2,3] for the stationary case, Using the method 

described in [1] and [3], and regarding the term (2 /v) (8/0t) as a perturbation, we shall write the new integral equa- 
tion equivalent to Eq. (1) in the form 


G(A, 0, ty, t)=Gse(1, 0)8(t,—8) — (2, 0) 2 G(A, 2, ty t), 


where Gy (1, 0) is the solution of the stationary equation. 


It was shown in [1] that 


<1 (204-4) (26-44 
a, b=0 


(4, 0, 


where kp = k/k, @ap(k, u) are the elements of the matrix fully discussed in [2,3], and r and ry are in units of J. 
On expanding Eq, (2) into a Fourier integral in (r-%) and a double Legendre series in P, (ko) and Pj, (ky), and 


carrying out the Laplace transformation \ e “'~)d (t — t,), we obtain the following matrix integral equation which 


is fully equivalent to Eq, (1): 


G(k, 6, uy, u) = O(k, uw — uy) — \ du, xX O(k, u—uy)CG(k, 6, Uy, Uy), 


Uy 


c+ico 
1 
ak x do exp [ik (r + 0 (t —t,)] x 


c—too 


X [G(k, 0, Pa Py 


G (A, 0, ty, t)= > (2a + 1) (26-1) 
a, b=0 


and C is a matrix with the following elements: 


(C)an = (2a + 1) 


Let us now replace the exact matrix ® by the approximate expression [3] 


u—uy) = exp [(u —u,) B (a— V + X T + Ty (k) 8 (u—u,), 


where the matrices T and Toy have the elements 


The first element with the Legendre polynomial can be put equal to zero for convenience, 


Equation (8) holds both for u = uy andforu-u, > €,i.e,,inthe whole range of values of u except for a small re- 
gion near uy. Moreover, it was shown in [3] that & is practically independent of u, a and b. 


(4) 
q 
u 
a 
where 
6) 
(1) 
| 
(8) 
4 and 


Substituting Eqs. (9) and (10) into Eq, (8) we have 


[H+ E(k, 0, uy, = 7 (k) + (w— wy) — 


u 


_ \ exp [(u —u,)B(a—V k? +a?)] x T (k) CG (k, 6, uy, Uy) dug. 


This equation is the final integral-matrix equation which we have to solve, 


Thus, the problem is reduced to the solution of Eq. (11), which is equivalent to Eq, (1) for u not too close to wy. 
The solution will therefore be distorted only for the first few collisions, 


Solution of the time-dependent transport equation 
Equation (11) can be solved by the method of successive integration: 


G(k, 6, (a— Vk? + 0°)] x 
exp | — 27 (k) To (k) Ine | T (k) [#+Zer ny] "+ 


+ (RC ] (um), 


1 


where E is the unit matrix. This solution is not formal, It can be computed by taking logarithms and multiplying 
the obtained matrix from the left by 2TT9’. 


By taking the rank of the matrices successively as 1, 2, etc., one can use the Lagrange -Sylvester formula to ex- 
press Eq. (12) in terms of functions of numbers rather than matrices, 


The limitation of matrix T to unit rank does not mean the limitation of the scattering function fan’, u-u') 
to the first term in the expansion in the Legendre polynomials, It was shown in [2] that the first element of the 
matrix includes contributions due to all the coefficients of this expansion, There exist other, simpler methods of 
evaluating Eq. (12),and these will be analyzed below. They are based on the use of certain specific properties of the 
matrices T(k) and Tk). 


Neutron energy distribution 


If one integrates Eq. (6) over all the configuation space and velocity directions, all the sums and integrals will 
vanish and only the function (12) with k = 0 will remain, 


Using the results of [1-3] it is easy to verify that the matrix T becomes diagonal with the first element equal 
to &~' where € is the average logarithmic energy loss, and the matrix TgC becomes a unit matrix, They can there- 
fore be regarded as numbers and Eq. (12) will assume the form 


UNE 
G(k =0, 6, u,,u)= 


Assuming wy = 0, ty; = 0, vy = Vo, and transforming Eq. (13) with the aid of the formulas given in [4], we find that the 
neutron energy distribution is 


(11) 
13 
(144+) v 
3 


2 
where May) is the Whittaker function and ,F,(p, Vv, x) is the degenerate hypergeometric series, 


When t > 2 /(vo-v), i.e., for large t, it is necessary for the first few collisions (for neutrons with initial 
energy of 1 Mev, t >~ 107° sec) to use the asymptotic expression for (cf., for example [5)): 
BY ymp 


if, (—a, B, —) 
z>0, a>0O. 


We then obtain the following formula: 


We have thus reached a surprising result, namely, if one ignores time intervals which are small compared with t = 

= l/vo, then the time distribution is exactly the same as the well-known Poisson probability distribution for the 
appearance of 2/€ neutrons during a time t, subject to the condition that the probability of appearance of neutrons 
is the same for all instants of time and is equal to (v/Z )dt. . The time distribution is quite independent of the initial 


energy[ to within the time independent factor (1 ee, ] and is entirely determined by the neutron parameters 
% 


of the medium. 


Marshak was very close to this result in [6]. He guessed the existence of the formula exp(- b/x —x)x x M+ 1)° /4M 
for large x = vt/1 >> 1, However, since he did not succeed in obtaining a formula of the form of Eq, (14),he also 
did not obtain the simpler formula given by Eq. (16), which is already correct for vot/l > 1, i.e., for vt/l > v/vo. 
The latter is a very much less stringent limitation, 


Kazarnovskii [7] has obtained a numerical solution of the differential equation for the neutron density for large 
M. Numerical calculations based on our formula are in good agreement with his data, He also missed the relation 
given by Eq. (16) which has a clear probabilistic interpretation. 


Integrating Eq. (16) with respect to t between 0 and ~(which is equivalent to integration with respect to t, be- 
tween — © and t) we obtain the stationary neutron energy distribution 


Character of time changes in the spatial and energy distribution; mean square slow- 
ing-down length 

Equation (12) can be used to determine the space-energy distribution of neutrons, The properties of the ma- 
trices Tand To which enter into this formula have been described in detail in [1-3], The matrix ToC is not very 
different from the unit matrix qo(k)E for any k. The only large element in the matrix T is the element in the top 
left corner, We can therefore write these matrices down in the form of a main matrix and a small addition. Thus, 
for example, second-rank matrices can be written down in the following form: 


(14) 
) 


~ 3ikq, (k) 


o(k)C= B+ — ikq, ( k) — 4k2q, (k) E— 


(18) 


Yoo: — tkyo, _ —aikyo, 0 1 
=|, 


where 


If we now substitute these expressions into the matrix functions and expand them into the corresponding series 
(bearing in mind the noncommutative property), we shall obtain in the zero-order and first approximations 


0 @(Yoo—3Yo1) 
Vial 


oly"! 4 0 (19) 


_ 
ry qo(k) | 


Inverting the integrals as in the preceding section we obtain the zero-order Green function harmonic 


__ (t—t,)a sin k|r—r,|exp [(u—u,) B(a—V 


For larger instants of time one can again use the asymptotic expansion for ,F, and Eq. (19) can be integrated 
numerically, However, in order to obtain a clearer insight into the basic regularities, and to investigate the exist- 
ence of correlation between the space and time distributions (if any), we shall compute it approximately for t > 

> lAvo-v) and v << vo. 


(20) 


We shall take into account the zero-order harmonic only, since the angular distribution has been discussed in 
uetail in [3]. 


Assuming ty = 0, ry = 0, vy = Vo(uy = 0); t > Z/vo(s 107° sec) we find that 


tC sin kr exp [ ul 21) 


r 2a 
0 Jo (k V Lar | 


V(r, u, t)=5- 


| 

—1 
) 


The time-dependent factor reaches a maximum when t & 22qd¢k)/€v and then falls off exponentially, There- 
fore, for all values of t for which the integrand is practically different from zero we have 


pu(a—V 4-a2) > » k>a(+0.5). 


Conversely, for very small k the time-dependent factor varies more rapidly than the energy term, We shall use this 


in the following way. We shall replace the I’-function by the Stirling formula and rewrite the argument of, the ex - 
ponential in the form 


apu(1— 9/145 -)— In qo (k) +- (22) 


vt vt 2a aN 2a 2a 
Tao (k) Igo (ky EV + tV 


Let us rewrite this in the form 


B--A/ k? +03 
V +a, (23) 


4 2 vt 2a* 
2 vt— 2a3 
A=u+< (int 1-1) 2a 


For small and intermediate k Eq. (23) is identical with Eq, (22), while when k + © it slowly accumulates a 
relative error of a few percent but not exceeding 
ag 


(m2) 
2 ut 


pu+— 


when No(u,t) = 0. 


Let us now integrate Eq, (18) with the aid of the well-known expression (cf. , [8]) 


uo 


) Kal } 2 p V q-P-1 


The final result is 


_ ate” Ki (aV +P) 
N(r, u, eV (24) 


Since A2 2puA (2); <1 for all t for which N(r, u, t) is practically different 


from zero, this root can be expanded into a series and second-order terms can be neglected. Next, using the asymp- 
totic expansion for the Bessel function we obtain the space-energy distribution of neutrons from a pulsed source: 


N(r,u, Ng w) Ny (u, t)[1-4 e(r, 
where No(u, t) is given by Eq. (16) and 


(25) 


(26) 


where 
| | 
6 


is the stationary distribution discussed in detail in [3]. The values of coefficients & and 6 are also given in [3] for 
multicomponent media, The function [1 + €(r, u, t)) is not very different from unity for times of the order of the 


21 
effective pulse width (= = tmax + Ve —* ) and distances r < Bul (which roughly corresponds to the applicability 


of the age approximation): 


3 


(27) 
2 


Integrating Eq. (24) with the weight of r’, and using the formulas given in [8], we find the following expression 
for the mean square of the pales anes length 


2 3 


where iE is the slowing-down length for the stationary distribution and the second term represents a typical diffu- 


sion correction, When u >> 1, r is not very different from ie for all t for which No(u, t) differs from zero; T°/6 
increases monotonically with increasing t. 


The following important conclusions may be drawn from Eqs, (25), (26), and (28): 


1, The space-energy distribution of neutrons up to distances of the order of Bul is found to take the form of 
the product of the space (stationary) distribution and the energy-time distribution, In other words, the slowing-down 
of a neutron to a velocity v over a distance r from the source, and the time necessary for the slowing-down to the 
velocity v, are independent events and their probabilities can be multiplied. 


2, The mean square of the slowing-down length given by Eq. (28) is equal to the sum of the mean square of 
this quantity in the stationary case plus a diffusion type time-dependent addition which is proportional to t, The 
latter is small compared with r,, when v << vo for all instants of time for which the time distribution is different 
from zero, The presence of this additional term gives rise to a monotonic increase in the mean square slowing-down 
length with increasing t. 


3, Thus, the physical interpretation of the slowing-down process is as follows: At the end of a certain time 
interval (which is characteristic of the given medium) after the emission of fast neutrons, neutrons of a given energy 
appear simultaneously in all space except forr >> Bul. 


Their numbers reach a maximum at 


which is slightly dependent on r and then falls to zero, 


oVE 
The height of the maximum is roughly equal to 5) IV x © Nee(", u)s and the area under the curve is equal to 


the stationary neutron density, Thus the effective width of the "pulse" of neutrons with velocity v is 


x 
oVE 

The complete space, energy, angular, and time distribution of neutrons from a source of unit strength, located 
at an arbitrary point r, and emitting neutrons having a logarithmic energy wy in the time interval between ty and 
ty + dt;, is given by Eqs. (6), (12), (19), and (20), These formulas are required in the development of the perturbation 
theory. 


At eff = 


| 

| 
| 


Solution of a specific problem with a variable mean free path 

Equations (16), (25), (26), and (28) are so simple that they do not require a graphical interpretation. Their 
clear physical interpretation enables one to predict space-energy neutron density for media with variable mean free 
path. 


In fact, if the total density is represented as a product of the stationary space density and a time factor, then in 
the case of a variable path it is necessary to take the stationary density for the variable 2 as discussed in [3]. The 
value of 1 corresponding to the velocity v, i.e., 1(u), should be substituted into the time factor since for v << vo the 
latter is quite independent of vo and therefore completely “forgets” the initial energy and the first collisions, This 
will also apply to the first approximation, 
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Fig. 1, Distribution of neutrons in space and time for neutron 
energy of 1.5 ev in sand with various water concentrations, 


This approximation can be further improved with the aid of the formalism of integral equations for Green 
functions, 


Nonstationary neutron prospecting is beginning to be widely used in nuclear geophysics both in the USSR and the 
USA [9]. The stratum under investigation is irradiated with neutrons from a pulsed neutron source for a short interval 
of time and the neutron density in the stratum is measured after a further short interval of time, The density is very 
dependent on the amount of hydrogen in the stratum, Figure 1 shows computed curves representing the variation in 
the neutron density in wet sand (SiO: + m%H,O) at various distances from the pulsed neutron source, It is clear from 
this figure that for each porosity there is a characteristic time interval during which neutrons of given energy are 
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present. The dependence of tyyay, i.e., the time to the appearance of the maximum, on the effective neutron pulse 
width for the energy corresponding to the indium resonance (1.5 ev, vjng ™ 12 cm/sec) is shown in Fig, 2 for r = 0, 


For neutrons having different energies the values of tpyay and tog¢ 
given in Fig, 2 should be multiplied by ving/v. The beginning 
of the pulse (the instant at which the neutron density reaches e~ 
of the maximum value) corresponds to the point t & trax© aad 


4 


These curves can be used in the design of the corresponding 
apparatus in the choice of the distance between source and de- 
tector, and in the interpretation of experimental results, 
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FRAGMENT YIELD IN THE FISSION OF U?3* AND Pu?3® BY FAST NEUTRONS 


E. K, Bonyushkan, Yu. S, Zamyatnin, V. V. Spektor, 
V. V. Rachev, V. R. Negina, and V. N. Zamyatnina 


Translated from Atomnaya Energiya Vol. 10, No, 1, 
pp. 13-18, January, 1961 
Original article submitted April 16, 1960 


A radiochemical method has been used to determine the absolute yields for a number of fragments 
in the fission of U* and Pu? by 14.5 Mev neutrons and fission-spectrum neutrons, Some regulari- 
ties in changes in the fragment mass distribution are discussed. The dependence of the difference 


Am between the average masses of the heavy and light fragments on the mass number A of the di- 
viding nucleus is considered, 


The present work is a continuation of [1] which was concerned with the determination of the fragment yield in 


the fission of heavy nuclei by fast neutrons, It was carried out in order to augment existing data on the fragment 
yield in the fission of u™ and Pu, 


Most of the information on the fragment yield in the fission of U™ and Pu” refers to the fission of these iso- 
topes by thermal neutrons [2]. There is practically no information on the fragment sane in the fission of U** by 


fast neutrons, while data on the yield of some of the fragments in the fission of Pu” are available in a limited num- 
ber of papers [2-4]. 


The relatively limited amount of data available for U™$ and Pu’? is due to their high (as compared with a 
and U™*) specific -activity and the resulting difficulties in chemical treatment. In view of this, and in distinction 
to [1], in the present work we have used small specimens of fissile isotopes, and as large fast-neutron fluxes as 
possible, 

Measurements and results 

We have measured the yields of a number of fragments Ag, ca", 
cs! cs!87, Bal4? Ce!) in the fission of U** and Pu by fission-spectrum neutrons and 14,5 Mev neutrons, The 


experimental part consisted in neutron irradiation of us and Pu? specimens, separation and radiochemical purifi- 
cation of the required fission fragments, and the measurement of their 6 -activity. 


The us specimens (in the form of U3Ogs) and the Pu”*? specimens (in.the form of metal foils) had a weight of 
120-150 mg and were placed in hermetically sealed brass cassettes ~ 10 mm in diameter, 


The fission-spectrum neutrons were obtained from a us breeding system without a moderator, The cassettes 
were placed inside the system, The 14,5 Mev neutrons were obtained from the target of an accelerating tube, which 
was saturated with tritium and bombarded by ~ 150 kev deuterons, The uranium and plutonium specimens were 
placed in the immediate neighborhood of the target. The irradiation time both in the breeding system and in the 
accelerating tube was 7-10 hours, and the total number of neutrons entering a specimen during this time was ~ 5- ". 


The absolute number of U™* fissile events occurring in the specimen during its irradiation by fission-spectrum 
neutrons was determined with the aid of a small ionization fission chamber with a suspended layer of aad placed in 


the position of the specimen, In this way it was possible to obtain simultaneously the absolute yields of all the frag- 
ments under investigation, 


In the case of the specimens irradiated with 14,5 Mev neutrons, the absolute fragment yields of U™* were ob- 
tained by determining the absolute yields of Ba*° and Ce'! and the relative yields of the remaining fragments (rela - 
tive to Ba’*® and Ce), The specimens used in the absolute measurements were placed at a constant distance from 
the target (48.5 mm) during the irradiation, and the number of fissile events in the specimens was determined from 
the known neutron flux, measured with the aid of an &-particle counter. 


‘Fragment Yields in the Fission of u* and Pu? 


233 Pu239 


Fragments | Fission -spectrum 14,5 Mev neutrons Fission-spectrum 


neutrons 


14.5 Mev neutrons 


6,30-+0. 60 
4,754+0.35 
0.41340 .045 
0.16-0.02 
0.0837+.0.008 
0,056-0.006 
0.602+0, 050 


6. 1740.60 


5.1+0,8 
4.3540. 40 


In the case of the fission of Pu’*®, the absolute fragment yields were found by equating the total fragment 


yield as determined from their mass distribution curve to 200%, In the case of U™ this procedure served as an addi- 


tional check on the values obtained for the absolute yields, 
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Fig, 1, Fragment-mass distribution in the fission of Vv, oO: experimental,y, @: 


reflected points, 
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In order to separate the fragments from the irradiated u*™ and Pu”, use was made of the isotopic carrier 
method, All the fission fragments were separated from a single sample because of their relatively low activity. The 
separation of Pu’*® from the carrier was carried out by the precipitation of its peroxide from 3N hydrochloric acid 
with a preliminary removal of silver chloride. The separation of u™$ from the carriers was carried out by extraction 
from 6N hydrochloric acid, using a 15% solution of tributyl phosphate in benzene (also with preliminary removal 
of silver chloride), The separation and the radiochemical purification of the carriers was carried out with the aid of 
methods described in handbooks on analytical chemistry and radiochemistry [5-7]. For physical measurements, the 
specimens were taken in the form of PbMoO,, RuO2, AgCl, Cd2P,07, Te, CseSO4, BaCrO4, CeO2, and SrSO,, 
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Fig. 2, Fragment-mass distribution in the fission of Pu, V,O: experimental; V, @: 
reflected points; X; data from [3]; @: data from [2]. 


The number of fragments per sample was determined from their 6 -activity, measured with a counter calibrated 
against a 49 -counter, 


The fragment yields for U™> and Pu’? are given in the table and represent averages of sixty independent 


measurements, These values include the yields of all the members of the radioactive chain which precede the given 


isotope, and are practically total yields of isotopes of the given mass (Te'?9™ and Cs!* are the only exceptions), 


Figures 1 and 2 show the fragment-yield curves in the case of the fission of U™ and Pu? by 14.5 Mev neutrons 
(curve 1) and fission-spectrum neutrons (curve 2), For comparison curve 3 [2] shows the results obtained for the fission 
of U3 and Pu®® by thermal neutrons, The figures also show the yields of additional fragments: the reflected points 
found for Vy233 = 2.7 and Up,239 = 3,2 in the case of fission by fission-spectrum neutrons, and V = 4,5 in the case of 
fission by 14.5 Mev neutrons (the dependence of v on the fragment mass was neglected), 


Discussion of results 


The dependence of the relative probability of a symmetric fission of nuclei on the excitation energy for a 
dividing nucleus and its Z7/A parameter, has already been discussed in [1]. The latter paper contains data for ad 
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pom and as well as and (Fig. 3). The dependence of the position 
U _ a ih of the maxima in the fragment-mass distribution curve on the mass num- 
ber A of the dividing nucleus in the case of fission by fission-spectrum 
J neutrons was discussed in [1]. Analogous curves for the fission of uu". 
“ 7) s U8 and Pu? by 14.5 Mev neutrons, shown in Fig, 4, indicate that in this 


case also, the main effect is the change in the average mass of the light 
fragment, 


It is necessary to investigate in somewhat greater detail the depend- 
ence of the difference between the average mass of heavy and light frag- 
ments 4m (degree of asymmetry) on the mass number A and the excitation 
energy of the dividing nucleus, An attempt to establish a relation between 
Am and Z?/A was made in [8] and led to the formula 


Am = 0.090 [ (40.2 + 0.7)-4 4, 
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which was deduced from the linear dependence between (Am/ A)? and z?/A, 
This relation does not appear to be well founded since it follows from ex- 
Fig. 3, The relative probability of ond different isotopes of a (for ex- 
symmetric fission as a function of ample, , and that the average mass of the heavy fragment in 
22/A: O: fission by 14.5 Mev neu- these isotopes is approximately constant, 

trons, (®: points taken from [12)); 
V : fission-spectrum neutrons; @ : 
thermal neutrons [2], 


The present results, as well as other published data, show that there 
is a simple linear relation between 4m and A, which follows directly from 
the fact that the average mass of the heavy fragment remains constant. In 


this connection, the effect of the parity of the dividing nucleus on the mag- 
nitude of 4m remains somewhat obscure and it is difficult to explain it on the usual versions of the liquid-drop model. 


It is possible that this may be explained by the effect of the separate nucleon states on the fission process, and in 
particular, on the mass of the heavy fragment [9]. 
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Fig, 4. Fragment-mass distribution in the fission of u™ (light curve), us (heavy curve) 
u? ® (dashed curve) and Pu’? (dot-dash curve) by 14,5 Mev neutrons, 
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Figure 5 shows 4m as a function of A for the dividing nucleus in the case of fission by thermal neutrons and 
fission-spectrum neutrons, Line 1 is drawn through points corresponding to the fission of nuclei with even A, Points 
corresponding to nuclei with odd A are displaced from line 1 by two or three mass units (line 2) with the average 
masses of the heavy and light fragments measured to an accuracy of 0,2-0.5 mass units, In the mass range under 
consideration, the dependence of 4m on A is well represented by the expression 

Am = 288 — 1.04A+ 6 
where 


= { 0 for even A 
~ 2 for odd A. 
The small difference between the coefficient of A and unity is due to the fact that the average mass of the heavy 
fragment tends to decrease with increasing A. 


As was pointed out in [1], the degree of asymmetry in the fission depends also on the excitation energy of the 
dividing nucleus, The values of 4m in the case of fission-spectrum neutrons, by reactor neutrons and also thermal 


neutrons, are equal to within experimental error, owing to the small variation in the energy of the primary neutrons, 
An increase in the energy of the neutrons to 14.5 Mev leads to a reduction in 4m by 3 or 4 units. On the other hand, 
in the case of spontaneous fission the values of 4m are greater by approximately 3 units as compared with fission by 
fission-spectrum neutrons. 
It should be noted that in the case of fission by 14.5 
am Mev neutrons, and also in the case of spontaneous fission, it 
55 is found that Am depends on the parity of A, 
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‘ 

The “base point" method is used to determine the optimum temperature PP fora regenerative 
water-heating system which corresponds to a minimum calculated cost c, (kopeks/kwh) for elec- 
trical energy. It is shown that for an atomic power station having a water-cooled and water-mo- 


derated power reactor (WWPR) the value of Saas is near the maximum defined by the steam con- 
ditions after the first turbine stage. 


The optimum temperature er of a regenerative water-heating system is determined by engineering -economic 
calculations and corresponds to the minimum calculated cost for electrical energy, Ce (kopeks/kwh). An increase in 
t, and an increase in the number of heating stages reduces the fuel component of ce, but increases the cost of the 
regenerative system and the investment component, which may nullify the gain in c,. 


In atomic power stations with a constant reactor thermal power Q the temperature of the regenerative heating 
system in general influences the construction of the boiler, determining its surface F and the relation between the 
economizer and evaporator sections. However, the boilers used at atomic power stations with WWPR differ from the 


usual type in not having a separate economizer section [1-3], and the heating of the water from the temperature t, 
to t, takes place basically as a result of the partial condensation of steam which forms on the lower tube banks and 
of turbulent mixing. Therefore, in such boilers, in determining the heating surface F, the water temperature in the 
second loop is taken as constant and equal to the saturation temperature t,. Consequently, for constant reactor power 
Q, coolant inlet and outlet temperatures ty, and ty, and minimum temperature drive Atpjn, the calculated boiler 
surface F does not depend on the temperature t, of the regenerative water-heating system. This is the difference, in 
principle, between a WWPR boiler and a boiler with an economizer section. 


The author of [4] started with the incorrect assumption that the boiler in a WWPR atomic power station has a 
separate economizer section. This is what led to the incorrect conclusion as to the dependence of the temperature 
ty, on the temperature of the regenerative heating system. 


In a turbine at limiting power with constant steam consumption by the last stage, an increase in tr causes an 
increase in power as a result of the larger steam consumption by the first stages and the subsequent discharge of steam 


into the regenerative system. This is very important, since an increase in the power of a saturated-steam turbine is a 
problem of great interest at present [5,6]. 


The influence of the regenerative system temperature t, on the calculated electrical energy cost c, was inves- 
tigated by the base-point method [5]. All parameters referring to the base point with fixed regenerative system tem- 
perature t,' are marked with a single apostrophe, and a double apostrophe is used for other regenerative system tem- 
peratures. The choice of a particular base point has no effect on the final results of the investigation. 


The Effect of Regenerative Heating System Temperature on Power Station Efficiency 
and Unit Turbine Power 


In saturated-steam cycles the maximum possible regenerative heating system temperature e * is determined 


by the steam conditions after the first turbine stage , while the minimum temperature depends on the deaerator pa- 
rameters, 


| | 
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In order to extend the limits of variation of t, we took 
a scheme with low-pressure deaerators, 


i, Cal/kg The regenerative heating system temperature was in- 

675 vestigated for a turbine with an expansion process similar to 
that of the AK-70 turbine, the parameters for which are 
shown in the i-S diagram in Fig.1. 


650 


It is known that the heating system temperature t," 
m Vr should be related to the number of regenerative heaters z,". 
600} 599.0\% 4 For t, = constant, we have a definite minimum number of 
heating stages corresponding to the maximum of the 
575 “x-0. efficiency curve n =f (t,» Zp). When Z, < an increase 
\ } in t, leads to a decrease in 7, which is inadmissible. There- 
550 , fore, for an increase in t, the. minimum number of heaters 
should be increased. In real schemes the number of heating 
stages z, is always greater than the minimum z,"", since 
increasing z, to infinity increases the efficiency of the power 
Fig. 1. Steam operation in a type AK-70 turbine station. However, the investment cost for the regenerative 
(a.at.: absolute atmospheres). system limits the increase of z. When faced with the prac- 
tical selection of the number of heating stages, one should 
keep in mind the fact that when the number of heaters is greater than three or four, the addition of each new stage 
adds less and less to the efficiency. However, for a fixed feed water consumption and for t, = constant, the fewer the 
heating stages z, the greater must be the surface of each heater. The number of heating stages z,, therefore, must 
not be taken smaller than the number for which the heater surface F exceeds the value for heaters manufactured by 
industry at the present time, For F = constant, an increase in t,, as well as an increase in feed water consumption, 
leads to an increase in the number of heating stages. To ensure uniformity of equipment, all high-pressure heaters 
(HPH) and all low-pressure heaters (LPH) must have equal surfaces. 
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Fig. 2. Heat flow schematic for the second loop: 1) boiler; 
2) high-pressure cylinder; 3) low-pressure cylinder 4) steam 
separator; 5) generator; 6) condenser; 7) pump; 8) condensate 
heater, using steam from the ejectors and turbine effluent; 

9) LPH; 10) deaerator; 11) HPH (a. at.: absolute atmospheres). 


Considering the maximum surfaces of the regenerative manufactured by the Khar’kov Turbogenerator Factory 
(KhTGZ) (F = 200 m? for LPH and F = 490 m? for HPH) a regenerative water-heating scheme was computed using a 
constant water subcooling value of 5°C and a condensate supercooling value of 3°C, for the volume of steam to the 
last AK-70 turbine stage equal to G; = 285.6 tons per hour (Fig. 2). The temperature t,' = 108°C, corresponding to 
the water temperature after the deareator, was taken as the base point. 
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The heating temperature was analyzed for a change in t, from 108.7 to 206.4°C; the number of LPH's was kept 
constant. Such a scheme satisfies the general thermodynamic principle that with an increase in t,” the optimum num- 
ber of heating stages is increased, 


The volumes of steam alloted to the LPH's 4+ az, the deaerator ag, and the HPH's a, + ag + Og + Og, the sepa- 
rated moisture content d,¢,, and the surfaces of the heaters are shown in the table. The table also gives the values 


of unit steam consumption & (kg/kwh) calculated for a generator efficiency of 0.98 and a turbine mechanical effi- 
ciency of 0.93. 


For a constant exhaust area the relative change in turbine unit power N-'/Ny" may be expressed [6,7] by the 
formula 


Ni dg(i1—Xa") 
Ny dg (1) 


where 1—f ca is the relative volume of steam arriving at the last turbine stage. 


Power Station Characteristics for Various Regenerative 
Water Heating Temperatures t, 


Characteristics | | 155 206.4 
/ 
No. of HPH's o} 3] 4 108 | 
a, 2.22) 5.16 | 4.91} 4,63] 4,33 7 
ay 3,36] 3.02 | 2.87] 2.22) 2.7 > | 
10,7 |10.15 | 9.15] 8.6 % 1.04 
ag 3.61] 3.12] 2.66] 2.22] 1.73 — 
— | 5,39] 5.14) 4.87) 4,59 nin 
ay | 5.4 | 4.78] 4.43 
77.1 |72.1. | 69.6 | 65.9 | 61.9 
de, kg/kwh 5.44 - et rp By Fig. 3, Relative changes in power station efficiency, 
— |340 [352° [394 in turbine unit power, and in unit steam 
— |350 |369 {392 consumption, d3/d3, as functions of regenerative water 
Fs, m, heating temperature ty. 
Fy, m 435 | 


Note, The’values of and are given in | 


The relative changes N_"/N_' in turbine power and 
percentages, 


erp. 8: in power station efficiency are shown in Figure 3. 


Relative Calculated Cost for Electric Power Sco and Relative Unit Investment Cost 
6q, as Functions of t 


The basic equation for the analysis of relative calculated cost Ce = Gt St Sp by the base point method for 
constant reactor thermal power is of the form 


(2) 


where CAP and oF - Cop are the relative components of calculated cost at the base. point; 6n = 
=7 "5/0 "st ~ 1 is the relative change in power station efficiency: Sq, = q"./a\ — 1 is the relative change in unit 
investment cost; Scop is the relative change in the operating component. At atomic power stations the operating 
component is small, and in calculation by the base method any change in it can be neglected. Taking 5cop = 0 
and Cop = cyp/ce = 0.1, using Eq. (2) we obtain the following simplified relation: 


cf dn + (0.9 —c}) (3) 


The relative change in the fuel component dc, = s /¢ -1 in the case under consideration is equal to the rela- 
tive change 6n in efficiency (Fig. 4). 
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The relative change in the investment component 6q = 
or/¢ - 1 is equal to the relative change in unit investment cost, 
which for constant investment cost for the first loop (C,, = con- 
stant} may be represented in the following manner; 


+ 


¢= [1 — bmr(1 — ner — + 
+ bine (S¢mr + 


(4) 


where and 6 are the relative changes in unit invest - 

ment cost for the machine room and for the industrial water sup- 

a ply: bh, =a’ / q is a basic constant equal to the ratio of the 

.. machine room investment cost to the investment cost for the en- 

tire station; = qi /q' and ny, = qa /q'__ are basic com- 

r ~ ponents equal to the ratios of the industrial water supply and elec- 
00 120 140 160 180 t,,°C trical equipment investment costs to the investment cost for the 

machine room. 


Relative change, % 


' 


Fig. 4. Relative change in the fuel component 


5c, andthe unit investment cost 6q,,, for the The ratios of industrial water supply and electrical equip- 
power station machine room as functions of t,. ment investment costs to the machine room investment cost have 


average values of ni, = 0.25 and n2, = 0.15, and we therefore 
have from Eq. (4) 


= (1 — + (89mm + 9-25 (5) 


The unit investment cost for the machine room, calculated for t, = 108°C, may be subdivided into a number of 
components (for an upper limit to the cost of the regenerative system); 


Component of investment cost Amount 
n' (turbine) 0.365 
tur 
n gen (generator) 0.315 
n' (equipment 0.085 
eq equip ) 
(structural parts) 0.105 


n (regenerative system) 0.025 


n' (condenser) 0.105 
con 


n (total) 1.00 


Then the relative change in unit investment cost for the machine room may be expressed by the formula 
” n q; 
q i 
The unit investm ent cost for the condenser Yon iS Proportional to the unit surface {,,,, required for 1 kw of elec- 


tric power. From Eq. (1) for the case of a constant vacuum we obtain the relation 


fon _ ‘oon, 20) Ne 


Icon Toon 


In calculating the unit investment cost for the turbine the exhaust area is taken as constant. Therefore, if 
there is a change in t_, there is a corresponding change in turbine power (Fig. 3). 
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The saturated-steam turbines manufactured by the KhTGZ have a high-pressure cylinder HPC and a low-pres- 
sure cylinder LPC. Since the length of a turbine blade in the last stage is constant and the amount of steam ay + a, 
separated for regeneration from the LPC changes little (see table), it may be assumed that the costs for the LPC are 
constant (C, = const). According to the data supplied by the KhTGZ the cost for the LPC is 1.5 times the cost for 
the HPC, and the cost for the turbine blades amounts to about 20% of the cost for the whole turbine. Then, assuming 
that the cost of the HPC blades is directly proportional to the blade length in the first stage, we obtain 
four = “Np 
Tur pt A hp ) Ny 
= (0.08 = Oi 


As can be seen from Fig. 5, an increase in t, pecan little fade in the turbine investment cost. 
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Fig. 5. Relative change in unit investment cost 
for the system and the turbine “3.05 720 140160. 780 t,, 


as f 
as functions of t,. 


Fig. 6. Relative change of unit investment cost 

5q_ as a function of t, and of the fraction of invest- 
The changes in unit investment costs for the ment cost spent for the machine room b;,, 
generator, structural parts, and auxiliary equipment 
can be neglected. 


The cost of one HPH as a function of its surface area, equal to 300-450 m’, is defined by the approximate for- 
mula 


Cypp = 32 + 0.245F thousand rubles. 


The data on the effects of changes in t, on the surface area of the HPH are shown in the table; in the calcula- 
tion the coefficient of heat transfer k was taken equal to 2100 Cal/m?- hour: °C. The unit investment cost for the 
regenerative system changes with an increase of HPH'sfrom 0 to 4 as follows: 


4 
[ApH + 3) (32-4-0.245F)] y. 


(9) 


K LPH Nr 
The dependence of the unit investment cost for the heaters on the regenerative system temperature for G pH = 


290,000 rubles is shown in Fig. 5, from which it can be seen that an increase of t, from 108 to 206°C increases a, 
more than 2.5 times. 


The relative change in the unit investment cost for the entire machine room tn , calculated from Eqs. (6-9), 
has a minimum value for t, = 130°C, and for higher t. values the unit investment cost for the machine room increases 
as a result of the greater expense for the regenerative system (see Fig. 4). 


The unit investment cost for the industrial water supply decreases with an increase in t,: 
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The dependence of 64. and 6c, on the temperature t, is shown graphically in Figs. 6 and 7. 


Evaluation of the Results 


From Fig. 7 it can be concluded that 6c, decreases for each version with an increase of t, up to the maximum 
possible value. Consequently, for the WWPR atomic power station under discussion the optimum temperature ~ 
any fuel component corresponds to the maximum temperature of regenerative water heating — max defined by the 


steam conditions of the first turbine stage. In the present article we consider an upper limit for the investment cost 
for the regenerative system n= 


4/43 ; obviously, for a cheaper regenerative system it is even more profitable to in- 
crease tr opt 


for 
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Fig. 7. Relative change in calculated cost for electric power 6c, as a function of t,, relative fuel 
component G and fraction of investment cost spent on the machine room bypr: 
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The conclusion thus reached is not at all typical, and for this reason we have devoted considerable space to the 
reasons for such a result. As has already been indicated, the optimum regenerative heating temperature t.*  corres- 
ponds to the minimum electric power cost c,, that is, to the minimum of the sum c, + ¢,. Ast, increases, the fuel 


component ¢, always decreases, and therefore the minimum of c¢, depends on the nature of the change in the capital 


investment component q.. In the case under consideration, an increase in t, essentially also produces a decrease in 
the investment component om (see Fig. 6). 


The decrease of G and q, with an increase in t, caused the optimum heating temperature t, Pt to coincide with 
max 
the maximum t, 


Therefore, in this case it was profitable to use the maximum possible regenerative water heat- 
ing temperature, taking into account, of course, the engineering possibilities of the construction of a rational scheme. 
As a result of the concavity of the 5c, = f(t,) curve, shown in Fig. 7, it is possible to reduce somewhat the calculated 
value of t, OPE so as to decrease the unit investment cost (see Fig. 6) and the initial cost for the power station K, =q,Ne- 


20.3 

1 EN 

CON 

-2 

00 120 40 180 200 0 40 160s 2008, © 

| 

20 


LITERATURE CITED 
S. A. Skvortsov ,Atomnaya Energiya, 5, No. 3, 245 (1958). 
G. V. Ermakov, Teploénergetika, No. 10, 88 (1957), 


J. Simpson et al., Collection; "Atomic Power Engineering” (Russian translation] Moscow, Gosénergoizdat, 
1956, p. 487. 


1. A. Trub, Atomnaya Energiya, 4, No. 3, 286 (1958). 
Yu. D. Arsen'ev, Teploénergetika, No. 10, 27 (1959). 
Yu. D. Arsen'ev, Physics and Heat Technology of Reactors [in Russian]. Supplement No. 1 to the journal 


"Atomnaya Energiya” Moscow, Atomizdat, 1958, p. 176 [English translation, Consultants Bureau, New York, 
1960). 


L. I. Kertselli and V. Ya. Ryzhkin, Heat-powered Electric Power Stations [in Russian] Moscow, Gosénergoizdat, 
1956. 


2. 
4, 
5. 
6. 
21 


BOUNDARY CONDITIONS IN THE METHOD OF SPHERICAL HARMONICS 


G. Ya, Rumyantsev 


Translated from Atomnaya Energiya, Vol, 10, No, 1, 
pp. 26-34, January, 1961 
Original article submitted March 9, 1960 


The equations of the method of spherical harmonics are derived for use with the single-velocity 
transport equation, From an analysis of the equations, conditions are obtained for the boundary 
between two different scattering media, It is shown that these conditions completely define a solu- 
tion in any-order approximation (we are concerned here with the so-called Py approximations) in- 
cluding even-order approximations, The possibility of using even-order approximations is of major 
practical interest. Thus, for example, the relatively simple P, approximation may serve as a 
method for improving the elementary diffusion theory without introducing excessive complexity. 


Even-order approximations, and in particular the P; approximation, have not as yet been 
used in calculations, 


Introduction 


The transport equation can be successfully solved with the aid of the method of spherical harmonics, This 
method consists in that the required distribution function F (r, Q) is expanded into a series of spherical functions, The 
expansion coefficients which depend only on the space coordinates, are related to each other by an infinite set of 
differential equations, If one limits one's attention to spherical functions up to N-th order only in the expansion of 
F(r, Q), then the set of equations for the coefficients will be finite, and can often be solved analytically. This is 


known as the Py approximation, The uniqueness of the solution is ensured by providing a sufficient number of bound- 
ary conditions, 


The method of spherical harmonics is frequently used in the solution of one-dimensional problems [1, 2]. The 
derivation of the equations for the Py approximation has been given by Davison [3] for the general case, However, 
the problem of the formulation of the boundary conditions in the method of spherical harmonics cannot be considered 
as entirely solved. This is true not only for the general case but even for the one-dimensional problems with simple 
geometry. This is probably the reason why it is generally believed the P2, Py, etc, approximations, i,e,, the even- 


order approximations, are inferior and less convenient for the solution of the equations than the odd-order approxi- 
mations, 


Let us state the problem more carefully. In the selection of the boundary conditions on the boundary between 
two media, one must distinguish between two cases. In the first case, the Pyy approximation equations cannot be 
applied to one of the media (for example, a semi-infinite vacuum), or they can be applied to both media but do not 
hold on the boundary itself and in its neighborhood, i.e,, the Py approximation serves only for the description of the 
asymptotic behavior of the solution, In this case, the boundary conditions are known to be artificial and hence not 
unique. Their choice should be dictated by considerations drawn from some other, more accurate method for the 
solution of the transport equation, or by the physical interpretation of the problem, 


In the second case, the Py approximation equations are assumed to hold throughout both media, including the 
separation boundary. This assumption can be admitted if the adjacent media are good scatterers and their physical 
properties are not too different. Such conditions are relatively frequently encountered in heterogeneous systems, 
Here, the boundary conditions are required only because it is not possible to express the solution by a single analyti- 
cal formula which would hold for the entire heterogeneous medium, The solution is written down in the form of 
separate functions each of which is a solution within a single homogeneous region, The boundary conditions, i.e., 
the conditions for mating these functions, should satisfy the requirement that the general solution made up of these 
functions should also hold at the mating points on the separation boundary, Here the conditions at the mating points 
cannot be arbitrary and are a unique consequence of the form of the equations. 


: 


Let us consider this type of boundary condition, 


Derivation of the equations of the method of spherical harmonics 

Consider a single-velocity transport equation for a medium of arbitrary geometry. Suppose the function F(r, 
Q) describes the spatial and angular distribution of the neutron current, Let us introduce a vector distribution function 
defined by 


F(r, Q) = QF (r, Q), 
where Q is the unit vector in the direction of the velocity. 


Clearly, the function F(r, 9) is the vector neutron current at the point r, Q while the function F(r, Q) is the 
absolute magnitude of this current, i.e. 


F (r, Q)=|F(r, Q)|. 


Suppose further that W( p19) is the probability that the neutron will be scattered and undergo a change of direc- 
tion from ' to Q, and 


= (Q’Q). 
The density of neutron sources in the medium will be noted by S(r, Q). 


Using this notation, the transport equation is of the form 


div F (r, ZF (r, F (r, W dQ’ = S(r, Q). 


Let us expand the functions F(r, 9) and F (r, Q) into a series of spherical functions (cf., [4]): 


F (r, Q)= >) F, (r, Q); 


n=0 


F (r, Q)= (2n+ 1)Y,, Q), (3) 


n=0 
where F(t, Q) and Y(t, N) are,respectively, the spherical functions with vector and scalar coefficients depending 
onr. 


Each spherical function consists of 2n + 1 linear independent components: 


The choice of the form of Y,,,,(Q) is immaterial. 


The source density can be represented in a similar way: 


S(r, Q)= > S, (r, Q), 


n=0 


and the probability W(19) will be expanded into a series of Legendre polynomials: 


o 
(2) 
(5) 
| 
(6) 


WV (1g) = (2 + 1) (Ho): 


n--0 
From the normalization condition Cg = 1. 


Bearing in mind that the spherical functions have the following properties [4], 


Py (ty) dQ’ =0 for kn, 


an 


(2) Py (14g) d2’ = (Q), 


4n 


we find that 


Wats w d2’ = (2n-+ 1)C,Y,, (1, Q). 
n==0 


An 


Since the spherical functions of various orders are linearly independent, Eq, (1) can be rewritten in the form 


div (r, 2) + Bal (r, Si (r, 2), (8 
O<n< ) 


In this expression 


It is now necessary to establish the relation between F,, and Y,. In accordance with the definition of F(r, Q) we have 


F (r, 2) = >) (2n+ 1) QY, Q). (10) 
n=0 
It is known from the theory of spherical functions that® 


Y, (2) =U, 


where p is the radius-vector in the three-dimensional space 1, Hy Hz, i.e., 


where = y, Hy = pQ y’ =UQ, and is the n-th degree harmonic polynomial in Hy. Hy, Hy. 


Using the Hamiltonian operator 


it follows from the definition of the polynomial U,, that it should satisfy the Laplace equation 


(a) = 0. 


We can also write 


The product #U_, is a homogeneous polynomial of degree n + 1, although it is not a harmonic polynomial, It 


follows that when pt = 1 it does not become identical with a spherical function of degree n + 1, Let us represent it 
by the following identity: 


* This mathematical transformation is derived in a chapter of Davis’ book [3]. 


| 
| 
| : 
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It can be immediately verified that | wU,.- — it vU =| and VU, satisfy the Laplace equation® and at 


the same time are homogeneous polynomials of degree n + 1 and n~1, respectively. This means when yp = 1, the 
first expression forms a spherical function of order n + 1 and the second, a spherical function of order n~-1, 


On rewriting the expansion (10) in the form 


n=0 


(the arguments r and Q will be omitted from now on), and selecting from the above sum spherical functions of order 
n, we have 


F,, = (VU + (2a — 1) aU, (12) 


Next, let us introduce the further Hamiltonian operator 


The equations for the spherical harmonics will then assume the form 


((VV)U + (2n — 1)(Va) (VV) Uy Juni +2,Y,=8 (13) 


n° 


We note that (VV) and (Vp) are scalar operators which, from the formal point of view, are the scalar products 
of the corresponding vector operators, i.e., 


(VV) = (VV) = Tay on, 


(Va) = (WV) = 


In order to remove the condition p = 1 it is sufficient to multiply Eq. (13) by up", since each term of this equa- 
tion is a homogeneous polynomial in‘p,, Hy Hy. At the same time, let us.replace n by n - 1 so that 


(VV) = ,_1—(2n—3) (Va) U, +2? (VV) (14) 


where Qn = us, are harmonic polynomials of degree n. 


Using Eq. (14) for n—2, n—4, n-6, etc,, we can eliminate the operator (VV) from the righthand side of Eq. (14) 
and rewrite it in the following form: 


(15) 
(Vu) [(2n 3)U otp? (2n a5 7) U, at (2n 11) ] 


(terms with negative numbers are zero), 


* This can be verified with the aid of Euler's identity (pV )A,=0 
degree n. 


%, where %, is a homogeneous polynomial of 


@ 
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In passing on to the Py approximation we shall assume that (VV)Uy , 4 = 0. In that case, Eq. (15) with the 
maximum number n = N + 1 assumes the form 


n= [On + + p'Qn-4+ 
((2V — 1) +p? (2N —5)Un-at (16) 
— 4+ ...]. 
Let us now consider the method of solution of the set of equations derived above, and determine the number of 
arbitrary coefficients in the general solution, According to the theory of linear equations, the number of coefficients 


is equal to the number of proper solutions of the corresponding homogeneous system, For this reason we shall deter- 
mine the form of the solution on the assumption that the terms Q, are absent. 


Equation (16) gives the function U, in terms of lower-number functions, All the terms of the equation are 
homogeneous polynomials of degree N, and we shall therefore refer to it as an N-th order equation, If we apply the 
operator (VV) to the whole equation, the order of the equation will be reduced by one, It is easy to verify that the 
following operator relations will hold: 


(VT) (Vp)! = (Vey! (WV) + 


= (VV) Qhp2 (Ppa), 


and consequently 


(VV) = (Wp)! (VV) 
Thus, as a result of the application of the operator (VV) to Eq. (16), no new scalar operators other than Vv? = A can 


be obtained, and the operator (VV) can always be excluded from the final expression with the aid of Eq. (15). We 
note that 4 commutes with all the other operators in the equations, just as (Vp) commutes with Th : 


If we repeat the operation (VV) on Eq, (16) m times, we obtain a new equation whose structure can easily be 


deduced, Firstly, it should be of order N—m, and secondly, it will only contain functions of order N~m or lower, 
since each time the operator (VV) is eliminated with the aid of Eq. (15), the function with the largest number is 
also eliminated, In view of this, the form of the equation should be 


Ln—-mUNn-m = Qn—m [(Vp) Un—m-1; p?U etc, 1. (17) 


The operator Ly _ ,, cannot comtain Vp or u” since this would be in conflict with the order of the equation, At 
the same time, this operator should include differentiation with respect to r m times if m is even, and m+1 times if 
m is odd, This is due to the fact that the operator (VV)™ is equivalent to differentiation with respect to r m times, 
and the functions N\y.,, Uy_g, etc., in the original equation (16) are already being operated on by a multiple 
differentiation operator with respect to r, However, differentiation with respect to r in the operator Ly _ |, can be 


expressed through the operator 4, i.e., Ly _ |, should be a polynomial of degree m/2 in A if m is even, and of degree 
(m + 4) if m is odd. 


When m = N the equation contains the function Ug only: 
L, (A) U,=0. 
It follows that the characteristic equation of the homogeneous system of equations will be of the form 


L,(p) = 0. 


It is clear that the number of roots of the characteristic equation, i,e., the eigenvalues x, is equal to N/2 if 
N is even, and (N + 1)/2 if N is odd. To each x’ there belongs a pair* of eigenfunctions Up, These eigenvalues and 
the corresponding solutions will be referred to as principal. 


*For each x”, the special solutions should be formed from AU = %7U, which in general has an infinite number of 
solutions. However, assuming separation of variables, one can reduce it to an ordinary second -order differential equation, 
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Next, let us proceed in the reverse order and consider the equation for m = N~-1: 
U1 = [(¥m) 


On solving this equation we obtain 


U,= Lj" (A) 4 UG, (20) 


where Uf is the solution of the homogeneous equation, i.e., L,(A)Uy = 0, Here we come across new eigenvalues 


which are the roots of the equation Ly(p) = 0. We shall call them subsidiary roots, There are two subsidiary proper 
solutions to each root, 


The function U, should satisfy Eq. (15), according to which the quantities (VV)U, should be expressable in 
terms of lower-number functions only, It follows that 


(VV) Uf =0. (21) 


Similarly, the solution for the functions U, (for n < N) is of the form 


subject to the condition 
(WV) U*=0. 


The functions Uy can be found directly from Eq, (16). 


It is clear that it is possible to determine the number of arbitrary coefficients in the general solution. We re- 
call that each function U, has 2n + 1 linear independent components, i.e,, each equation of order n is in reality a 
set of 2n + 1 equations, This means that in the solution for U,, to each subsidiary root there correspond .2n + 1 pairs 
of different subsidiary solutions. The condition given by Eq. (23) represents an equation of order n - 1, i.e., equi- 
valent to 2n—1 equations, and since it should be satisfied by each linearly independent solution of the pair, it follows 
that to each subsidiary root there correspond 2n~1 pairs of conditions, Consequently, the number of arbitrary coef - 
ficients belonging to a single subsidiary root is equal to 


2 [(2n 1) —(2n — 1)] = 4. 


Solutions corresponding to the principal roots do not have to obey additional conditions, although the number 
of spherical components in Up is equal to unity. It follows that to each principal root there correspond two arbitrary 
coefficients. It is now quite easy to compute the tutal number of arbitrary coefficients. 


Let us divide Eq. (17) into pairs beginning with the N—1 term. In each pair the degree of the operator-poly- 
nomial Ly _,,, with respect to A is the same and equal to t (the number of the pair), i.e., the number of subsidiary 
roots in the pair is 2t, The following two cases will be considered separately, 


1, N odd. After subtracting the equation of number N we obtain N equations, The number of pairs is equal to 
(N~-1)/2; and the zero-number equation does not form a pair, The number of all the arbitrary coefficients is 


2. N even, Let us divide the N equations, with the exception of equation of number N, into pairs, The total 
number of pairs is N/2, Of these, (N/2)—1 pairs contain subsidiary roots only. In the last pair one-half of the roots 
are subsidiary and the rest principal, The number of arbitrary coefficients is 


(19) 3 
(23) 
| 8 w+). 
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8 = NN +1). 


t=1 


We thus see that the number of arbitrary coefficients in the general solution of the system is independent of the 
parity of the approximation and is equal to N(N + 1). 


In order to obtain the solution of the nonhomogeneous system of equations it is sufficient to add to the general 
solution of the homogeneous system a special solution associated with the presence of the righthand side, which in 
our case can be expressed in terms of the polynomials Q,. Proceeding in accordance with the above scheme, we 


shall find that the righthand part qn—mlQo; Qi; ee Qnl, will appear in Eq. (17), and hence the required 
special solution will be of the form 


(A) qn [Qo; Q, Qn]. (24) 


A more detailed analysis of the solution is not of interest in the present context. 


Boundary conditions in the Py approximation 


It is usually assumed as self-evident in the method of spherical harmonics that all the spherical harmonics 


retained in the P,, approximation should be continuous across the separation boundary between the two scattering 
media, 


In one-dimensional problems, the number of arbitrary coefficients in the solution is equal to the number of 
spherical harmonics if the order of the approximation is chosen to be odd, In this case the requirement that all the 
harmonics should be continuous can be satisfied, and it is probably for this reason that one-dimensional problems are 
always discussed in the literature in terms of the odd (most frequently Py and Ps) approximations, However, in the 
general case, the continuity condition for all the harmonics, right up to N, cannot be satisfied even in the case of the 
odd approximations, owing to the lack of coefficients, This can be easily verified in the case of the Py approximation 
for the three-dimensional problem, The number of spherical harmonics in this case is equal to four, while the num- 
ber of coefficients given by the formula derived above is two, It followed that the conditions for the mating of all 

the harmonics cannot, in general, be mathematically incorrect, The problem of the boundary conditions in the 
general case is treated in Davison's book [3], but in our opinion a satisfactory solution of this problem is not given, 

We shall try to derive the boundary conditions which, without loss of generality, should remain applicable for any N, 


We note that in the derivation of the initial equations (8) or (13) it is not necessary for the cross-sections to 
be constant in the medium, and hence the equations are applicable to the heterogeneous ensemble as a whole, pro- 
vided it is assumed that the cross-sections are piecewise constant functions of the spatial coordinates, It is important 


that the equations are assumed to hold also at points lying on the boundary itself. This will enable us to formulate 
the necessary boundary conditions, 


In a medium with continuously distributed sources, the terms S,, and Ep Y, in Eq. (8) are bounded everywhere, 
including points on the boundary, It follows that the magnitude of div Fyjshould also be bounded at all points, We 
shall derive a condition for the magnitude of divF,, to remain bounded on the separation boundary of two media. 


Let us construct three mutually perpendicular unit vectors v, 1, € at an arbitrary point fp on the separation 
boundary. Suppose that one of them, for example, v, is in the direction of the normal to the surface, It is clear 


that the two remaining unit vectors and € will be tangential to the surface, In the neighborhood of the point fo, 
he vector F, can be represented by 


in which case 


OF... OF 


The components of the vector F,, in the two media are finite at all points adjacent to the boundary, This means that 
on the boundary itself they are either continuous or have a finite discontinuity, 


7 


Let us assume that Fyn and F,, are discontinuous. This means that they are indeterminate on the boundary, 

OF 

The derivatives 7g and -— " along the 7 and ¢ directions, which are tangential to the boundary, will in 
% 

general also have a finite discontinuity and will be indeterminate at tp However, since the discontinuity is finite, 

the indeterminancy does not mean that the function will be infinite, i.e,, the discontinuity in the components F pn 

and F_, does not contradict the fact that divFp(r9)is bounded, If, on the other hand, the component F py has a dis- 

continuity at fo, then in accordance with the definition of the derivative, 


Fay (to OV) # Fry — Ov). 


It follows thatdivF, can only be bounded at fp, if F,,,, is continuous at this point, Since the point fp is arbitrary, the 
result is general and holds at all points on the separation boundary. The boundary conditions for the P,; approximation 


can be expressed in terms of the requirement that the quantities Fpy should remain continuous for all numbers between 
0 and N. 


Let us determine the number of equations which have to satisfy this requirement, According to Eq, (12), for the 
Py approximation 


Fa = (VU ngs + (2n — 1) pU — (25) 
for n< WN; 


Fy = [(2N — 1) — 
The boundary conditions can then be written down in the following form: 
[(2N — 4) (vp) — p? (27) 

(W)U ngs + [(2n— 1) 

are continuous for n < N, 
The components of the normal v,, Vy» Vz are determined by the form of the boundary. 


The condition given by (27) is equivalent to 2N + 1 algebraic equations containing only 2N—1 components of 
the function Uy =": It follows that the condition (27) can only be satisfied if the function Uy _, is continuous. 


Consider now the conditions (28) for n = N-2, in which case 
4-[(2N —5) (vp) — p? (vV)] — 


is continuous, 


Since the continuity of U,)_, is a consequence of the continuity of (vV) U1, it follows that the second 
term in Eq. (29) should also be continuous, and this means that the function Un-3 should be continuous, Proceeding 
in this way we find that all the spherical harmonics Un ~-1- UN ~3, Un 5, etc., are necessarily continuous, We 
shall call these the first group of conditions, 


This conclusion does not hold for the harmonics Uy, Uy 9 Un ~4, etc. For these harmonics the boundary 
conditions have to be written down in the original form (28), i.e., 


ar... (to+5¥ )—Fuv( ) g 
ON 

since 

and 

(29) 
29 


(v¥) Uy + [(2N — 3) (vp) — (WV)] 


are continuous, 


It is easy to see that the neglect of the function Uy in this group of conditions would lead to the continuity of 
the functions Uy _ », UN -4, Un ~6, etc. It follows that in our approximation the discontinuity in the above func- 
tions is of the order of Uy. If with increasing n, the spherical harmonics U, become attenuated, then in the limit, 
i,e., for an infinitely increased N, the magnitude of the discontinuity would become arbitrarily small, Let us call 


these the second group of conditions, We can now calculate the number of equations to which the above conditions 
can be reduced, 


The mating of each harmonic of number N~2t~1 (in the first of the above two groups) is equivalent to a set 
of 2(N~2t—1) + 1 equations, Each mating in the second group is also equivalent to 2(N - 2t— 1) + 1 equations, 
In summing up all the equations it is necessary to consider the limits of t in the following two cases; 


1) N even, Here, in the first group 
1<2t+1<N—1 


and in the second group 


O<2t<N—2, 


i.e., in both groups 0 = t = (N/2)—1, It follows that the total number of all the matings is 


2 [2(N —2¢—1)4 14] =N(N +1). 


t= 


2, N odd. In the first group 


and in the second 


0<2t<N—1, 


i.e., in both groups 0 = t =(N- 1)/2, Hence the number of matings is 


N—1 


2 
2 2, [2(N — 2¢— 1) 4-1] =N(N +1). 


The coefficients of the two media take part in the mating process, although one-half of the coefficients in 
each medium are used to satisfy the conditions on the opposite boundary, at infinity, at the origin, etc. Thus, the 
total number of unknown coefficients and the number of mating equations are equal, whatever the order of the ap- 
proximation, In special cases when, owing to symmetry conditions, some of the harmonics vanish, the number of 
matings is correspondingly reduced, but then some of the subsidiary solutions will vanish. In one-dimensional prob- 
lems all the subsidiary solutions should vanish, since otherwise condition (23) cannot be satisfied, Consequently, the 
number of coefficients is equal to twice the number of the principal roots, i,e., 2N if N is even and 2(N + 1) if N is 
odd. In this case each spherical function contains only one component and hence, the number of mating equations 
is also equal to N and N + 1, respectively, as can be easily verified by direct calculations, It is clear that with odd N 


the number of equations is equal to the number of spherical harmonics, and hence it is necessary to mate all the 
harmonics separately, 


Since the spherical functions are orthogonal it may be shown that the boundary conditions derived above are 
equivalent to the requirements that the following integrals should be continuous : 


1<2t4-1<N, 
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\ (r, 2) Yum (2) dQ, 


4n 


O<m<n, 0<n<N. 


In other words, these boundary conditions ensure the continuity of the function QF with the aid of spherical moments 
of order up to N inclusive, 


It is well-known that the continuity of QyF can be obtained directly from the original transport equation, but 
in the exact solution the multiplier Q,, is unimportant because it is common, As we have seen, Qy plays an impor- 
tant role in the method of spherical harmonics. The requirement that F rather than Q),F should be continuous leads 
in the Py approximation to incorrect boundary conditions, since it involves the mating of all the spherical harmonics 
of order 0 to N, which in general is not possible, 


Extrapolated Boundary of a Semi-infinite Medium SUMMARY 


Order of ap- Boundary with | Boundary with We have obtained the boundary conditions only for 
proximation vacuum * black-body the case where the Py approximation of the method of 
spherical harmonics is used in both of the adjacent media, 
Py 0.667 0.577 If, on the other hand, the solution of the transport equation 
P2 0.667 0.776 is required on one side of the boundary only, then the num- 
Ps 0.705 0.694 ber of free coefficients available for the conditions on the 
P, 0.702 0.730 boundary is reduced by a factor of 2, In that case, the 
Peo 0.710 0.710 boundary condition should be specified by a distribution 
With Marchal conditions function eatined within one-half of the total solid angle 
(if the medium is convex throughout), These quite new 


considerations prevent a straightforward extension of the 

above boundary conditions to the present case, However, 
one can probably maintain that any approximate boundary conditions must be based on the behavior of the function 
QyF and not F alone, For example, on the boundary with a vacuum, and in the case of the Marshak conditions [3,5], 
this consideration enables one to decide in favor of functions which are odd in Qy, since the expansion of F in the 
interval — 1 = Qy = 0 in terms of odd harmonics is roughly equivalent’ to the expansion of QyF in terms of even har- 
monics, while the expansion of F in terms of even harmonics cannot be reduced to the expansion of Q,F. The num- 
ber of odd functions of order between 0 and N is $N(N + 1), i.e., the number of boundary conditions is the same as 
the number of free coefficients. 


The vacuum is physically equivalent to a perfect black body, and hence, it is occasionally replaced by a 
medium in which the capture cross-section tends to infinity, This formal device enables one to use boundary con- 
ditions of the type discussed in the present paper even though, strictly speaking, the Py approximations of not too 
high an order are not applicable in a strongly absorbing body. 


The table gives a summary of the computed extrapolated boundary for a semi-infinite non-absorbing medium 
with a plane surface, It is known that with an isotropically scattering medium the extrapolated boundary lies at a 


distance of YA, from the real boundary, The table gives the values of Y in the first four approximations, and also 
the exact value. 


It should be pointed out that in the case of one-dimensional plane geometry, and in the absence of absorption, 
the contributions due to even approximations are the least significant. The above table can only provide limited 
information about the errors and the convergence of even approximations as compared with odd approximations, 


In multidimensional problems, concerned with the diffusion of neutrons in directions perpendicular and tangen- 
tial to the boundary, even approximations will have the property that all the components of the diffusion current will 
be continuous and not only the normal component as in the odd approximations, This is important in the averaging 
of the diffusion coefficients over various directions in a heterogeneous medium, The P2 approximation, which has 
so far remained practically unexplored, is of interest, 


The Pz approximation is in many cases more accurate than the Py approximation, and is capable of including 
certain effects which are omitted in the Py approximation, for example, the dependence of Y on absorption in the 
medium, the curvature of the surface, and other defects, The Pz; approximation gives a more correct description of 
the diffusion of neutrons in multidimensional problems, It gives good results in media with relatively large absorp- 


i 
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tion, At the same time, the P2 approximation involves roughly the same order of complexity as the Py approxima- 
tion, and is of course, much simpler than the Ps approximation, However, as in the case of the Py approximation, it 
can only provide the asymptotic solution of the transport equation, This must be borne in mind whenever the di- 


mensions of the medium are small, 


In conclusion, the author wishes to thank Candidate of Physicomathematical Sciences, V. V. Orlov, for his 
consideration, interest and frequent discussions of the results, Thanks are also due to E, I, Gladyshev for assistance, 


and to all other colleagues who took part in discussions, 
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THE RED COLORATION OF MINERALS IN URANIFEROUS VEINS 


Yu.M.Dymkov andB. V. Brodin 


Translated from Atomnaya Energiya, Vol. 10, No. 1, 
pp. 35-42, January, 1961 
Original article submitted June 15, 1960 


This paper describes different types of hematitization of ferruginous carbonates (siderite, ankerite 
and ferruginous manganocalcite) in uraninite -bearing veins. 


It has been ascertained that hematitization not spatially associated with uraninite may be pre-ore 
or post-ore and may differ in origin. It may be related to oxidation of iron carbonates as a result 
of chemical activity by hydrothermal solutions and may also be due to mechanical transfer of fine 
particles, by solutions, of previously formed hematite. Hematitization, systematically associated 


with uraninite and forming, in ferruginous carbonates, aureoles about uraninite lenses, is of radio- 
genic origin. 


In conclusion, note is taken of the practical value of the various types of hematitization as indica - 
tors in prospecting for uranium in hydrothermal vein deposits. 


Introduction 


For a long time miners and geologists have known that a reddish aureole ("red alteration”) is commonly found 
around uraninite in carbonates, barite, and feldspars, the color being due to fine flakes of hematite or finely-dispersed 
hematite in these minerals, Such hematitization is thought to be the result of radioactive irradiation [1,2] or of oxi- 
dation -reduction processes originating in the solutions between Fe” ions and uranyl ions (UO; ) [3,4]. It has been 


suggested that the hematitization may be associated with the oxidizing action of carbonic acid solutions devoid of 
sulfur ions [5]. Not long ago, a new view was advanced — that the cause of the red coloration in rocks containing 
uranium deposits is due to the oxidation of iron associated with ".., radiation-chemical reactions originating in hydro- 
thermal uranium -bearing solutions" [6]. It should be noted that the existing views give no consideration to the possi - 
bility of different origins for the hematitization found in uraninite-bearing veins. It is well known that hematitiza- 
tion is also observed at nonuraniferous deposits and in barren veins at uranium deposits. On the other hand, some dis- 
tricts containing rich uranium ores show no red alteration [7,8]. 


The authors of the present paper have studiedthe red coloration of minerals in uraniferous carbonate veins con- 
sisting of the paragenetic sequence of siderite, ankerite and manganocalcite and of calcite of several generations. 
The uraninite in these veins is closely associated with late nonferruginous calcite and a very small quantity of sulfides, 
most of which formed after deposition of the great bulk of the uraninite, 


In spatial relations the red coloration in the uraniferous carbonate veins may be divided into two types of hema- 
titization, varying in practical importance: hematitization spatially unrelated to the accumulations of uraninite and 
hematitization forming reddish aureoles about uraninite lenses, 


Hematitization Spatially Unrelated to Uraninite 


In relation to the time of deposition of the uraninite, hematitization may be divided into an early (pre-ore) 
and a late (post-ore) phase in the vein carbonates. The best defined is the early phase, being characterized by an 
intimate association of hematite and goethite. Both minerals arise from oxidation of ferrous iron in siderite, ankerite, 
and ferruginous calcite. The hematite-goethite association is commonly observed through extensive parts of the veins, 
measured along the strike or down the dip in tens and hundreds of meters, and occurring in places where there is ab- 
solutely no uraninite nor any other minerals of the uraninite stage. In the lower parts of the veins the hematitization 
locally increases in intensity, whereas the level of mineralization in the same vein is considerably nearer the surface. 


: 
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Fig. 1. Intersection of uraninite veinlets (1) with 

fragments of extensively hematitized ankerite (2) Fig. 2. Mineral level in a solution cavity of a carbonate 

and hematitized ferruginous manganocalcite (3). vein. The level consists of alternating layers of "alluvial" 

The hematitization occurred before the deposition hematite, clay material, and rock fragme nts cemented 

of grayish white granular calcite (4). Drawn from with calcite; 1) hematite, partly invading along growth 

a hand specimen. zones in crystals of manganocalcite; 2) colloform man- 
ganocalcite; 3) radiating manganocalcite; 4) partly hema- 


iti ‘ li and specimen. 
Of the iron oxides, the most widespread is dis- 


persed cherry-red hematite (identification confirmed 

by x-ray studies); well-crystallized hematite and radiating goethite are less abundant. The latter occurs as metacrysts 
in already hematitized siderite, partly replacing early hematite. Locally, the goethite forms extensive accumulations 

of radiating aggregates. It is characteristic that in such zones chalcopyrite, of the early sulfides, proves to be unstable, 
being replaced by bornite and then by cubic chalcocite. 


Thus, the early hematitization in the investigated uraniferous veins does not differ fundamentally from the he- 
matitization observed in a number of nonuraniferous leadezinc deposits [9-11]. 


The deposition of the uraninite occurred considerably after the early hematitization, and the formation of 
goethite was not genetically related to either process. Proof of this is found in widely observed intersections of hema- 
titized siderite and ankerite with fresh manganocalcite (containing about 1% FeCOs), its pale rose or white color un- 
changed. In individual veins fragments of hematitized ankerite and early manganocalcite are cemented by gray 
granular calcite and then cut by veinlets of uraninite (Fig. 1), Furthermore, the early hematitization in ore veins 
apparently had an adverse influence on the localization of uraninite. In zones where siderite has been extensively 
replaced by goethite and early hematite, uraninite is generally absent or present only in small quantities. 


Apart from the described early hematitization, the age of which is clearly established, red coloration of dis- 
persed hematite and later calcite is also observed locally. This hematite probably represents a residuum from the 
solution of the early hematitization of the carbonates. Dispersed and colloidal particles of hematite invaded crystal 
aggregates of calcite along the seams between grains, and also attacked calcite in the process of crystallization. On 
the floor of hypogene solution cavities of early carbonates (ankerite and early manganocalcite), one may find layered 
aggregates of dispersed hematite* (Fig. 2). The hematite penetrates along the contacts between grains of the ferrugi- 
nous manganocalcite to some distance from the layered precipitates, although the great bulk of the manganocalcite 
is not hematitized. These facts represent good examples of "mechanical" redeposition. 


It is possible that some of this hematite was formed by oxidation of ferrous iron in the carbonates in the solu- 
tions from which the minerals of the uraninite stage were deposited [13]. It is assumed that such solutions, during the 
transport of the uranium, were essentially oxidizing. The activity of these solutions was also responsible for dolomiti- 
zation of manganocalcite, producing dark gray borders about the grains of this mineral, and for the precipitation of 
hematite immediately before the deposition of uraninite. 


The late hematitization in the ore veins, occurring much later than the deposition of uraninite and the succeed- 
ing sulfides (Fig. 3). is of limited development. Apart from zones of red coloration in post-ore calcite, microscopic 


*Such formations have been termed "mineral levels" [12]. 
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inclusions of goethite or hematite spherulites are found in individual crystals of calcite and quartz. The origin of 
the late hematitization is unclear, 


Hematitization Associated with Uraninite (Radiogenic Hematitization) 


Hematitization systematically associated with uraninite is found only about uraninite lenses, forming distinct 
aureoles (Fig. 4). In carbonate veins the hematitization involves only the ferruginous carbonates. Nonferruginous 
carbonates, even near large accumulations of uraninite, remain unaltered. 


Near the ore lenses ankerite (9.27% FeO) shows various 
degrees of hematitization; manganocalcite (1.08% FeO) is 
very strongly hematitized, and white calcite (up to. 0.28% 
FeO) only weakly hematitized, The latter mineral commonly 
remains unchanged. 


The extent of the aureoles of hematitization along the 
strike and down the dip corresponds approximately to the di- 
mensions of the uraninite lens or nest, and the width is gener- 
ally limited by the thickness of the vein.” Below is shown 
the relation between dimensions of aureoles of reddened car- 
bonates and the dimensions of uraninite segregations. This 
Scale £9 ~~ relation is most clearly observed in small point-like segrega- 
tions: 
Fig. 3. Intersection of uraninite veinlets with 
later hematitized calcite: 1) grayish-white Diameter of uraninite © Width of hematitization aureole 


calcite; 2) uraninite; 3) rose-colored granular spherulites, mm in the carbonates, mm 
calcite; 4) late hematitized calcite; 5) hema- 0.05 2 ; 


titized ankerite; 6) chalcocite with inclusions 0.1 

of chalcopyrite grains. Drawn from a polished 0.3 

hand specimen. 0.5 25 
0.15- 1(incrustation) 1-8 
3.5 more than 180 


The morphology of the red-colored aureoles depends on the internal structure, the dimensions, and the position 
of the ferruginous carbonates in the veins. 


Manganocalcite shows the greatest alteration near the uraninite. The primary color of the manganocalcite is 
rose, pale rose, or white. In individual zones not associated with segregations of uraninite, the manganocalcite is 
zonally or spottily (near minute fractures) gray or dark gray, apparently because of disseminated manganese oxides. 


TABLE 1. Chemical Composition of Manganocalcite (analyzed by B. M. Eloev) 


Insolu- 
Color MnO Fe203 2 ble resi- | Totals 


analysis 


dues 


Red=brown, beige 5.3 49.55 0.37 43.10 0,20 99.85 
Dark gray 4.94 49.85 de- 43.72 0.20 100,36 

tected 
Rose 5.05 6 49.30 0.48 43.30 0.22 100.11 
Cherry -red 6 4, ; 50,02 1.20 42.67 0.20 99.76 
Cherry -red j 6 2.63 42.35 0.20 100.33 


N= 


a) 


Hematitization near uraninite rarely extends to the zones of gray coloration, Individual zones of earth manganocal- 
cite are grayish brown or light brown with rose-yellow tints. Admixtures of goethite generally impart such colors to 


*In places such aureoles of red coloration were observed superimposed on zones of early hematitization, creating a 
very complex picture of alteration in the vein carbonates. 
+ An inner dark-colored ring 0.75 mm wide can be distinguished in this aureole. 
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the carbonates. Around individual segregations of uraninite, locally in zones, one may see dark red, beige, and dark 
gray colors. The mangocalcite varies somewhat in composition in the different zones, as chemical analyses show. 


Thus, gray zones contain a high quantity of magnesium and are characterized by a rather low content of ferrous 


oxide (see analysis 2, Table 1)° 


Fig. 4, An aureole of red coloration about a nest of uraninite; 1) hematitized siderite (without 
goethite); 2) smoky quartz; 3) hematitized ankerite; 4, 5, 6) different generations of mangano- 
calcite; 5a, 6a) intensely hematitized manganocalcite; 7) columnar calcite; 7a) columnar calcite 
hematitized; 8) scalenohedral white calcite (nonferruginous); 9) uraninite; 10) scalenohedral 
late calcite; 11) stubby rhombohedral calcite, Center of a drusy cavity. Drawn from a zone in 


the roof of a mine working, 


A comparison of hematitized cherry-red (0.5 m from uraninite) and rose-colored (3 m from uraninite) mangano- 
calcite from the base of an earthy incrustation shows a connection between the hematitization and the oxidation of 


ferrous iron that was present at the site (see analyses 3 and 4, Table 1)."" 


A recomputation of the data of analyses 3 and 4 (Table 1) to metals indicates that the quantity of iron increases 
0.14% in the hematitized zones (analysis 3, Table 2). In zones where uraninite is segregated in fractures that cut 
manganocalcite, films of introduced hematite are observed along individual slippage plains, barren of uraninite; 
these are due to a substantial increase in the Fe,O, content in the manganocalcite (see analysis 5, Table 1). 


Partiai analyses of unaltered and hematitized carbonates (see analyses 1 and 2, Table 2) show that the amount 
of iron may increase substantially both in the aureoles of hematitization of manganocalcite and in the zones of early 
red coloration of ferruginous calcite. The increase in iron in the colored aureoles about uraninite lenses attests not 
only to the oxidation of ferrous iron of carbonates in place but also to some introduction of iron from the solutions 
that deposited hematite at the zones of oxidation. 


*Spectral and x-ray analyses of the ferruginous manganocalcite have shown that the content of microscopic 
admixtures did not change noticeably during hematitization. 


““No “blackening” of manganocalcites (by oxidation of Mn* ) was observed in the aureoles of hematitization, 
over Fe”, This may be due to the relatively high oxidation-reduction 


potential of the reaction + 2e (Mn* gives a rose color) in comparison with Fe“ + e [14]. 
Intense oxidation of Mn®’ in manganocalcite is known only in zones of exogene processes. 


despite a marked preponderance of Mn 


207 
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X-ray studies of the insoluble residue and the elutriation products of comminuted manganocalcite have shown 
that the great bulk of the staining substance, as also in the zones of early coloration, is hematite.” Brown oxides of 
iron (goethite) are less commonly present near the uraninite lenses than is hematite. In the aureoles of hematitiza- 
tion they formed mostly after the hematite, being related in time with the late (post-ore) hematitization. 


There is no reliable data indicating simultaneous formation of hematite and uraninite. Where nonferruginous 
pre-ore calcite is associated with uraninite, no aureoles of red coloration are observed about the uraninite. Hema- 


tite and goethite are also absent in the white nonferruginous scalenohedral calcite that was deposited simultaneously 
with uraninite in thin veinlets. 


TABLE 2. Change in Content of Iron Oxides during Hematitization of Pre-Ore Carbonates 


Number 


of Minerals 


Content in 


unaltered mineral, % 


Content in 


Quantity of 


hematitized mineral, %|FeO changed 


analysis 


FeO 


Fe,O, 


FeO 


Fe,0s 


to Fe,O, 


Total intro- 
duction of 
iron 


Type of 
hematiti- 
zation 


1 | Early columnar 
calcite 


1.58 


0.50 


2.90 


0.72 


1.12 


(*) 


Colloform 


manganocalcite| 1,00 0.12 | 0.79 


3 |The same 1.08 0.48 | 0.61 


(*) Early, not associated with the formation of uraninite. 
(**) Aureole of alteration about accumulations of uraninite. 


These facts indicate that the aureoles of red coloration in the carbonates about uraninite are not associated 
directly with the deposition of the uraninite, but were formed later. The hematite in these aureoles, as analyses 
have shown (see Tables 1 and 2) was formed not only by oxidation of ferrous iron in the carbonates but also by oxi- 


dation of iron introduced from outside in later stages. Different relative ages for the iron oxides attest that the 
oxidation of Fe” about the uraninite continued for a long period of time. 


The morphological features of the segregated hematite in the zones of red coloration in the carbonates and, 
above all, the isolation of hematite along the contacts of crystal grains of manganocalcite (Fig. 5) show that the 
hematite was deposited from hydrothermal pore solutions (the calcite rock containing the lenses and nests of uranin- 
ite has a porosity of 0.7-2%). 


The striking restriction of hematite to uraninite segregations, the later formation of the hematite, and the ab- 
sence of any connection between hematitization and the deposition of uraninite—all these point to a relationship 
between the formation of hematite and the effect of uraninite on the hydrothermal solutions. An increase in the 
width of an aureole of hematitization up to a certain limit (on the order of one meter) corresponding to an increase 
in the mass of uraninite in a lens suggests, in some degree, a dependence of the penetration of gamma radiation in 
the rocks on the mass of ore. This leads one to explain the phenomenon of the aureoles of hematitization by radia - 
tion chemical processes.** It may be assumed that the aureoles of red coloration were formed by oxidation of a fer- 
rous iron by pore and capillary solutions, the oxygen in which came from radiogenic-chemical reactions associated 
with prolonged radiation of uraninite. The extensive size of the aureoles of red coloration may be explained by the 
action of penetrating gamma radiation.*** Prolonged action of radon on the pore solutions that produced the expan- 
sion of the aureole of red coloration may have been of fundamental importance. 


* The x-ray analysis was made by N. G. Nazarenko. 

** However, if it is considered that the oxidation occurs because of alpha and beta radiation of uranium ions (or col- 
loidal particles) [6], then it is necessary to assume that these ions existed in solution an extremely long time and in 
very great concentrations. On the other hand, during such a process (even if the intensity of radiation were sufficient 
for development of hematite) uraninite should not be associated with the hematite, since the uranium would be oxi- 
dized by the same oxygen and would be found in the hexavalent state. 

***It is assumed that the inner aureole (0.75 mm) may be referred to the aureole of beta radiation. The aureole of 
alpha radiation may be observed only under the microscope. 


. 
0.21 1.40 (**) 
1.20 0.47 0.14 (**) 
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As a result of decay of the uraninite, there may have occurred displacement and redeposition of the uranium in 
a disseminated state in the carbonates and other minerals. This may have led to the development of numerous new 
centers of red coloration farther removed from the uraninite lenses. 


Existing data on the content of uranium in the hema- 
titized and nonhematitized carbonates do not permit us, 
however, to claim a connection between intensity of hema- 
titization and the presence of disseminated uranium in a 
mineral, 


The clearest increase in uranium content is found in 
hematitized (of the early coloration) siderite (0.00027%, 
0.0045% U), where the quantity of uranium increases a 
hundredfold in comparison with the unaltered (0.000052%) 
and tenfold in comparison with the slightly dolomitized 
siderite (0.00016% U). Inasmuch as the uranium mineraliza- 
tion occurred after the early red coloration, the increase in 
uranium content in the hematitized carbonates is apparently 
a consequence of alteration of siderite associated with the 
hematitization, and is not a cause of the hematitization. 


In the hematitized ankerite and manganocalcite such 
increase in uranium content is not observed. In one of the 
two analyzed samples of hematitized (early coloration) 
ankerite, the content of uranium shows a tenfold increase 
(0.001%U); no change is detected in the other (0.00017% U). 
The intensely hematitized manganocalcite from the aureoles 
(0, 00014-0.00052% U)shows no substantial increase in urani- 
um content over the slightly altered (0.00001-0.0001% U) 
varieties (in one sample the increase was on the order of 
0.00105% U). The nonferruginous (stained with later iron 
oxides) yellowish-grayish-white calcite accompanying 
uraninite generally contains 0.00n% uranium. 


Analyses of ten samples of later post-ore nonferrugi- 


Fig. 5, Distribution of hermatite (black) in hema- nous calcite show a content of uranium on the order of 
titized colloform manganocalcite from an aureole of —_-0.000n-0.0000n%, A comparison of the intensely hematitized 
red coloration about uraninite, Photomicrograph of (hydrothermally altered) and nonhematitized rocks alsoshows 
a thin section. Parallel nicols, X 170. no substantial change. Thus, no clear relationship between 


hematitization and disseminated uranium has been detected. 


C ONCLUSIONS 


1, Several types of hematitization of vein carbonates are discussed; these are distinguished both in time of 
formation (relative to the deposition of uraninite) and in localization in space; 

a) hematitization unrelated in space or time with uraninite segregations (pre-ore and post -ore); 

b) hematitization forming aureoles about accumulations of uraninite. 


2. Three genetic types of hematitization are found in veins; 

a) early hematitization of iron-bearing carbonates, locally accompanied by segregations of goethite. This 
hematitization is related to oxidizing conditions, normal for the end stages of mineralization in sulfide deposits; 

b) "alluvial" hematitization, appearing repeatedly at various stages of development of the carbonate-urani- 
nite veins. It is associated with solution of hematitized carbonates and mechanical transfer of colloidal and dispersed 
particles of iron oxide. Solutions containing suspended particles of hematite invaded along contacts between grains 
and through fractures of the carbonates, staining them irregularly; 


c) hematitization forming aureoles about lenses of uraninite, resulting from radiation-chemical processes. 
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3, Extensive hematitization of various types attests to the repeated manifestation of an oxidation environment 
under hydrothermal conditions, favorable for the transfer of uranium, The formation of uraninite was not accom- 
plished by hematitization of the earlier iron-bearing carbonates. 


4. Early (pre-ore) hematitization is not a direct indicator of uranium. The “alluvial” hematitization, appear- 
ing at various stages throughout the course of vein development, merely attests to hypogene solution of early hemati- 
tized carbonates, "Radiogenic” hematitization may indicate large accumulations of uraninite near it. 


The identification of hematitization with one of the indicated three genetic types may be made by analysis of 
the age relations of the altered or unaltered iron-bearing carbonates at a given part of the vein. 
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THE CONNECTION BETWEEN THE STRUCTURE AND ANISTROPY 
OF THE THERMAL EXPANSION OF URANIUM, 
AND PLUTONIUM 


NEPTUNIUM 


N. T. Chebotarev 


Translated from Atomnaya Energiya, Vol. 10, No. 1, 
pp. 43-49, January, 1961 
Original article submitted April 11, 1960 


A connection has been established between the anisotropy of the thermal expansion and the changes 
occurring in the structure of a number of modifications of uranium, neptunium and plutonium on 

heating. It has been shown that in all structures considered the anisotropy of the thermal expansion 
is due to a weakening of the four covalent bonds with increase in temperature and their nonuniform 
distribution in the various crystallographic directions of the lattice. The regularities given lead to 


a better understanding of certain specific properties of uranium, neptunium and plutonium which 
distinguish them from ordinary metals. 


Introduction 


IK is known that noncubic modifications of uranium, neptunium and plutonium have different values for the co- 
efficient of thermal expansion in different crystallographic directions, The anisotropy of the thermal expansion has 
an important effect on the properties of these metals, The different values of the coefficient of thermal expansion in 
the various crystallographic directions indicate that during heating or cooling, the crystalline lattice does not remain 

unchanged, This in its turn means that the mutual positions of the atoms in the crystalline lattice change. Therefore, 
for an analysis of this phenomenon it is essential to establish those relationships governing the change in the intera- 

tomic distances in the crystalline lattices with change in temperature. 


Structure and Thermal Expansion of &®-Uranium 


The alpha-phase of uranium has a rhombic structure with lattice parameters a = 2.853 A, b = 5.866 A and 
c = 4,955 A. The four atoms of uranium in the unit cell occupy positions 4(c) of the space group ])}7-—Cmcem: 


+ (000; 


where y = 0.105. 


On the projection of three unit cells on the (001) plane (Fig. 1) it can be seen that the atoms of zigzag chains 
threading the lattice in the direction of the c-axis are arranged in the form of layers with a distance between them of 
c/2. One of these layers is represented by the atoms ABCDEFO, the other by the atoms AyB)C,D,EyFyO,. On 
considering a tetrahedron with a base in the form of an isosceles triangle OAF and apex Oy, the following character 
ot coordination appears; 2d;= OO, = 2.77 A; 2 dp = OB = 2.85 A; 4 dg = OA = 3.26A; 4d, = OC, = 3.31 A. The mean 
interatomic distance d,, = 3.13 A. If the total coordination number is equal to 12 then the four shortest bonds 
are separated out with interatomic distances 2,77 and 2.85 A. 


If the triangle OAF was equilateral and the projection of the atom O, was exactly at its center, then the struc- 
ture of the a-uranium would be an ideal hexagonal packing. Table 1 compares the values of lattice parameters of 
such an ideal packing (with the condition r = dp, /2) with actual values of parameters for the structure of a-uranium, 
Here we give the coefficients of thermal expansion in the range 20-500°C according to the data of [1]. 


1 
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TABLE 1. A Comparison of the Lattice Parameters of «-Uranium and an Ideal Compact Packing 


Lattice Ideal values of Actual values of Extent of Coefficient of thermal ex- 
parameters parameters, A parameters, A deviation, % pansion along the axis, 10 


3.13 2.85 —9.0 32.9 
5.42 5.87 +8.3 -6.3 


5.11 4.96 -2.9 27.6 
0.167 0.105 -37 ~ 


It can be seen from these data that under the action of the above four covalent bonds the lattice suffered con- 
siderable deformation, being stretched along the b-axis and shortened along the a-and c-axes. 


The interatomic distances in the unit cell are determined from the following expressions: 
1 
V (4yb +e ; 
d,=a; 
b?: 


dy=5 Vai dy). 


Using expressions (1) we can plot the dependence of the interatomic distances on the temperature since the 
temperature dependence of the lattice parameters a, b and c are now well known [2], 


The results of the calculations are shown in Fig. 2. This 
figure also gives values of the parameter y which (taking into 
account the absence of noticeable changes with temperature in 
the difference in interatomic distances d4-d,) can be determined 
from the following expression: 


y= (2) 


With increase in temperature the short bonds (see Fig. 2,a) in- 
crease considerably, whereas the long bonds (see Fig. 2, b) donot 
change noticeably. The value of y also increases from 0,105 at 
room temperature to 0.112 at 640°C, which is in good agreement 
with the experimental data given in [1]. 


From equations (1) it follows that a = dy; b = |/ 4d? — d?; 


Fig. 1. The projection of three unit cells of 
a-uranium on the (001) plane. Height in 9 = 
fractions of a parameter c; @) 1/4; 0)3/4 


at a 


i.e.,that an increase in the short bonds d, and d, with temperature in the absence of noticeable changes in dg and dy 
leads to an increase in the parameters a and c and a decrease in the parameter b. 


The effect of short bonds on the anisotropy of thermal expansion can be seen from a consideration of their ori- 


entation with respect to the axes of the lattice. The values of the direction cosines for these types of bonds are as 
follows: 


Type of bond COS Gq cos Gp, COS My 


0 0.44 0.90 
2d,=2.85A 0 0 


b 
y 
: a 
8. 8 
| 
4. 
: 
a 
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It can be seen from these data that short bonds act mainly in the direction of the a-and c-axes, causing an increase 
in the coefficient of thermal expansion in these directions with respect to its mean value. On the other hand, in the 
diection of the b-axis, the coefficient of thermal expansion is much less than its mean value. 
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Fig. 2, Temperature dependence of interatomic dis- 
tances in the structure of &-uranium and the param- 
eter y: a) short bond and parameter y; b) long bonds, 


Fig. 3. Projection of four unit cells of a-neptunium on 
the (001) plane 
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This analysis shows that the strong anisotropy of the 
thermal expansion of a-uranium is caused by the consider- 
able weakening of the four covalent bonds with increase in 
temperature and their uneven distribution along the axes of 
the crystalline lattice. 


Structure and Thermal Expansion of 
a-Neptunium 


According to the data of [3], a -neptunium has a 
rhombic lattice, the temperature dependence of whose 
dimensions is given in Table 2. The structure corresponds 
to the space group piS -Pmen. In each unit cell there 
are eight atoms occupying the following positions: 


1 
4Npr — 4(c): + 
where y, = 0.208 + 0.006; z, = 0.036 + 0.006: 


1 
where y, = 0.842 + 0.006; z,=0.319 + 0.006. 


Figure 3 shows the projections of four unit cells of 
a-neptunium on the (001) plane, From a consideration of 
the diagonal cells ABCD, which consist of eight cells of the 
type EFGH and FIKG with dimensions a' = 3,399 A, c' = 
3.382A (c'/a* = 0.98) and an apex angle equal to 88° (the 
cells are on two levels) it can be seen that the structure of 
a-neptunium is the result of deformations of a body-cen- 
tered cubic lattice. 


Like the previously considered structure, in the a- 
neptunium structure four covalent bonds appear clearly 
with distances from 2.60 to 2.63 A. The mean interatomic 
distance is dp, = 3.10 A. 


In accordance with the data given in Table 2, the 
values of the coefficients of thermal expansion for a- 
neptunium in the range 20-275°C are equal to: 


a, = 24-10°* 1/°C 
a, = 25-10" 1/°C }, a, = 1/°C. 
a, = 1/°C 


It can be seen from these data that the anistropy of ther- 
mal expansion for a-neptunium is expressed rather weakly. 
Since the anistropy is connected with nonuniform distribu - 
tion of short bonds along the axes of the lattice, as was ob- 
served for the structure of a-uranium, we will consider the 
projections of these bonds for the a-neptunium structure 
(Table 3). 


It can be seen from these data that the covalent 
bonds are relatively uniformly distributed along all three 


axes of the lattice, which is in good agreement with the 
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absence of clearly expressed anisotropy of the thermal expansion. 


At a temperature of 278 45°C, a tetragonal structure of 8 -neptunium forms [3], the result of deformation of the 
body-centered cubic lattice, but with a degree of distortion much less than for a-neptunium. Of the eight bonds in 
this structure, four correspond to a minimum interatomic distance 4d =2.79 A and the other four to 4dmax = + 3.26A. 
Hence it can be seen that in the § -neptunium structure, four covalent bonds appear, although in correspondence with 
the increase in temperature the difference d_, - dyjjn here becomes much less than in the a-neptunium structure 
(0.23 and 0.48 A respectively). At temperatures of 550-570°C a body-centered cubic lattic of y -neptunium forms 
with eight identical bonds at a distance 2.97 A. 


Structure and Thermal Expansion of y-Plutonium 


According to the data of [4] y-plutonium has a rhombic lattice, the temperature dependence of whose dimen- 
sions is given in Table 4. 


TABLE 2, The Temperature Dependence of the Lattice TABLE 4. The Temperature Dependence of Lattice 
Parameters of a-Neptunium Parameters of y -Plutonium 


Lattice parameters, A 


Tem Lattice parameters, A 
Temp., °C a | b c 
a b c 


213 3.16052 5.76275 10.1442 
20 4.72340.001 | 4,887-+0.002 | 6.663--0.003 233 3.15909 5.76769 10.1615 
212 4.7460.002 | 4.909-+-0.003 | 6.704--0.004 258 3.15622 5.77371 10.1834 
275 4.7524-0.002 | 4.920-0.003 | 6.722-0.004 312 3.15397 5.78574 10.2290 


TABLE 3. Projections of Short Bonds on the Axes of the TABLE 5. The Temperature Dependence of the Inter- 


a-Neptunium Lattice atomic Distances in the Structure of y -Piutonium 
Projections on the axes of — 
the lattice, A é Interatomic distances, A 
Type of bond - — Temp., °C 
dg | dp d, ad | | 
1d,=2,598 A 0 1.79 1.87 | 213 3.022 3.160 3.286 
{dy—2.632 A 0 2.20 | 1.55 = 
2ds=2.634 A 4.72 | 1.30 | 1.94 312 33,052 3.154 3.295 
= 4.72 5.29 5.36 


The structure corresponds to the space group D3, —Fddd. 
position 8 (a): 


4.8: 
144;7+( 000;05>5;: 505; 550). 
wey ( ) 


The structure of y-plutonium is an alternation of several distorted pseudohexagonal layers placed perpendicu- 


lar to the c-axis. Figure 4 shows the projection of three unit cells on the (001) plane and one unit cell on the (100) 
plane. 


The interatomic distances in the y -plutonium structure and their temperature are given in Table 5, The mean 
interatomic distance at 213°C is dp, = 3.16A. 


From these data it can be seen that, similar to the previously considered structures, four covalent bonds appear 
clearly in the y-plutonium structure, the interatomic distances in the direction of which increase considerably (from 
3.022 to 3.052 A) with increase in temperature. The interatomic distances in the direction of the long bond change 
much less, the distance d, even being reduced somewhat with increase in temperature, 


| 

| 
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The values of the direction cosines for all three types of bond are as follows: 


Type of bond cos 


4d,=3.022A 0,26 0.48 
2d,=3.160A 1 0 
4dy=3 08 | 08 


In accordance with the data given in Table 4 the coefficients of thermal expansion of y-plutonium in the 
range 213-312°C are equal to: 


a, = ) 
a, = 39.5-10°° 1/°C , am = 34.5-10°8 1/°C. 
a, = 84.3-10°8 1/°C j 


On comparing these values with the direction cosines given above, it can be seen that the strong anisotropy of 
the thermal expansion is due to the short bond 4d,, mainly acting in the direction of the c-axis. In accordance with 
this, the coefficient of thermal expansion in the direction of the c-axis has a maximum value, in the direction of 
the b-axis it is much less,and in the direction of the a-axis it has a negative value. 


Fig. 4. Projections of three unit cells of y-plutonium on the (001) plane 
(a) and one unit cell on the (100) plane (b). 
Height in fractions of a parameter c: @) 0; ©) 1/4; @)1/2; ©) 3/4. 


Structure and Thermal Expansion of a-Plutonium 


According to the data of [5], a-plutonium has a monoclinic lattice, whose parameters at 21°C are: 


a = 6.1835 + 0.0005 A; b = 4.8244 + 0.0005 A; 
c= 10.973 + 0.001 A; B= 101.81 + 0.02°. 


The structure corresponds to the space group C3, — P2,/m. In the unit cell there are 16 atoms, occupying 
eight double positions: 
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Figure 5 shows the projection of four unit cells of a-plutonium on the (010) plane. The numbers near the at- 
oms indicate the number of the position, the atoms with a minus sign are placed at a height equal to 


2 (e): 


The approximate values of the x and z parameters are as follows; 


Parameters 


Puy 


0,767 
0.169 


Parameters 


Puy; 


0.459 
0.642 


parameter with respect to the remaining atoms 
lying in the plane of the diagram. The structure 
of the &-phase can be represented as an alterna- 
tion of very strongly distorted pseudohexagonal 
layers placed perpendicular to the b-axis at a 


distance equal to a half of the parameter. 


Fig. 5. Projection of four unit cells of plutonium 


on the (010) plane, 


TABLE 6, 


Structure 


half of the 


Interatomic Distances in the &-Plutonium 


Fosition of 


atoms 


Short bonds 


Long bonds 


interatomic 
distances 


number of | 


interatomic 
distances 


Coordination 


~ 


< ~ 
si<< 


52— 


2. 
2. 
2. 
2.5. 
2. 
2. 
2 
2. 


58—2.73 
52—2.70 
68—2.70 


3.40—3.51 
3.33—3 .62 
3 20—3.67 
3.20—3.54 
3.30—3.56 
3.33—3 .62 
3.34—3 .56 
3.30—3.67 


TABLE 7, Projections of d-dyp on the Main Crystallo- 


graphic Directions of the &-Plutonium Lattice 


Direction 


Projections of d-d_ 


for short bonds!for long bonds 


— | | Puyry | Puyy 4 
0.332 0.138 0.654 
z | 0.152 | | 0.337 | 0.456 
| Puy | | | Puy] | 
x 0.013 0.335 | 0.885 
z | 0.617 | | | 0.924 | 0.897 , 
Ge 
07 07 pe 5 2.69 7 12 
o 4 2—2.69 | 10 14 
| 4 ,—2.67 | 10 14 
4 i—2 70 | 10 14 
oh 4 {2.73 | 40 14 
-3 5 4 10 
i 4 
5, o, 4 10 14 
P4 VIIT| 3 13 16 
O» -5 Of 
1 O 
1 
°F “7 
if 
te —8.8 8.8 
45 


The interatomic distances in the a-plutonium structure for atoms at different positions are given in Table 6. 


The top part of Figure 5 shows the short bonds for atoms in all positions, the continuous lines representing bonds 


lying in one layer and the dotted lines representing bonds between the atoms of neighboring layers (each dotted line 
corresponds to two bonds). 


Atoms at the positions I and VIII (for example, atoms M and N) are joined by a short bond. Considering these 
atoms as a single whole, we find, as at the other positions, a mean number of bonds per atom of 4, and a coordina- 


tion number of 14. Therefore, in the a-plutonium structure, as in all the structures considered above, four covalent 
bonds clearly appear. 


TABLE 8. Dependence of the Coefficient of Thermal 
i i nevertheless vary over certain limits (from 2.52 to 2.73 A) 
Expansion on the Main Crystallographic Directions of ‘ . 
c.-Plutonium and the Temperature in order to characterize their distribution in three mutually 
; perpendicular directions in the lattice, calculations were 
Temperature range, °C made for the projections of the difference d-d_ in these 


_ directions. The results of the calculations are given in 
21-104 | —196--20 Table 7. 


Coefficient of 
thermal expansion 


Table 8 gives values of the coefficient of thermal 
a, | 4 7 expansion along the directions shown in Table 7 in the 
99 14 temperature range 21—104°C [6], -196—20, 20to 100°C[7]. 


It canbe seen from the data of this table that the pro- 
jection of d-d_, both for short and for long bondshas a consid - 
erable value along all three directions of the lattice. There- 
fore, for e-plutonium there isnot such a clearly expressed 
anisotropy of the thermal expansion as for y-plutonium. 


As in the transition from §- to a-neptunium, the difference in the distances 4-4... increases from 0.23 to 
0.48 A, on transition from y-to a-plutonium this difference increases from 0,12 to 0.56 A. In agreement with this, 
the covalent character of the bond on transition from high temperature modifications of plutonium to the lower tem- 
perature modifications increases considerably and the crystalline structure deviates to a greater extent from the struc- 
ture of ideal metallic lattices assembled on the principle of compact packing of the atoms. 


The saturation of covalent bonds in the low temperature region and their relatively uniform distribution along 


the various crystallographic directions means the the anisotropy of the thermal expansion for a- plutonium is expressed 
comparatively weakly. 


With increase in temperature, the weakening of the covalent bonds is accelerated and their distribution in the 
various crystallographic directions becomes more uneven, which finds expression in the strengthening of the anisotropy 
of the thermal expansion. In fact, the anisotropy of the thermal expansion for 6-plutonium, as follows from [6], is 
expressed more strongly than for a-plutonium,and the anisotropy for y-plutonium is expressed more strongly than for 


6-plutonium. The problem of the therma! expansion of 6- and n-plutonium is more complex and requires a special 
study. 


The author is very grateful to S.T. Konobeevskii for discussing the results of this work and for his valuable com- 
ments. 
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THE STRUCTURE AND THERMAL EXPANSION 
OF 6- AND n-PLUTONIUM 


S. T. Konobeevskii and N. T. Chebotarev 


Translated from Atomnaya Energiya, Vol. 10, No. 1, 
pp. 50-57, January, 1961 
Original article submitted April 11, 1960 


Based on the increase in the covalent component of a bond with reduction in temperature in crystal- 
line modifications of uranium, neptunium and plutonium, a new idea is put forward for the structure 
of 6- and n-plutonium. The crystalline structure of these modifications is considered as the result 
of a comparatively small deviation from the ideal face-centered lattice under the action of four 
covalent bonds. This idea makes it possible to explain a number of anomalous properties of 6- and 


n-plutonium, in particular the negative values of the volume coefficient of thermal expansion and 
the high values of the atomic volume. 


Introduction 


The modifications of 5- and n-plutonium have certain specific properties not found in ordinary metals. Of 
these properties, the first which should be mentioned is the negative value of the volume coefficient of thermal ex- 
pansion and the anomalously high atomic volume of these phases. 


According to the existing data the 6-phase has a face-centered cubic lattice with parameter a = 4.637 A (at 
320°C), and an 1 -phase — a face-centered tetragonal lattice with parameters a = 4,701 A and c = 4.489 A (at 490°C). 


The density of these phases is 15,9 and 16,0 g/cm*, respectively, which is much less than the density of the high- 
temperature body-centered ¢ -phase (16.5 g/cm’). 


In accordance with the data of [1] the mean value of the coefficient of thermal expansion of the 6 -phase in 


the temperature range 320-440°C is equal to -8.6-10°% and the coefficient of thermal expansion of the n-phase in 
the various axes has the following values; 


a, = (305 + 35)-10°8 according to the 
a, = (—659 + 67)-10°* | data of [1); 


a, = (444.8 + 12) .10°° | according to the 
a, = (— 1063.5 + 18.2).107% } data of (2). 


In a number of papers, attempts have been made to explain the unusual properties of these phases. In particular, 
in [3] it was suggested that due to the small difference in energy levels of the 5f- and 6d- electrons in plutonium, 


two forms of the atoms could exist simultaneously with a different electron structure and with different dimensions 
respectively, 


With increase in temperature, the number of atoms with a smaller diameter increases due to a reduction inthe 
number of atoms of larger diameter, which leads to a negative value of the coefficient of thermal expansion. The 


tetragonal nature of the n-plutonium lattice, according to this concept, is due to the formation of an ordered struc- 
ture of these two forms of atoms. 


An attempt was made in [4] to explain the negative coefficient of thermal expansion of 5-plutonium by the 
possibility of sharp increase in density of the states in the region at the top of the filled part of the Brillouin zone. 
However, most authors consider that a satisfactory explanation of the mechanism of the anomalous thermal expansion 
of §- and n-plutonium is yet to be found. 


In a paper by one of the present authors [5] it was shown that the anisotropy of the thermal expansion of a num- 
ber of modifications of uranium, neptunium and plutonium is due to a weakening of the four covalent bonds with in- 
crease in temperature and their uneven distribution in the various crystallographic directions. The weakening of the 
covalent bonds means that the minimum interatomic distance din on transition to higher temperature modifications 
increases abruptly and the difference dp) - dyin decreases accordingly until it is equal to zero in the higher tempera- 


ture modification (in the body-centered cubic lattice). The values of d,,j;,, for a number of neptunium and plutonium 
phases are given in the table. 


or the phase - - 
The Value of d,,;,, for Some Phases of Neptunium the phase of 
. dmin is 3.28 and 3.25 A respectively. A comparison with 
and Plutonium 
the data of the table shows that these values of are 


. much higher than those for the highest temperature modi- 
fication — e-plutonium. This lack of correspondence in 
the value of dryjpn for 6- and n-plutonium contradicts the 


Ele- mas dmin: A Ele- 4d min dinin 


above considered regularity. Hence it follows that the 
Np  p Po ly 3.05 structure of these phases is actually more complex. If we 
3.15 


assume that the four covalent bonds appearing in plutonium 
in the a- and y-modifications also occur in the 5- and: n- 
*An approximate value from [6]. modifications, where dj, corresponding to these bonds lie 
— within the limits between these values of dmin for y- and 
€ ~phases (between 3.05 and 3.15A) then it is possible to 
explain not only qualitatively but in the first approxima - 
tion quantitatively the anomalous character of the thermal 


expansion of 6- and n-plutonium. The anomalously high value of the atomic volume and a number of other proper- 
ties of these phases can also be explained. 


y 3.03 | 


Since the structures: of the 6- and n-plutonium are considered as the result of distortion of an ideal tightly 
packed structure under the action of the covalent bond forces, the analysis should start with the higher temperature 
modifications — the n-phase, in which these bonds should appear more weakly. 


Structure and Thermal Expansion of n-Plutonium 


It is well known that a body-centered cubic lattice with a parameter a), can be considered as a tetragonal 
face-centered lattice with parameters ag, and Cs, and a ratio of the axes Cr Jag, = 0.707. The deformation of this 
face-centered cell, accompanied by approach of the ratio of axes Ch /agy to unity, will be a transition from a body- 
centered to a face-centered cubic lattice, On transition to the latter the initial body-centered cubic lattice changes 
to a body-centered tetragonal lattice with ratio of axis cy,/ap, = 1.41. 


The volumes of the initial and final cells will be respectively 
3 


If we assume in the first approximation that the volume of the unit cell does not change during transformation, i.e. 
we neglect the change in degree of compactness of the packing, then ah, = 42a).. Hence it follows that 


a= 73 0.89 a6; “Kr= V2 1.264)... 


Therefore, for complete conversion of the body-centered cubic cell to a face-centered cell it is essential to have 
expansion of the initial cell along the z-axis by gue? and compression along the x-and y-axes by 11%, In accordance 


with this ratio the increases in axes <& = - Dirt = ~2-4. If the transformation stopped at some intermediate stage, 


then the values of Ac and Aa respectively decrease. 


This transformation of the lattice can be accomplished by displacement of the atoms under the action of four 
covalent bonds. A crystallographic analysis shows that the tetragonal deformation of a body-centered cubic lattice 
can be accomplished in two ways, schematically shown in Fig. 1. In the first case (see Fig. 1, a) of four atoms at 
the centers of the unit cells, two atoms are displaced upwards (the + sign) and two others downwards (the - sign). As 
a result of these displacements, the atoms under consideration approach the four atoms of the upper or lower bases of 


the cell and open the hole which is formed by them. A tetragonal structure of the 6 -neptunium type then forms 
with a ratio of axes 5 < 1. 
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In the second case (see Fig, 1, b) the atoms at the centers of the unit cells are displaced in the horizontal plane 
in diagonal directions (the directions of displacements are represented by continuous arrows), In order to preserve 
the symmetry, the cell should be doubled not only along its base but also along its height. The centering atoms of 
the second stage of this cell should then be displaced in directions opposite to the displacements of the atoms of the 
first stage (the directions of the displacement are shown by dotted arrows), The considered displacements of the atoms 
cause an expansion of the unit cell in the direction of the z-axis with corresponding compression in the direction of 
the x- and y-axes. The result is a tetragonal cell with a ratio of axes c/a > 1, 


This process of transformation of a body-centered cubic cell of the €-phase to a tetragonal cell of the n-phase 
n only occur according to the second method of deformation. 


If we take a doubled cell of the n -phase in a body-centered aspect, then its dimensions will be (at 450°C) 


ay = 2-3.326 = 6.652 A 
Cy = 2-4.463 = 8.926 A 


= 1.34. 
It was mentioned above that in deriving the relationship Ac ~ -2.4 no account was taken of the change in 
atomic volume during the transformation of the lattice, A decrease in the degree of compactness of the packing due 
to the considered transition € — 7 should lead to an increase in the atomic volume, It is known that the transition 
€ + n is accompanied by a 3% increase in atomic volume, which corresponds to a 1% increase in linear dimensions, 


In this case, the effective changes in the parameters c and a should be 27 and -10% respectively (instead of 26 and 
-11%), i.e., 


In fact, if the value of the €-plutonium lattice parameter is reduced to its value at 450°C, then, using a value 
of the coefficient of thermal expansion equal to 36,5: 10 © we obtain a, =2a" = 7.261 A, and 


Ac _ Sn —% _ 8.926—7.261_ 
Aa a,—a,  6,652—7.261 


The projection of the n-plutonium unit cell on the (001) plane, in accordance with the ideas put forward, is 


given in Fig, 2, The structure belongs to the space group l4mmn-Djh. The number of atoms in one unit cell is 16, 
The coordinates of the atoms are as follows: 


— 4(d):054; £04 

— 4(e): 002; 00z +( 000; )- 

8Puyy —8(h): 220; 220 
220 3 


min 
the value d,,jp should be closer to the upper limit since there is also a 6-phase between the n- and y-phases, It can 


be assumed that this distance is equal to about 3,12 A. Then x = 0,265; z = 0.255, The deviation of these values 
from 0,250 characterizes the displacement of atoms from the positions of an ideal face-centered lattice, It should 
be mentioned that the change in parameters very sensitively affects the value of d,,;,. For example, with decrease 
in the parameter x by 0,002, dj, = 3.15 A which corresponds to d.,;, of the €-phase, 


It was mentioned above that d__._ for the n-and 6 ~-phases should be between 3,05 and 3,15 A. For the 1-phase 


> —> 


| 
ai 


In accordance with the mean value of the minimum 
interatomic distance drpjin = 3.12 A,we obtain the following 
interatomic distances in the structure of the N-phase (in A): 


) Puy — 8Pu = 3.25; Puyy — 2Pu = 3.14; 


| 

— 4Pu = 3.33; — 2Pu = 3.13; 
| — 4Pu = 3.25; 
— 2Pu = 3.37; 
— 2Pu = 3.52. 

4Pu = => 

Fig. 1. Two ways of tetragonal deformation of a 4 
body-centered cubic lattice under the action of — 4Pu = 3.37; 


four covalent bonds, Here dy, = 3.27 A. 


y2_| | | 3 
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The above considered displacement of atoms from ideal positions of the face-centered lattice should lead to 
the appearance of additional superlattice lines. The structural factor will be determined by the expression 


4f [ cos (2h +l) 4 cos 2ulz + 2cos hx cos 20 ke |. 


The intensity of the superlattice lines should increase with the order of reflection, For example, we will cal- 
culate the relationships between the intensities of two back lines 820 and 822 of the x-ray pattern placed side by side, 
taken in copper radiation, The first of these lines is superlattice, the second is structural: 


4f (14-1 — 2-0.73.0,98) = 4f-0.57; = 4f (— 1 — 1 — 2-0.73.0.98) = — 4f.3.43: 
Te20 (PF _ 8 (4f-0.57)7 
Therefore, the maximum intensity of the superlattice lines is ~1% of the intensity of the main lines, The ob- 
servation of such weak reflections on the high-temperature powder diagrams of plutonium is very difficult. This also 
applies to the structure of the 6-phase considered below. 


We will now explain the behavior of the n-phase with change 
in temperature, 
<a y When the 1)-phase is heated in the region of its existence, the 
rs above considered structure transformation (in the reverse order) will 
occur only partially, however the ratio between the changes in the 
parameters c and a, appearing formally as a process of thermal ex- 
pansion should remain about the same, i,e.,- 2,7. Since 


Ac =ca,AT; Aa = aa, AT 


(a, and a, are the coefficients of thermal expansion in the direction 
of the a-and c-axes, AT is the difference in temperature), 


Ac @c 


Aa 


In accordance with the data given above [1,2] for a, and a, 
the ratio Ac/Aa is obtained equal to -2.9 and -3,5, which agrees 


airly well with the theoretical ratio of -2.7. 
Fig. 2, Projection of the n-plutonium unit 


cell on the (001) plane, Height in fractions * If, in the determination of coefficients of thermal expan- 


of a parameter c: sion of the -phase we use the latest data [2],then in the range of 
O-0 and ,o- 0.005 and 0,495; existence of this phase(451-480°C) 
and -0,005; o- %. 
Ac = — 8.926. 1063.5. (480 — 451) =— 0,275 A; 


Aa = 6,652 - 444.8. 10°* (480 — 451) = 0.086 A. 


Since for the case of total transformation (7 + ¢) 


= Ae (480) — Cn (451) = 7.269 — 8.926 = — 1.657 A; 
Adeotal = Ve (480) — An (451) = 7.269 — 6,652 = 0.617 A, 


the values given above comprise from Ac, ,,, and Aa, . 4), 16.6 and 14,0% respectively, Hence, it follows that 
heating the 7-phase in the range 451-480°C is accompanied by its structure approaching that of the €-phase by.about 
15%, The remaining part of the lattice transformation occurs in the process of the phase transition n >€. 


Structure and Thermal Expansion of 6-Plutonium 


We will consider a unit cell of the n-phase in a face-centered aspect. For this purpose on the projection of 
the unit cell of the n-phase doubled along the x-and y-axes on the (001) plane, we will divide the diagonal cell with 
a base ABCD (Fig. 3). The dimensions of the cell are a= 9,408 A and c = 8,926 A, The cell consists of eight sym- 
metrically equivalent cells of the type BEFG with dimensions a* = 4,704 A and c' = 4,463 A, Each of the atoms of 
the type I-IV has four short bonds at a distance 3,11-3,13 A, With reduction in temperature these distances continue 


— 
he 
1 
D, D 
~3.13 
| v7) 
y 
: 
| D, ° 3 D; 
3,52 
| 
4%, 
as 83 
50 


to decrease under the action of the covalent bond forces, It can be assumed that during the transition to the region 
of the 5-phase djjjn = 3.09 — 3.10 A, 


A decrease in the distance of the type I-II from 3,13 to 3.09 A, i.e.,by about 1.3%, should lead to a decrease 
in the parameter a‘ by the same value, i.e.,from 4,70 to 4.64 A, If in the first approximation we neglect the change 
in volume during the transformation, then the parameter c’ should increase by 3,1% (based on the relationship dc, 
~2,4), i.e., from 4.46 to 4.60 A. If we bear in mind that further displacements of atoms of the type I-IV shout also 
lead to further reduction in the degree of compactness of the lattice, then the value of c' should be obtained some- 


what higher, As a result the considered transform ations lead to the condition c’ = a’, i.e, to the development of a 
cubic lattice, 


The unavoidability of satisfying the condition c’ = a’ 
also proceeds from the following. The atoms V and VI, like 
the atoms I-IV, center the upper and lower bases of the unit 
cells BEFG, However, while the unit cells remain tetra- 
gonal, atoms of the type I-IV and the type V-VI are not 
structurally equivalent, The displacement of these atoms 
considered above leads to the elimination of their structural 
nonequivalence, as a result of which the distances of the: 
type II-III should become equal to the distances of the type 
III-V and the distances of the type II-B equal to distances 
of the type V-B, The satisfying of these conditions means 
the presence of a cubic cell. 


As a result of this transformation, we arrived at a 
structure of the cubic cell of the 5-phase of plutonium with 
parameter a‘ = 4,632 A (at 450°C), As mentioned above, a 
certain decrease in the compactness of the packing should 
lead to a small decrease in density during the n > 6 transi- 
tion, It is known from experimental data [7] that this re- 


Fig. 3. Projection of four unit cells of n-plutonium duction is 0.4%, 
The projection of the unit cell of 5-plutonium onthe 
meter c: @- 0 and 
+ (001) plane in agreement with the considered ideas is given 
®— 0,005 and 0,495; @— 0,505 and 0,005; 
@-%; o- ¥, in Fig. 4. The structure corresponds to the space group F23 
45 OC” 


—-T’, The lattice parameter a = 2a" = 9,274 A; the number 
of atoms in the unit cell is 32, The positions of the atoms 
are as follows: 


4P uy — 4 (a) : 000 
444 


24P — 24 (g) 


/ oO The value shown above = 3,09 A corresponds to 
O 


the value of the parameter x = 0,014, which means the dis- 
placement of atoms from the centers of the faces of an 
x ideal face-centered unit cell by 0,13 A, i.e. by about 4% 


+ ( 000; 440). 


from din’ With this value of x the interatomic distances 
Fig, 4, Projection of unit cell of 6-plutonium on siniecntaaa.: 


the (001) plane. Height in fractions of a parameter 
c: @-0 and O-'4; @— %; and 0,486 Pu; — 12Pu = 3.28; Puy; — 12Pu = 3.28; 
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Puyyy — 4Pu = 3.09; 
— 4Pu=3.28: — 4Pu = 3.46. 


The value dry = 3,28 A, i.e., is equal to that value which was assumed constant for all atoms, 
We will now consider the behavior of the 5-phase with change in temperature. 


The values d_;,, and dja, can be expressed by coordinates of atoms in the following way: 


= 2(4-2)a; (1) 


In order to establish the character of the temperature dependence of the lattice parameter we will consider 
the relationship 


(2) 


The derivative from this expression with respect to temperature will be the coefficient of thermal expansion 
da 06d min 


9a _ dmin?= 


With reduction in temperature accompanied by an increase in the covalent bond forces there will be a reduction in 
dmin 2nd an increase in the parameter x, The first term of the right side of equality (3) will be a positive compo- 
nent and the second term will be a negative component of the coefficient of thermal expansion, We will designate 


} 2( )or 2( 7-2) a7 


in accordance with which 


where B is a positive component and y is a negative component of the thermal expansion coefficient. The ratio of 
these components is 


admin 


min 
In connection with the fact that the value 4 -x differs little from Y we can assume 
Ox 


od min’ 
d-min 


If the parameter x is expressed in absolute units, then the expression assumes the form 
dx 
a dmi 
@dmin 
1 


Since a = 9,274 A, and dmin = 3.09 A, then dmin i] 3 
a 


ence — = ————._, for 7 ==6 . Substituting the value of Y in expression (4) we obtain 
B 3 Odmin 3° 8dmin 


4 Ox 


min 


= 
52 


and for finite increases 


4 Az 
B [ 3 Admin (5) 
From a comparison of the n- and 5-phases it follows that the value of the displacement Ax in the n- 6 transition 
for eight atoms is approximately equal to 0.03 A, and for four atoms to 0,09 A, ' In the temperature range of existence 
of this phase, as follows from an analysis of a number of uranium, neptunium and plutonium phases, the value of the 
displacement should be less, It can be assumed that for the regions of the 5-phase Ax = 0,02-0,03 A. The change in 


dmin is also equal to a value of this order, Hence it can be assumed that in the temperature region of the 5-phase 


= 1, In this case 
Admin 


The positive components of the thermal expansion coefficient can be taken equal to the mean value between the co- 


efficient of thermal expansion of the €-phase (36,5- 107) and the coefficient of thermal expansion of the y-phase 
(34.5°107§), i.e.,35.5 x Then 


The obtained value agrees well with the experimental value of the thermal expansion coefficient of the 6-phase, 
equal to -8,6° 107°, 


We will now consider in what way the thermal expansion coefficient of the 6-phase should change at a tem- 
perature less than 320°C, i.e.,in the case of supercooling of the phase, With further reduction in the temperature, the 
value of X cannot increase to any great extent since then, in accordance with the second equation of expression (1), 
the value of d,,,, would increase strongly. In [5] it was shown that the value of dja for a number of structures 
changes very little with temperature, On the other hand, the value of d_,;,, during transition in the region of low 
temperature modifications, as can be seen from the table, changes to a considerable extent, As a result, with super- 
cooling the value of Ax/Ad__._ should tend to zero and the value of & should tend to the value of 8. Hence, it 
follows that with supercooling of the 5-phase, its thermal expansion coefficient should become at first equal to zero 
and then to a positive value, This conclusion agrees well with the data of [8] on the determination of the thermal 
expansion coefficient of the 5-phase of alloys of plutonium with aluminum, 


The addition of an alloying element, weakening the covalent component of the plutonium bond, should lead 
to a reduction in the values of Ax, As a result, according to equation (5) the thermal expansion coefficient of the 
6-phase during addition of the alloying element should change from negative values to positive values in the high 
temperature region also, The reduction in the difference between d,,;,, and d,,,, caused by the decrease in the value 


of x will also weaken the covalent character of the bond in the structure and facilitate fixations of the alloyed 5-phase 
at room temperature, 


SUMMARY 


1, In accordance with the concept of increase in the covalent component of a bond with reduction in tem- 
perature, we have considered the most probable character of change instructure during the aleotropic transformations 
€ +7 and n-6 of plutonium. It has been suggested that the ideal face-centered lattices for N-and 6-plutonium are 
only the first approximation of their structure, Actually the structure of these phases is more complex and is charac- 
terized by a displacement of atoms from positions of the ideal face-centered lattice by a value of about 3-4% from 
dyy+ This displacement leads to a reduction in the compactness of the packing in these structures and to an increase 
in their atomic volume, 


2. It has been shown that the strong increase in the parameter a and the still greater reduction in the para- 
meter c of the 1-phase which occur with increase in temperature, formally taken as the result of thermal expansion, 
in actual fact should be considered as a process of continuous transformation of the structure, bringing it nearer to the 
structure of the €-phase in accordance with weakening of the covalent component of the bond, The difference inthe 
lattice dimensions of n- and €-plutonium in the range of the n-phase region (451-480°C) is reduced by about 15%, 


3, It has been shown that the change in the lattice parameter of the 6-phase with reduction in temperature 


can be considered as the result of simultaneous superposition of two opposite processes, The first of them is the re- 
duction in d,,,;,, under the action of increase in the covalent component ot the bond and a reduction in the vibrational 
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amplitude of the atoms and leads to a reduction in the lattice parameter (the positive component of the thermal ex- 
pansion coefficient), The second process is the increase in displacement of atoms from the ideal position and leads 
to a reduction in the degree of compactness of the structure and to an increase in the lattice parameter (the negative 
component of the thermal expansion coefficient), It was shown that in the region of high temperatures the predomi- 
nating role isplayed by the second process and the overall value for the coefficient of thermal expansion is negative. 
During supercooling of the 5-phase the first process begins to play a more important role and the thermal expansion 
coefficient becomes at first equal to zero and then to a positive value, 


As a result of alloying of the 5-phase, accompanied by a weakening of the covalent component of the bond, 
the role of the second process should decrease and the coefficient of thermal expansion of the 5-phase during 
addition of the alloying element becomes a positive value not only in the region of low temperatures, but also in the 
region of high temperatures, 


LITERATURE CITED 


F, Ellinger, J. Metals, 8, No. 10, 1256 (1956). 

A. Coffinberry et al, Proceedings of the Second International Conference on the Peaceful Use of Atomic 
Energy (Geneva, 1958), Selected Reports of Foreign Scientists, Vol. 6, Nuclear Fuel and Reactor Materials 
[in Russian}, Moscow, Atomic Energy Press, 1959, p, 157. 

M. Waldron, Atomics, October, 383, (1957), 

J. Varley, Proc, Roy. Soc,, A237, 413 (1956), 

N. T, Chebotarev, Atomnaya Energiya, 10, 1, 40 (1961). 

W. Zachariasen, Acta crystallogr., 12, No. 3, 175 (1959). 

E, Jette, J, Chem, Phys.,,23, No, 2, 365 (1955). 

M. Waldron et al, Proceedings of the Second International Conference on the Peaceful Use of Atomic Energy 
(Geneva, 1958), Selected Reports of non-Soviet Scientists, Vol. 6. Nuclear Fuel and Reactor Materials [in 
Russian] (Moscow, Atomic Energy Press, 1959) p, 142. 


1, 
4 3. 
4, 
5. 
6. 
1, 
8. 
4 54 


SOME PROBLEMS IN THE LOCALIZATION OF RADIOACTIVE 
ISOTOPES IN CONNECTION WITH THEIR SAFE BURIAL 


P. V. Zimakov and V. V. Kulichenko 


/ 
Translated from Atomnaya Energiya, Vol, 10, No. 1, 
pp. 58-63, January, 1961 
Original article submitted April 8, 1960 


It is shown that radioactive fission fragments can be safely buried by converting the liquid wastes 
into solid materials, From the point of view of providing maximum localization of the fragments 
in a minimum volume of materials such as glass there are undoubted advantages, An analysis is 
given of the physicochemical principles for preparing these materials and certain features are con- 
sidered in the behavior and state of radioactive fission fragments included in them, By using 
methods of fractional chemical and x-ray structural analyses we have shown the heterogeneity of 
the structure of the fused materials and the related features of behavior of fission fragments during 
leaching. It has also been shown that the ionizing radiation helps to increase the crystalline frac- 
tion in the materials, 


A deciding factor in the safe burial of radioactive wastes is their reliable localization at a suitable place for 
a long period, regardless of the place chosen or the means of burial, The main danger in any “cemetery” is the 
threat of possible delocalization — the buried radioactive isotopes being scattered from it into the surrounding medium. 


The most important problem at the present time is to bury radioactive fission fragments which form in large 
quantities in many countries in various nuclear power units, 


Considering in the wastes of nuclear installations the presence of isotopes such as Sr™, Cs'87, pu®®, it has 
been calculated [1-3] that the time of strict localization needed is several hundreds of years,while the danger of 
radioactive contamination of the surrounding medium hardly ever disappears, 


Most workers have now recognized that the burial of liquid radioactive wastes in containers of any material, 
regardless of the place and means of burial, is dangerous, Nor is it safe to pour liquid radioactive wastes into the seas, 
rivers or lakes [4]. 


A container of any material containing radioactive solution cannot withstand corrosion and hence cannot pre- 
vent radioactivity passing through the surrounding medium, Gaseous products of radiolysis, unavoidable in liquid 
mediums, can increase the pressure and accumulation of explosive gaseous products, In some radioactive solutions 
the gas evolution per day can reach 10 cm’ per 1 curie, Furthermore, liquid radioactive wastes, especially highly 
active wastes, can dry up during self-heating; this leads to steam or pyrochemical ejections and explosions [5]. Putting 
liquid wastes into deep wells is no guarantee that the radioactive isotopes will not come into contact with the out- 
side layers of the soil and the soil waters, It has been shown that the safest burial of radioactive fission fragments 
is to convert the liquid wastes into solid materials which can firmly retain the radioactive fragments included in 
them [1,3,6-12]. 


A comparison of the characteristics of various solid preparations containing fission fragments (Table 1) shows 
that the glass type preparations are the most suitable for safe burial since this material can reliably localize the 
radioactive isotopes within itself. This material is radiation- andthermally stable and gives the greatest reduction 
in volume of the radioactive material which has to be buried. The direct preparation of this material is preceded 
by the concentration of radioactive fragments of liquid wastes by sorbing them on difficultly soluble precipitates 
(hydroxides, salts of phosphoric, oxalic or other acids) or by concentrating to obtain a mother liquor usually con- 
taining a large quantity of salt (for example, nitrates), The obtained concentrates (in the form of salts or hydrated 
oxides) in which the active fragments are distributed contain an excess of water ~ free, crystallization and consti- 
tutional water, The complete removal of water is the first and main condition for giving a radiation-stable prepa- 
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Table 1, 


— 


ration, The second condition is the thermal decomposition of the radiolyzed salts (for example, nitrates) and the al- 
most complete conversion of inert and radioactive elements to the oxides. 


The Characteristics of Some Solid Preparations Containing Radioactive Fission Fragments 


Characteristics 


Oxides of 
metals dried 


at 200-250°C 


[1,8,9] 


Oxides of 
metals 
roasted at 
500-700°C 
[1] 


Roasted clays 
(1, 7, 9, 13) 


Cemented 
blocks 
(3, 12] 


Molten Materials such 
salts as glass 
13] [1, 9, 14) 


Difficulties in 
preparation 


Change in volume 
compared with 
initial volume 


Liberation of gaseous 
products due to 
radiolysis on storage 


Thermal stability, °C 


Thermal conductivity, 
keal/mehr* °C 


Washing out of frag- 
ments in soil waters 
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strong dust 
formation 


100-fold 
reduction 


0.05 


considerable 


dust 
formation 


100-fold 
reduction 


0.03 


noticeable 


limited exchange 
capacity of clays 


50-fold 
reduction * 


up to 900 


0.05 


small 


necessity 


30-50% 
increase 


up to 100 


no data 


small 


* Depends on the exchange capacity of the clay, degree of salt saturation and acidity of the solutions, 


Figure 1 shows a thermogram for the dehydration of a hydroxide precipitate (pulp) obtained in the treatment 
of solutions containing fission fragments, Free and adsorbed moisture is removed at low temperatures (up to 150°C) 
and is accompanied by a noticeable endothermic effect, At higher temperatures (up to 300°C) this effect is caused 
by the removal of capillary adsorbed water and the dehydration of some less thermally stable hydroxides (for example 
chromium hydroxide), The exothermic effects observed in the temperature range 300-600°C should be explained by 
recrystallization of some hydroxides, The dehydration of hydroxides continues to temperatures of 700-800°C; on the 
thermogram this is represented by several endothermic inflexions, The iron and manganese oxides in the preparation 
dissociate thermally at higher temperatures, which is shown up on the thermogram by the endothermic effects [14]. 


for mixing 


corrosion of 
apparatus 


high tempera- 
tures ( >1000°C) 


50-70-fold 
reduction 


350-fold 
reduction 


~1 to 2 


soluble very small 


Whereas the thermal treatment of the precipitated hydroxides of the metals at temperatures up to 800°C 
gives radiation-stable preparations, the same treatment of, for example, nitrate mother liquors does not give com- 
plete decomposition of the nitrates, In this case, displacement reagents should be used. Thus, if the appropriate 


amount of silicic acid is added to the mother liquor, there is almost complete decomposition of the nitrates with 
the liberation of oxides of nitrogen. 


This treatment still does not provide sufficiently stable localization of the fragments, especially those such 
as Sr® and Cs**’, the free oxides of which are very soluble in water. 


As can be seen from Table 1, the chemical stability of the preparation can be increased by including the 
fragments in a solidified mass such as glass, For this purpose, the obtained oxides are mixed with various inert addi- 
tions which can form comparatively readily fusiblechemically stable compounds (for example, silicates, on the addi- 
tion of SiO) or can help the solution of high melting components in the liquid phase of readily fusible additions (for 
example B,Os). Inert additions are chosen so as to give the preparation maximum localization of the fragments; 
easily fusible compounds such as B,Os, NagCOsg, etc,, are therefore added in small amounts. The main components 
giving minimum solubility of the fused preparations are aluminum oxide, silicon dioxide and some others, 


: 

4 a Yes No No Yes Yes No 

> 1000 > 1000 | upto 300 | >1200 


It should be borne in mind that the high-temperature treatment of radioactive materials converts a part of the 
radioactive fragments into the steam-gas phase and aerosols [10-12], As shown by our investigations, the main rea- 


son for the formation of radioactive aerosols is the sublimation of readily volatile compounds, mainly isotopes such 


as and Nb™, 


Fig, 1, Thermogram of the dehydration of 
hydroxide pulp. 
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Fig, 2, The formation of aerosols during the 
thermal treatment of radioactive pulps, 


The liberation of radioactive aerosols increases sharply (Fig, 
2) at temperatures above 1200°C, Care should therefore be taken 
when preparing high- melting preparations, especially in the treat- 
ment of highly active wastes, Suppressing the formation of aerosols 
is the object of later investigations, Success can be achieved by 
firmly locking the fragments in asolid melt, by selecting the fusion 
conditions and composition of inert additions [11,12]. 


Figure 3 shows the behavior during fusion of one of the mix- 
tures based on a radioactive hydroxide precipitate, As can be seen 
by comparing the thermogram, the curves for the change in weight 
of the mixture and solubility of the obtained material, at the fusion 
temperature the decomposition of thermally unstable compounds 
continues and the obtained preparations do not have the maximum 
chemical stability for the given composition. The working tempera- 
ture therefore usually exceeds the melting point and is determined 
separately for each actual mixture [14], 


The preparations obtained by high temperature fusion are 
silicate structures with inclusions of various compounds (such as 
ferrochromates and other simple and complex oxides), 


The washing out of radioactive fragments from these struc- 
tures can vary between the limits 10*— 10°*% (depending on the 
composition of the flux additions) from 1 cm” to a depth of 2 mm 
per day. The reduction in washing out of radioactive fragments 
compared with the washing out of roasted oxides is due to the fact 
that as a result of the high temperature treatment in the presence 
of silicic acid in the mixture various silicate compounds of sepa- 
rate fission fragments form, the chemical stability of which is much 
higher than the oxides, 


The oxides of other radioactive isotopes, not subjected to 
serious chemical changes at the selected temperatures due to the 
presence of a silicate medium surrounding them, become less 
accessible for the solvents, For example, at 1100°C and a “boiling 
time" of 2 hours, the formation of strontium silicate could hardly 
be expected, For this formation, higher temperatures and much 
longer boiling time for the melts would be needed, If the washing 
out of Sr® or Sr® from the preparation obtained at 1100°C is 


6+ 107%, then for this preparation "boiled" at 1350°C for two hours, the washing out is reduced to 3- 10™“% and when 
boiled for four hours, it is reduced to 5- 107%, Boiling at 1350°C for six hours reduces the washing out of Sr® and 
se” by one order, In addition, the washing out of Sr® and Sr® from molten preparations obtained even at 1100°C is 
several times less than the washing out of these isotopes from a simple roasted powder of the oxide, However, in- 
creasing the temperature above 1200°C causes a sharp increase in the change of radioactive isotopes to aerosols and 
to the gaseous phase; therefore operation at 1350°C and excessive increase in the time of boiling of the melts are 
undesirable, It would be better to strengthen the structure of the silicate medium containing the fragments, 


Figure 4 shows the time dependence of the washing out of B-activity from one of the fused fragment concen- 
trates by acid and alkali salt waters, The sharp reduction in the transfer of fission fragments to acid water with time 
is explained by the formation on the surface of the preparation of a film of silicic acid, protecting the specimen 
from further destruction, Alkali waters distrub the silicon framework, causing an increase with time in the washing 


out of the fragments from vitreous concentrates, 
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Fig. 3, The behavior during fusion of one of 
the mixtures based on a radioactive hydroxide 
precipitate: 1) thermogram; 2) curve showing 
change in weight during thermal treatment; 

3) solubility in water of a preparation obtained 
at corresponding temperatures [14], 
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Fig. 4. Time dependence of the washing out of 
B-activity from a fused fragment concentrate by 
acid and alkali soil waters: 1,2) washing out at 
pH of 8 and 4 respectively, from a preparation 
obtained without additional annealing; 3) washing 
out after 2 hr annealing in water at pH = 4. 
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It should be mentioned that for vitreous fused preparations containing active fragments, the value of the chem- 
ical stability of the preparation cannot be identified with the value of washing out of radioactive fragments from it. 
As shown by the results of successive treatment of a powdered preparation by various reagents at the temperature of a 
boiling water bath (Table 2), compounds containing fission fragments behave somewhat differently from the inert com- 
ponents which make up the main mass of the preparation, The heterogeneity of the structure of fused fragment pre- 
parations is confirmed by the results of an x-ray diffraction study (Fig. 5). 


With rapid cooling of the melt, a brittle vitreous pre- 
paration is obtained which contains no crystalline inclusions 
(see Fig. 5: 1), Its strength can be increased considerably by 
additional “annealing” at 400-600°C, Together with the 
reduction in the structural stress, there is crystallization of 
some compounds (mainly iron oxides) contained in the fused 
preparations (see Fig, 5: 2), 


The following conclusions can therefore be drawn, 
The fused preparations obtained in the treatment of highly 
active wastes cannot exist in the purely amorphous state; 
they will always contain crystalline inclusions, mainly com- 
pounds of iron, which are basic fragment carriers, In the 
specimens investigated, most of the silicon compounds re- 
main in the amorphous state, 


An inevitable propeity of radioactive isotopes is their 
generation of heat due to radioactive decay, Figure 6 shows 
the extent of the contributions of the separate fission frag- 
ments to the total heat liberation of the mixture [6]. The 
temperatures developed in the preparation, including such 
isotopes, also depend on the thermal conductivity A and 
the specific heat of the preparation C, The dependence 
of the self-heating temperature of cylindrical fused frag- 
ment concentrates on their specific activity, diameter of 
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Fig. 5, Results of an x-ray study of fused preparations 
obtained in the rapid cooling of the melt: 1) before 
annealing and radiation (lines are not noticeable, 
amorphous structure); 2) after annealing at 400°C; 

3) after Co™ radiation; 4) after radiation with simulta - 
neous annealing, I) SiO, (quartz); IT) y(Fe,O3); III)FegO,. 
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TABLE 2, The Fractional Solubility of Vitreous Radioactive Preparation in Powder in 
Various Mediums 


Transfer of fragments into 
- 
ime of treat Lenin went solution, % 
Reagent ment at100°C, % 
ne regar 
B -radiation y -tadiation 


Water 13.4 0,03 0,04 
NaOH (20%) 42 11.3 2.5 1.4 
CHsCOOH(1: 1) 12 49,2 93.7 39.8 
HC1 (1:1) 18 5.8 3.3 41.0 
Insoluble residue = 20.3 0.5 17.8 


: TABLE 3, Self-Heating of Fused Fragment Prepara- 
10 tions (X = 1.1 keal/m-hr- °C, C = 0.2 kcal/kg + °C) 


Self-heating temp. , °C 


Specific 
activity, | Diameter.|”” with age of fragments 


curie/liter| 75 180 
10 days days year years 


Tor al 30 000 45 ; 
5910] 6670 | 2580 | 490 
490| 206} 78| 45 


229] 96 38 9 


— 


0,277 Cal/hr. ton U 


0 


Time after unloading reactor, years 


Fig. 6. Change in heat liberation of separate fis- 
sion fragments with time (enrichment 1.2%, power 
of reactor 5,5 Mv/ton; radiation level 2750 Mw/ton). 


the cylinder and age of the fragments is shown in Table 3, The data of Table 3 were obtained by calculation, the 
thermophysical characteristics C and \ being determined experimentally. 


The calculation was checked by an experimental determination of the self-heating temperature of the natural 
fused concentrates of the fragments, 


It can be seen from Table 3 that by selecting the appropriate dimensions of the preparations, the danger of ex- 


cessive self-heating can be avoided even with the presence of high levels of activity. The self-heating decreases 
strongly with time. 
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A METHOD FOR DETERMINING DOSES IN THE INHALATION 
OF RADON DECAY PRODUCTS 


I, I. Gusarov and V. K. Lyapidevskii 


Translated from Atomnaya Energiya, Vol. 10, No. 1, 
pp. 64-67, January, 1961 
Original article submitted January 12, 1960 


This article gives a new method for determining the absorbed energy in the inhalation of radon decay 
products from the number of a-particles emitted during the total decay of daughter products of radon 
separated from 1 liter of air. A method is proposed for determining the concentration of radon from 


the number of a-particles emitted due to total decay of the daughter products collected on a filter from 
a continuously purified volume. 


The determination of doses in the inhalation of radioactive material is a very complex problem. In fact, its 
solution requires information not only on the quantity of each separate radioactive substance introduced into the or- 


ganism but also its distribution within the organism. The latter depends on many widely fluctuating factors and it 
can therefore hardly be studied in a general form. 


In addition, it is obvious that since the amount of material taken by the organism during inhalation depends on 


the quantity of these materials in the air, the study of the whole complex of problems should rest on accurate data as 
to the content of radioactive substances in the air. 


Calculation show that during the inhalation of radon decay products the dose is mainly determined by the a- 
radiation, since the energy and biological activity of the 8-radiation compared with the a-radiation is negligibly 
small, The main attention is therefore paid to the study of a-active products [1]. 


The methods now used for the determination of radon decay product concentrations in the air are based on the 
determination of RaA, RaB and RaC. The latter is achieved either by analyzing curves for the decay of radioactive 
substances collected on a filter from a given volume of air [2] or by analyzing the amplitudes of impulses given by 
the a-particles leaving the surface of the filter [3]. These methods are complex, which prevents their widespread use. 
Furthermore, with low contents of radioactive substances in the air, close to the natural background, the analyses of 
the decay curves and the impulse amplitudes are hindered due to the high statistical scatter of the experimental 
points. The method considered in this article is largely free from these faults. 


We will determine the energy formed during the total decay of a-active atoms separated from 1 liter of air: 


E,=a (ERaA +E RaC’) + (b-+c) ERaC’, 


where Ey is the total energy; a, b, c are the number of atoms of RaA, RaB, RaC (in 1 liter of air) respectively; FRaA 
and ERaC are the energies liberated during one act of decay of the appropriate element. 


Formula (1) can be written in a different way if we average the energy of a-radiation of RaA and RaC’. We will 
represent this value by E,. Then 


E, Eq (2a b+ c) = Ean, (2) 


where n is the number of a-particles emitted during total decay of radon products separated from 1 liter of air. 


Calculation shows that if we assume Eg = 7.35 Mev, then the deviation of E, from E, with a change in the 
ratio aA, : bAg: cAg (Ay, Az, and, Ag are the decay constants of RaA, RaB and RaC, respectively) from 1:1;1 to 
1:0,1:0,01 varies within the limits # 3.5% from E,. In view of the fact that in the important practical cases the 


deviation from the equilibrium concentrations varies by not more than the limits shown above, the replacement of 
E, by E, can be considered valid. 
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We will assume, as is usual, that in the inhalation of decay products there is no selective absorption of radio- 
active materials, i.e., the coefficient of capture 7 is the same for each of the radon decay products, Then the total 
energy of radioactive decay liberated within the organism on the inhalation of 1 liter of air is 


E= mEq. (3) 


These discussions show that the initial value for the termination of doses in the inhalation of daughter products 
of radon can be n, i.e.,the number of a-particles emitted during complete decay of radon daughter products, separa- 
ted from 1 liter of air. 


We will deal with different methods for determining the value n. 


1, A fixed amount of air is drawn through a filter with a known filtration efficiency %,, The a-particles leav- 
ing the surface of the filter are recorded by an apparatus. Depending on the geometry of the count and the self-ab- 
sorption in the filter, the efficiency of recording 2 can be different. Then, establishing a continuous count of par- 
ticles from the start of filtration and up to complete decay of the radioactive materials collected on the filter, it is 
possible to determine n: 


Np 


where N, is the total number of a-particles recorded by the apparatus; v is’ the rate of evacuation (liter/ min); tis 
the time during which filtration takes place (min). 


2. The second method is based on the assumption that if, with a fixed rate of evacuation v, the number of a- 
particles n, recorded in unit time does not change over a large interval of time, then during time At within this in- 
terval, the number of radioactive atoms remains constant. In other words, during the time At the number of radio- 
active atoms settling on the filter is the same as those which decay during this time. 


We will assume that during time At one liter of air passesthrough the filter and that it collects ?, (atb+c) atoms 
of RaA, RaB, RaC, According to the condition,the same number of atoms should decompose during this time. The 
corresponding total number of decays ¥4(2a+b+c) = yn. Hence, it can be seen that if the count rate is equal to Np 
impulses per minute and the rate of evacuation v liters per minute, then the total number of decays is 


To use the above method for determining doses, in collaboration with A. M. Konstantinov, we developed an 
apparatus which could continually record a-particles leaving the filter. One form of the apparatus is a somewhat 
modified scintillation attachment to the B type apparatus. 


A fixed volume of air recorded by a gas counter was drawn through the filter, the a-particles being recorded 
directly from the moment of start of filtration. The rate of evacuation was measured with a rheometer. The number 
of impulses was counted every minute, The total number of impulses was recorded simultaneously. 


In the first series of experiments, the filtration was carried out during time t and the number of decays from 
the start of filtration and up to complete decay of the radioactive material was recorded directly by a mechanical 
counter. In the second series of experiments, air was drawn for a long time through the filter with a constant rate 
so that the count rate Np was constant. 


The first case corresponds to a short stay of the organism in the contaminated atmosphere, the second corres- 
ponds to a long period when a more or less constant concentration of radioactive materials is established in the or- 
ganism, different from their concentration in the inhaled air. 


The amount of suspensions in the atmospheric air varies with time; the measurements were therefore made 
with two identical instruments, operating simultaneously. The air samples were taken from the same volume. 


It was established that the number of a-particles emitted during the total decay of daughter products of radon, 
separated from 1 liter of air is almost independent of the method of determining the radon daughter products. Thus, 
despite the fact that the relationship RaA :RaB :RaC on the filter depends on the rate and time of filtration, the total 
number of decays (on the same filter) is the same in all cases, if the same amount of air is drawn simultaneously 
through both filters. 


n 
a= 
UPiP2 
x 
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As already mentioned, when determining doses from the radon decay products, the initial value is usually the 
total activity of RaA, RaB and RaC in the air (curie/liter), With the same activity, the total number of a-particles 
and hence the dose obtained by the organism can vary 2-3 fold depending on the relationship RaA :RaB:RaC. This 
relationship depends to a large extent on the conditions of the air exchange, the amount of dust in the air and other 
factors; therefore, the total activity of the decay products does not accurately determine the dose. The determination 
of the relationship RaA :RaB:RaC in the chain of decay products from the decay curves or by spectrometry is complex. 


In view of this, in the determination of the dose from 
radon decay products the initial value should not be taken as 
the total activity of these substances in the air, but the total 
number of a-decays n, which is given by the decay products 
separated from 1 liter. 


A similar method can also be used to determine the 
dose from daughter products of radon forming in the organ- 
ism itself. For this purpose the air which has first been puri- 
fied by filter 1 (Fig. 1) from decay products is continuously 
drawn through the container A. The radon daughter products 
which are formed in the container A are trapped by the filter 
2. By means of scintillation counters 3 and 4, there is con- 
tinuous recording of the a-particles leaving the surface of 
the filters 1 and 2. 


Fig. 1. Apparatus for determining concentration of 
radioactive materials in the air (the arrows show the 
direction of movement of the air). 


With regard to the total number of a-particles leaving the filter 1 and 2 it is possible to decide on the ratio of 
the does obtained by the lungs in the inhalation of radon daughter products contained in unpurified air to the dose ob- 
tained by the lungs due to daughter products formed from radon in the lungs themselves, It is assumed that the effi- 
ciency of the filters 1 and 2, the volume of air drawn through the filters and the volume of the vessel A are known. 
The coefficient of capture of the daughter products by the walls of the vessel A is also known. 


Figure 2, a and b, gives curves for the collection and decay of radon daughter products collected on filters 1 
and 2, It can be seen that the total number of decays on the filter 1 is much greater than on filter 2 (at a rate of 
evacuation of 15 liters/ min and with the volume of the vessel A equal to 1 liter). These data correspond to the state- 
ment that'the main danger in inhaling radon and its decay products is connected with the action of radon daughter 
products in the inhaled air. 


The radon concentration (curie/liter) can be determined from the number of o-particles emitted due to total 
decay of the daughter products collected on the filter 2 during time t. In fact, each atom of radon which decays in 
the volume A will give two a-particles on the condition of complete collection of daughter products on the filter. If 
under these conditions the number of impulses in the total decay of daughter products is equal to N, and the efficien 
cy of recording », then the number of radon atoms N = ap. . The activity (curie/liter) can then be expressed by the 
formula 2 


7’ 


10718 
“= "44QatV ’ 


where t is the time of collection on the filter 2 of daughter products formed in the vessels A(min);V is the volume 
of the vessel A (liter). 


The given formula holds when the number of decays of RaA in the volume A and the amount of daughter pro- 
ducts settling on the walls of the vessel during the evacuation are negligibly small. In the opposite case for the same 
radon activity, the total number of decays will depend on the rate at which air is drawn through the vessel A. An 
experimental study showed that at a rate of evacuation greater than 15 liter/min, there was no dependence on the 
rate, i.e. under these conditions the decay of RaA in the volume A and settling of the decay products on the walls of 
the vessel can be neglected. 


The described method for determining the activity of radon by measuring the total number of oa-decays of 
daughter products settled on the filter from a measured volume can also be used to measure short lived radioactive 
gases, for example rhoton or actinon. 


| 
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Fig. 2. Curves for collection and decay of radon daughter products: Np is the count rate; t is the time 
of evacuation, 
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LETTERS TO THE EDITOR 


AVERAGE NUMBERS v AND 1 OF NEUTRONS IN THE FISSION 
OF AND BY 14 Mev NEUTRONS 


N. N, Flerov and V. M. Talyzin 


Translated from Atomnaya Energiya, Vol. 10, No, 1, 
pp. 68-69, January, 1960 
Original article submitted May 17, 1960 


The fast neutron detector (graphite prism) described in [1] has been used to measure the absolute intensity of 
neutron sources, the average numbers v and 7 of neutrons produced in the fission by, and in elastic interactions of, 
14 Mev neutrons on U™ and U™ nuclei [2], and also to measure the cross-section for the (n, 2n) reaction in beryl- 
lium, lead and bismuth [3]. 


In 1954, this detector was used to measure V and 1 for 14 Mev neutrons interacting with u™ and Pu nuclei. 
These quantities were determined from the increase in the neutron current when 14 Mev neutrons are passed through 
a specimen of the fissionable element, 


Denoting the detector counting rate with and without the specimen by Ny and Np, respectively, and assuming 
that 


Ny 


we have the following formula for a spherical specimen with n nuclei per cm’: 


where 0j, is the inelastic cross-section, kg and ky is the sensitivity of the detector to 14 Mev neutrons and neutrons 
produced in the inelastic process, respectively, and 1 + 6 is a coefficient allowing for the increase in the number of 
neutrons as a result of fission by secondary neutrons, The first term on the righthand side represents the fraction of 
primary neutrons which have passed through the specimen without interacting, and the second term represents the 
fraction of neutrons produced as a result of the interaction between primary neutrons and the target nuclei, From 
Eq. (1) we have 


= ) 
a( 7 


k 1 
where a=-° 


ky 1 +8" 


Assuming that the basic processes were fission and the (n, 2n) reaction, we have 


vO; +26n; on 


Gin 


where 0 £,0p, 9p are the fission and (n, 2n) reaction cross sections, respectively, Hence 


o 
q=— 

|| 
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Spherical specimens made of metallic uranium and plutonium (outer diameter 4 cm, inner diameter 3,4 cm, diam- 
eter of aperture 2,1 cm) were placed in 0,9 mm thick cadmium jackets, These were then inserted into a thin-walled 
brass container which was located at the center of the graphite-prism cavity so that the center of the specimen 
coincided with the center of the target forming the source of neutrons, The neutron current was determined from the 
number of &-particles from the T(d, n) He‘ reaction. 


Results of measurements The quantity (1 + Y) was determined in a series of 
measurements with and without the specimens, In experi- 
Isotopes i+y ments without the specimens, the empty brass container 
with the cadmium cover was retained in the working posi- 
233 A554 | 4. tion, The table gives the average results of a large num- 
ber of measurements of (1 + Y). 


Pu2s9 
. In the determination of the coefficient Y we used 
Ko/ky = 0.874 # 0,030, as reported in [1], and an experi- 
mentally determined value of (1 + B), The latter was ob- 
tained with a source which imitated the spectrum of fission 
neutrons, a: comeotion was introduced = fission by slow neutrons, It was assumed that the correction coefficient for 
U5’ and Pu is the same as that for U (1,002 0, ag [2}), Assuming that o,,, is equal to the value of ojn ex- 


trapolated from the data given in [4] for u* and Pu " (O%n = 2.854 0,10 barn), we obtain the values of 1 given in 
the table, 


The fission cross section of U"** and Pu? is not very different from ojp. 


Assuming that ¢ = 0,2 0,1 barn for and 0,1 0,1 barn for and replacing in Eq, (3) 
by [1 —(ojn - 0 5)/o;,)", we obtain the values of v given in the table. 


Owing to the large excitation energy of the nuclei, the values of v for U, pu, u* and U™ [2, 5] in fission 
by 14 Mev neutrons are roughly equal, and according to our measurements approach the value 4.4 # 0,2, 


Our values of v are in good agreement with the results reported in [6-9]. 
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EFFICIENCY OF WAVEGUIDES USED AS ACCELERATING SYSTEMS 
IN ELECTRON SYNCHROTRONS 


A. N. Didenko and E£. S. Kovalenko 


Paper [1] is concerned with the possibility of using waveguides as the accelerating systems of electron synchro- 
trons, However, no estimate of the proposed system efficiency has been given. The following is an attempt to correct 


Let the wave LE, propagate at the azimuth angle ¢ in a closed waveguide with diaphragms. We are interested 
in its phase and group velocity, attenuation a, and the coupling impedance Rcoy. Together with the Q-factor, these 
parameters will determine the branch resistance Rpr, which represents the measure of the high-frequency system effi- 


Translated from Atomnaya Energiya, Vol. 10, No. 1, 
pp. 69-71, January, 1961 
Original article submitted August 31, 1960 


The exact calculation of the above quantities would 
be cumbersome. These calculations can be greatly simpli- 


Reoy Ohm: V feisty fied if we neglect the wave three-dimensional harmonics 
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Fig. 1. Dependence of a, Q,R 
the wavelength. 
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and rpy on 


Meg/cm and the effect of the curvature on the field structure in the 


waveguide cross section. The first of these assumptions is 
usually made in such calculations [2], while the validity of 
the second assumption follows from [1]; it is obvious that, 
for large synchrotrons,its effect is negligible. By calculating 
the magnitude of attenuation a, we obtain 


a. (1) 
Here, 


2WBox? (1 sh 2kyg-+-2kyg ) 
2 sh 2kyg+2kyg 


ch? kyg 
2wk3 sin? 


kBy = sh 2kyg—2kyg 


O is the conductivity of the walls; w is the angular frequency; 
A is the wavelength; v, h is the wavelength velocity and Dis 


the period of the system (the remaining data are given in 
Fig. 1). 


— 
0.18 
0 
é | 
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In a similar manner, we shall find the following expressions for the group velocity Vor and the coupling impe- 
dance 


2k2KBoc 
gr Ban? sh 2kyg--2kyg ch2kyg 
2A2 2 2 
ps )+ 2kyg — 2kyg Ach 2kyg — (2) 
where 
2k3 
A= “Ke sin® sin ) 
and 


Ky kx { 


coir “80 Bi sh 2kyg—2kyg 


(3) 


By using expressions (1) and (3), we readily find Q and rp; — the branch resistance per accelerator unit length: 
@ 4n Vor 


Figure 1 shows the dependence of &, Q, Reoy, Vor and 1, on the wavelength A for vp, = c, 2g= 6 cm, a= 
5.85 cm and D = A/4, In Fig. 2, Riou, Vorlc, and r, are given as functions of 2g for A = 9,3 cm, It follows from 
these diagrams that r,,, which characterizes the high-frequency system efficiency, slightly changes with changes in 
A, and that it varies considerably with changes in g: if 2g varies from 2 to 6 cm, ft, decreases from 0.5 to 0.1 meg/ 
cm, It is of interest to compare these theoretical 1, values with the branch resistances of actual linear accelerators, 
It is known that, for the Stanford 600-Mev accelerator, the interaction space of which is equal to 2 cm, ft}, = 0.47 
meg/cm. This confirms the fact that the theoretical t,,, Values which we obtained are within the limits of technical 
possibilities and that high branch resistance values are characteristic of high-frequency systems based on waveguides, 


Reoy OhM; tp, Meg/cm Vor Ic From Eq, (4), it follows that the expression for Rj, is given by 
06 Ry = ‘ (5) 
where Qegf is the effective waveguide quality factor, which takes into ac- 
1000 {2:2 0,10 count all the possible losses in it. 
; Let us compare the efficiency characteristics of a high-frequency 
ro. 4 system consisting of a closed waveguidewith diaphragms with the efficiency 
characteristics of the better and the most modern accelerating systems of 
+ 0.3 large synchrotons. With a sufficiently good manufacturing technology, it 
500 : 0,05 is possible to obtain a waveguide with the following approximate parame- 
ters: Rooy = 50 ohms, = 104, and = 0.lc. Under these conditions, 
J a waveguide that was used in such an accelerator as the German 6-Bev DESY 
would have the value Ry}, = 2000 meg for A = 10cm. In order to obtain 
04 such a branch resistance by means of a resonator system, it is necessary 
ie that each of the 80 resonators in DESY has the value Rp, = 25 meg. This 
| greatly exceeds the limits of technical possibilities (the DESY resonators 


0 2 4 6 g 


2g, om have Rp, = 2 meg). Thus, from the point of view of attainable branch 


resistances, a waveguide accelerating system offers greater possibilities 
Fig. 2. Dependence of Reoy, Vor and than a resonator system. 
2g. 
“br 00 “8 It follows from Eq. (5) that, all other conditions being equal, an in- 


crease in multiplicity m leads to an increase in the high-frequency system 
efficiency, and therefore, from the point of view of high-frequency power, greater m values are advantageous. How- 
ever, in comparison with ordinary resonator systems, the application of waveguide accelerating systems is unavoidably 
connected with a q-fold increase in multiplicity, In this case, due to radiation quantum fluctuations, the value of the 
root-mean-square spread with respect to phase will be ¥q times as large. In order to avoid this, if we wish that the 
particle losses remain unchanged, it is necessary to increase the high-frequency voltage amplitude by suitable means, 


By using the equation proposed by Christie and Robinson [3], it can be demonstrated that the particle losses for two 
accelerators are equal if 
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Fig. 3. Results of the numerical solution of 
Eq. (6). 


1 
(tg Psi— = (tg (6) 


It follows from Eq. (6) that, with an increase in multiplicity, 
the phase ¥, must increase to such an extent that Eq. (6) remains 
valid. Since the radiation losses are equal in both cases, this will 
result in a voltage increase given by 

Vp COS 

Vi COS 
However, with an increase in multiplicity, the value of Ror for 
system (4) will also increase, and therefore, if we take into account 
the radiation fluctuations, there is no evidence indicating that an 
increase in multiplicity would be undesirable. Figure 3 shows the 
results obtained by solving Eq. (6) numerically for different multi- 
plicity ratios m,/my. The region to the left of the straight line MN 
corresponds to phase values for which the transition to higher harmor 
ics leads to an energy gain. It is obvious that this region contains 
all the values of large accelerator equilibrium phases which are 
most interesting from the practical point of view. 


In conclusion, we hereby express our deep gratitude to Prof. A. A. Vorob'ev for his continued interest in the 


work and the discussion of the results. 
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TRANSFORMATION OF THE ENERGY OF SHORT-LIVED 
RADIOACTIVE ISOTOPES 


M. G. Mitel'man, R, S, Erofeev and N. D. Rozenblyum 


Translated from Atomnaya Energiya, Vol. 10, No. 1, 
pp. 72-73, January, 1961 
Original article submitted April 22, 1960 


The a- and 6-active isotopes that are formed as a result of the interaction between neutrons and the atoms of 
the material can serve as emitters of charged particles, which, as they accumulate on a collector, create a poten- 
tial difference. On the basis of this principle, it is possible to design a converter consisting of an emitter and a col- 
lector which are separated by a solid dielectric or vacuum. 


The number of charged particles that are produced by the emitter in unit time is proportional to the current 
obtainable in the converter: 


—0.693t 
), 


where Ng is Avogadro's number, o is the neutron capture cross section, n is the neutron flux, G is the emitter mass, M 
is the emitter material atomic weight, T is the evolving isotope half-life, and t is the emitter irradiation time.* 


For small half-life periods (t > T), thenumber of charged particles is independent of 
time. For large half-life periods (t < T), the number of charged particles is proportional to the 


irradiation time. 

Ni x In order to secure stable and time-independent operation of the converter, an emitter ma- 
N? 39 2 terial with the shortest possible half-life should be used. For a more efficient utilization of neu- 
N: EN trons, this material must have a large neutron capture cross section. 
\N: TS ad We performed experiments where a rhodium isotope (Rh'®*) with a capture cross section of 
N: \ 150 barn was used as the emitter. As a result of interaction between the neutrons and rhodium, 
Nz : \ isotope Rh'™ was formed, which emitted B-particles with an energy of 2.5 Mev and a half-life of 
41,8 sec, 


The diagram of the energy converter is shown in the figure. A rhodium wire with a diame- 
ter of 0.8mm and a weight of 0.42g, which was covered with insulating lacquer and a 1.5mm 
thick polyethylene layer, was placed in an aluminum shell which served as the electron collector. 
In assembling the converter, special attention was paid to the quality of insulation. 


Diagram of the 
energy conver- 
ter: 1) emitter; 
2) dielectric; 
3) collector. The thus prepared converter was placed in a channel of the research reactor at the I. V. 

Kurchatov Atomic Energy Institute, Academy of Sciences, USSR. An RK~1 television cable was 
used for the lead-out wires. The current and the voltage were measured by means of an electrostatic voltmeter and 
a string electrometer. 


By placing the converter in a neutron flux of 10” neutrons/cm*sec, a current of 4.2: 10 amp was obtained 


for an external resistance of 10!° ohms and a voltage of 420 v. 


The amount of 8-particlesthat are produced by the emitter for the above neutron flux can provide a current of 
6-107® amp. The observed lower current intensity in comparison with the maximum possible intensity can be ex- 
plained by the absorption of 6-particles in the emitter and the dielectric materials, 


After the converter was moved to a reactor zone where the neutron flux was 10°— 10"! neutrons/cm?sec, its 
current was reduced to 1.6. 107° amp in 2 min. This result does not contradict the theoretical assumptions. 


*R. Murray, Introduction to Nuclear Techniques [Russian translation] (IL, Moscow, 1955). 
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Converters of this type can be used for constructing high-voltage de voltage sources, Moreover, the converter 
can be used for determining the neutron flux and consequently, the reactor power level with respect to the intensity 
of the current supplied by the converter. 


The possibility of securing a steady-state current for a given neutron flux depends on the half-life of the isotope 
that is formed as a result of the emitter irradiation. 


In switching over from a lower neutron flux value to a higher one, accumulation of the radioactive isotope 
takes place. In practice, equilibrium is attained for t/T = 6.65 (deviation 1%), which for rhodium, corresponds to 
280 sec. A certain time is required for the transition from a higher neutron flux value to a lower one while the ac- 
cumulated radioactive isotope disintegrates, If the neutron flux value changes by a factor of 10, equilibrium is estab- 
lished almost for t/T = 10 (deviation 1%), which corresponds to 418 sec if rhodium is used. A device with such inert- 
ness can be used only in the case where the neutron flux changes slowly. 


In the case of fast-changing processes, we must determine the current derivative with respect to time, which 
is also proportional to the neutron flux: 


NqonG 0.693 _ 0.6931 
e 
For small t/T values (<0.0145), we can consider that 


dA NqoG0.693 
a = : n=const n, 


i,e., the derivative is independent of time. Thent <0.6 sec if rhodium is used. 


Thus, a converter with a rhodium emitter can be used for measurements that last either less than 0.6 sec or 
more than 418 sec, For a suitable emitter material, we can vary the inertness of the device and use it for processes 
that differ from each other with respect to their development rate and duration. For instance, for the Co™ (n, p) Fe® 
reaction, the B-active Fe™ isotope with the half-life T = 46 days is obtained; it can be used for measurements over 
periods t < 11 hr. 


In all the cases considered, the measuring circuit inertness was not taken into account. 


DATA ON THE SEPARATION OF BORON ISOTOPES 
IN THE FORM OF VOLATILE COMPOUNDS’ 


I. Kiss, I. Opauszky and L. Matush 


Translated from Atomnaya Energiya, Vol, 10, No. 1, 
pp. 73-75, January, 1961 
Original article submitted June 21, 1960 


Various methods developed for the separation of stable boron isotopes suffer from certain drawbacks from the 
practical point of view [1-7]. An investigation of the change in the isotopic composition of various boron compounds 
is therefore of theoretical and practical interest, Experiments with different boron compounds showed that isotopic 
effects arising in phase equilibria were contrary to classical ideas on these processes in many cases, Among the vola- 
tile compounds studied, the boron compound containing the heavy boron isotope was found to be more volatile than 
the corresponding compound with the light isotope. This phenomenon has been observed with other elements besides 
boron [8,9]. 


No general theory explaining this phenomenon has been developed up to now. However, it may be surmised 
that it is explained by both intramolecular and intermolecular interactions, depends on the association or dissociation 
of molecules in one or both paases, etc, [10], Due to the lack of a quantum mechanical theory for the liquid phase, 
certain difficulties arise in the theoretical solution of the problem and therefore it is necessary to supplement experi- 
mental data available at the present time by systematic investigations of various types of molecules, 


The purpose of the present work was a comparison of the practical efficiency of various oxygen-containing or- 
ganic compounds with BF, when used for enrichment of isotopes, and a study of the differences in nature and degree 
of the volatility of isotopic compounds of boron, 


Investigation of Complex Compounds of BF, 


Boron isotopes are separated most frequently by distillation of organic complexes of BFs and also by chemical 
exchange between a liquid complex and gaseous BFs. 


We investigated the systems anisole — BF, ethyl ether — BF, acetic acid — BF, ethyl acetate —BFs, etc. The 
boron isotope separation factor for the first two systems is known from [4, 5, 7, 8], The separation factor determined 
by a static method is not the only criterion of the applicability of a given isotope effect for enrichment in a certain 
system, and therefore experiments for comparing these compounds were carried out on the same enrichment column, 
1800 mm high and 12 mm in diameter, The column was packed with three-turn spirals of silver wire 0,2 mm indi- 
ameter. The length of the spirals was 1.3 mm, 


The complex of anisole with BF; decomposes to anisole and BF; when heated, Experiments were carried out at 
20°C with counter-current flow of the liquid complex and gaseous BF; liberated from it, Experiments with the other 
systems were carried out by the method in [4]. The absorption of neutrons was used to determine the isotopic compo- 
sition of the fractions obtained, The experimental apparatus and working procedure and also the neutron absorption 
method were described previosuly [6]. The experimental results are given in Table 1, 


The enrichment was calculated from the formula: 


Bu 
) at column head 


Aun—— 


Bu 
Bw i column still 


Table 1 shows that the complex of BF, with ethyl ether was enriched most, This compound is most stable, 
Though the separation factor in the system anisole— BF, is greater than in the system ethyl ether—BF;, the enrichment 


* article received from the Hungarian Peoples Republic. 
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in the first case was less than in the second under the given conditions, In addition, the complex of BF, with anisole 
decomposed considerably, The method of continuous purification of the anisole by distillation [11] complicates the 
technological process, in which the economy of the method must be considered, Though there was a considerable 
change in the isotopic composition of the other compounds during distillation, they were unsuitable for the enrich- 
ment of isotopes in practice because of irreversible decomposition of the complex. 


Table 1, Separation of Boron Isotopes in Various Complex Compounds 


Temperature, 
Compound Enrichment Comments 


CgHsOC H° BF 20 2.0 The complex formed tar after prolonged 
heating 


(CyH),0 BF; 60 2.4 


2H,0 * BF 80 1.9 The complex decomposed irreversibly 
during distillation 


CH; OH: BF 92 1,7 The complex decomposed irreversibly 
during distillation 


CH3sCOOH: BF; 96 1.9 Considerable decomposition of the complex 
was observed 


CH3yCOOC BF, 15 The complex decomposed on heating and 
the degree of enrichment was not deter- 
mined 


Investigation of Trialkyl Borates 


Boron halides and their volatile complexes have been used for the separation of boron isotopes up to now. 
Among the other volatile compounds of boron, it is profitable to investigate the alkyl esters of boric acid first of all, 
These compounds are stable and their boiling points are suitable for simple distillation. 


The distillation of methyl borate was proposed for the separation of boron isotopes in the USA in 1942 [12], but 
the separation factor found was considered low for the practical application of this method at that time, Since BF; is 
unsuitable for the working gas for the separation of boron isotopes in a thermal diffusion tube [13], attempts were 
made to use methyl borate for this purpose [14], The experiments showed [13] that this compound cannot be used 
for separation by thermal diffusion, but methyl borate was found to have a very high thermal stability. 


We investigated isotope effects in the liquid— vapor equilibrium with methyl, ethyl, and n-butyl borates, Ray- 
leigh distillation was used to determine the separation factor a [15], The experimental procedure was analogous to 
that described in [3]. The isotope effect of butyl borate could not be determined by Rayleigh distillation, For deter 
mining the separation factor, we used a fractionating column 12mm in diameter with a bed of packing 1000 mm high. 
The packing consisted of three-turn spirals 1.2-1.3mm long of silver wire 0.2mm in diameter, In the distillation of 
the complex of BF, with ethyl ether at 20mm Hg and 60°C with a liquid flow rate of 0.6 ml/min, the efficiency of 
the column corresponded to 25 theoretical plates with a separation factor of 1,028 [4]. The change in isotopic com- 
position during the distillation of methyl and ethyl borates was also determined under the previous conditions at 740 


mm Hg. a factor was calculated from the enrichment equivalent to 25 theoretical plates for each sys- 
tem: @ = 


The alkyl borates were prepared by the method described in [16], The azeotropic mixture of methyl borate and 
methanol was separated with the aid of CaCl, and LiCl. The physical constants of the materials corresponded to 
values given in the literature. 


The isotopic composition of the boron in the various fractions was determined on an MS-1305 mass spectrome- 
ter. The spectra of the alkyl borates were very complex, Peaks of individual ions were superposed on each other, 
for example, 42 (BY°O, + B!!OCHs) or 58 (B!°O3 +B"0,CH,), and the peaks of ions containing B without one proton 
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were superposed on peaks of ions containing B** etc, B’ ions were obtained with electrons with an energy of ~100ev, 
but they had a low intensity, Therefore. it was advantageous to obtain BF, from the esters and Goverinine the isotopic 
ratio in it, In this case, the ratio was readily determined by a comparison of the peaks from BF; or B* [17]. The error 
in the determination of the isotopic ratio B“/B"° was #0.01, 


The BF, was obtained from the boric esters in the following way. The ester was hydrolyzed with water and H,F, 
and NH,OH were added to the boric acid liberated. A solution with pH 8 was obtained and evaporated to dryness, The 


addition of H,SQ, to the compound (NH,),O(BFs) and heating the mixture liberated gaseous BF. These reactions pro- 
ceed almost to completion, 


The results of experiments with the three esters are given in Table 2, 


Table 2, Separation of Boron Isotopes in Trialkyl Borates 


— 


Pressure, Enrichment | Separationfactorby | Separation factor 


mm Hg Temp. , °C etermined on a 


ices oncolumn A |Rayleigh's method. | column with 25 plates 


B(OCH,)s 740 56 4.0725 10035-00005 4.0028 
B(OCsH,)s 740 56 1.0420 4.00230. 0008 1.0017 
B(OC4Hy)s 20 128 1.0170 1/0007 


The table shows that as in the distillation of boron halides, in the distillation of alkyl borates the distillate was 
enriched in the heavier isotope. This change in the isotopic composition might be expected from a comparison of 
the structure of boron halides and borates, The structures of the two types of compound are similar: the boron atom 
lies at the center of the plane of a right triangle formed by halogen or oxygen atoms, 


The separation factor falls with a decrease in the size of the molecules, It is interesting to compare this pheno- 
menon with an experiment in [18]. The separation factors Bp" in ether— BF, and thioether—BF, complexes hardly 
decreased or even increased with an increase in the molecular weight of the complex: 


Separation 
factor 


(CH,),0-BF3 1.025 
(CgH,).0- BF 1,027 
Cgl,OH-BF, 4.027 
CoH! ICH: BF; 1,030 
(C4Hy)2S- BF, 1.033 
BF, 1.036 
BF, ~1 . 040 


During the distillation of complex compounds of BF with ether there is chemical exchange between BF, bound 
in the complex and gaseous BFs passing into the vapor phase on dissociation of the complex. In this case the separa- 
tion factor is determined primarily by the exchange equilibrium and the reason for the isotopic effect we observed in 
the distillation of alkyl borates is the difference in the vapor pressure of the two different isotopic molecules, 


The degree of isotopic exchange during the distillation of methyl borate approximately equals that during the 
distillation of BFs and BCls. The advantages of esters are their thermal stability, low affinity for water, and the 
simplicity of their recovery. The distillation may be carried out at low temperature and atmospheric pressure, Boric 


esters and their decomposition products do not corrode glass or metal, Very diverse compounds of boron may be ob- 
tained in high yield from methyl borate. 


On the basis of all that has been stated above, it may be considered that boric esters are of practical use in the 
separation of boron isotopes, 
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EXPERIMENTAL DETERMINATION OF AXIAL SELF-ABSORPTION 
IN CYLINDRICAL Co®® SOURCES 


K. K, Aglintsev, G, P, Ostromukhova, and E, A, Khol'nova 


Translated from Atomnaya Energiya Vol, 10, No, 1, 
pp. 75-76, January, 1960 
Original article submitted August 18, 1960 


Co® sources are widely used in various branches of science and technology as sources of ¥ -radiation, and 
frequently have relatively large dimensions so that they cannot be looked upon as point sources, The dose rate due 
to such sources is usually calculated from the relationship 


p=4 Tye pir, (1) 


where [lo is the linear absorption coefficient of the y -radiation in air, r is the distance from the source in cm, A is 
the activity of the source in me, I is the Y-constant of the isotope which is numerically equal to the dose rate (r/hr) 
at a distance of 1 cm from a specimen having an activity of 1 mc, and f is the self-absorption coefficient which 
takes into account the absorption of ¥ -radiation in the source itself. The latter coefficient indicates the difference 


between the Y -radiation outside the source and the total y -radiation which is determined by the activity and the 
spectral characteristics of the isotope, 


Measured Values of the Self-absorption Coefficient f The correction for self-absorption of the radiation 


— leaving a cylindrical source along its axis is usually calcu- 
Length of Diameter of j lated (separately for each Y -line) from the formula [1] 
source, mm source, mm 

2 2 0.96 f =, (2) 
4 5 0.92 
ine : aa. where y is the linear absorption coefficient of the y -radia- 
15 5 0.826 tion of cobalt and J is the length of the source in cm, 
However, the calculation of f from Eq, (2) cannot be con- 
f, rel. units sidered as entirely correct since the derivation of this for- 
16 mula is based on the assumption that the radioactive iso- 
tope is uniformly distributed in the source [2] and that 
the radiation scattered within the source does not give 
rise to additional ionization. These two effects are 
R neglected in the calculation. 
0.9 


In order to improve the value of the correction for 
~ self-absorption, an attempt was made to determine f for 
Co® sources experimentally, This was done by measuring 


mae the dose rate and the absolute activity of sources of various 
08 N lengths (between 2 and 15 mm), The dose rate was meas- 
aS ured at the end of the specimen, using a narrow collimated 


beam and a standard apparatus incorporating a free-air 
chamber [3]. The absolute activity of the sources was 
07 measured with the aid of a differential y -calorimeter [4]. 
"0 5 1” 15 The results obtained are shown in the table and figure 
Length of source, mm (curve 1), Curve 2 shows the values of f based on Eq, (2). 


Dependence of the self-absorption coeffi- Comparison of the two curves shows that Eq, (2) is 
cient on the length of Co® sources; 1) ex: correct for sources whose length does not exceed 2-4 mm, 
perimental, 2) calculated. As the length increases, the discrepancy between the ex- 


perimental and calculated data also increases, This can apparently be explained by the effect of radiation scattered 


within the source [5, 6]. 
In order to avoid errors in calculations of the dose rate P, the coefficient f should be determined from the 
experimental curve (curve 1), 
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ATTENUATION OF GAMMA-RADIATION BY CONCRETE 
AND SOME NATURALLY OCCURRING MATERIALS 


P, V'yugov, K. S, Goncharov, S. Dementii, 
and A, M, Mandrichenko 


Translated from Atomnaya Energiya, Vol, 10, No, 1, 
pp. 76-79, January, 1960 
Original article submitted September 5, 1960 


The problem of shielding of the personnel and the general population from dangerous radiation is of great im- 
portance in the design of linear accelerators, Owing to the considerable length of these accelerators, the shielding 
structure is relatively expensive, and hence it has become customary to replace concrete by earth (for example, iri the 
case of the large linear electron accelerator at Orsey, France) or other less expensive materials, 


The shielding properties and the cost of various types of concrete are well-known [1], while for other building 
materials there is practically no published data [2], The present paper reports results of studies of the attenuation of 
Co” y -rays by materials whose chemical composition and density are shown in Table 1, 


TABLE 1, Characteristics of Shielding Materials 


Losses on Composition, 


Material | calcina- 
tion, / SiO0g Al2O Ti0g Fe203 CaO 


Earth 78.85 11.94 
Sand 96.16 1.48 
Clay 64.40 | 12.92 | 0.90 


*Refers to freshly deposited materials, 


Design of the experiment. The geometry of the experiment is illustrated by Fig, 1. A square-shaped aperture 
(1.04 m‘) in a concrete wall 2m thick was employed. A spherically symmetric Co™ source having an activity of 
0.57 C was placed on the axis of the aperture and in the plane of the wall, The detectors, which were in the form 
of stilbene crystals mounted on FEU-19M photomultipliers, were located on the axis of the aperture at 2 and 4 m from 
the source, respectively, All the measurements were thus carried out under wide-beam conditions, 


Tin boxes and concrete blocks 1 X 1 X 0,1 min size were employed, In the case of loose materials, the boxes 
were inserted into the channel in the wall, and the absorption was determined with and without the materials so that the 
effect due to the wall of the box could be eliminated. The gaps between the boxes and the walls of the aperture were 
filled with the material under investigation, No attempt was made to ram the material into the boxes, The thick- 
ness of the absorber was increased in the direction towards the source until the indications of the detector placed at 


2 m from the source were of the order of the background, which corresponded to 120 cm of concrete and 170 cm for 
the other materials, 


The results obtained with the detector at 2 m from the source are shown in Fig, 2, A similar graph was ob- 
tained with the detector at 4 m from the source, For comparison, Fig, 2 also shows the attenuation curve for 1,25 
Mev 7 -rays by concrete having a density of 2.3 g/cm’, This curve was plotted on the basis of the data given in [3], 


The attenuation curves can be used to calculate the absolute absorption coefficients of freshly deposited 


materials for Co Y-rays. For absorber thicknesses in excess of 0,3-0.4 m, the attenuation coefficients remain prac- 
tically constant. 


: 

Alka- 

3.86 1.54 1.00 | 2.42 0.25 1.46 

z 0.50 0.46 0.10 | 0.50 0.20 1.64 

Z 5.00 4.30 1.48 | 2.90 0.22 1.39 
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Concrete blocks or 
boxes filled with 
shielding material 


Detector 


Fig. 1. Experimental arrangement, 


The data now reported can be used to estimate 
the thickness of earth, sand or clay which will be as 
effective as concrete in reducing the level of radiation 
to a given figure. Thus, for example, 1 m of concrete 
is equivalent to 1,36 m of sand, 1,52 m of earth, and 
1,61 m of clay. 


Discussion of results. It is known that the y-ray ab- 
sorption coefficients for Co™ y ~rays are proportional to 
the atomic number Z, while for y -rays in the region of 15 
Mev they are proportional to z*. For concrete and the 
other materials investigated, Zef¢ is very clearly constant. 


The energy spectrum of the lateral radiation 
from linear accelerators has not as yet been investi- 
gated although it is to be expected that relatively 
low-energy Y -rays (10 Mev or less) should be present 
at large angles to the beam. 


TABLE 2, Calculations of the Cost per 1 m’ of Shielding 


6.0 


Q 


~ 
3 
4.0 
3 
E 30 
oO 
x 
10 


50 100 150 
Thickness of absorber, cm 


200 


Fig, 2, Attenuation curves for Co Y -rays for concrete, 
sand, earth and clay. Source to detector distance 2 m, 
@ — Concrete, — coarse-grained sand, O —-earth, 
clay,- = — curve computed on the basis of data given 
in [3]. 


Character and construction of shielding 


Cost per 1 m’ of shielding 
equivalent to 5 m of 
concrete 


Cost per 1 m? of shielding 
equivalent to 2 m of con- 
crete 


rubles %o rubles %o 


Solid continuous concrete wall... 
Two reinforced -concrete walls separated by sand*..... 
Reinforced -concrete supporting wall with a bank of earth 
on one side covered by turf: 
to a height of 5m.. 
to a height of 7 m. 


*3 km from the pit. 


with all the other cases, 


Two reinforced-concrete walls separated by ordinary earth. . 


Shielding wall made of M-25 earth blocks** .........-. 


515 1285 
315 61 425 33 
ee 320 62 435 34 
ee 250 49 320 25 
ee 260 90 335 26 


* * The attenuation coefficient for earth blocks was assumed to be the same as for freshly deposited soil, which in the 
present case leads to an overestimate of the cost, The building of the shielding wall from earth blocks ensures com- 
plete industrialization of the operation, Constructional difficulties are reduced by a factor of 2,3-2,5 as compared 


\\ | 

Source 

YY fy | 

Y Yj 

| 375 900 70 


TABLE 3, Estimated Cost of the Shielding for a 2 Bev Linear Accelerator 


Cost of shield- | Consumption of | Consumption 


Character and design of the shielding ing, % % 


— 


Solid continuous concrete walls with prefabricated reinforced- 
concrete covering of the accelerator tunnel .........6. 100 100 100 
Inner (between the accelerator tunnel and working space) 
screening wall made in the form of two reinforced-concrete 
walls separated by earth; the outer wall made in the form 
of a reinforced -concrete support covered from outside by a 
earth-bank and turf; accelerator tunnel covered with M-50 
earth-blocks on a prefabricated reinforced -concrete plate . 47 35 58 
The same, with sand replacing the earth*®,.... ee 47 35 58 
Shielding made up of soil blocks M-25; accelerator tunnel 
covered by M-50 soil blocks on a prefabricated reinforced - 


*3 km from the pit. 
** The attenuation coefficient for earth blocks was assumed to be the same as for the freshly deposited soil, which 
in the present case leads to an overestimate of the cost, The building of the shielding wall from earth blocks ensures 
complete industrialization of the operation, Constructional difficulties are reduced by a factor of 2,3-2,5 as com- 
pared with all the other cases, 


In the present measurements the shielding materials were not compressed since this would only improve their 
shielding properties, It follows that the difference between the attenuation of Y-rays by concrete, sand and earth at 
energies corresponding to the Co” y -rays should not be appreciably greater when the shielding concrete walls are 
replaced by walls made of sand and earth. 


Such a replacement should improve the neutron shielding properties of the walls, owing to the lower Zo¢¢ and 
greater thickness and moisture content, 


Costing problems, An estimate was made of the change in the cost of the shielding when the concrete walls of 
a linear accelerator are replaced by sand and earth in the middle parts of the USSR, The cost of sand and earth in- 
cludes the cost of ramming [4]. 


Table 2 gives the estimated cost per 1 m? of the shielding for different materials, Table 3 gives the estimated 
cost for a 2 Bev linear accelerator. 


As can be seen from Tables ‘2 and 3, replacement of the concrete by earth and sand leads to a considerable 
reduction in the cost of the shielding of a linear accelerator, 


In conclusion the authors wish to express their gratitude to V. V. Katrich and V. S, Poryatui for assistance in 
this work, 
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NEWS OF SCIENCE AND TECHNOLOGY 


TENTH INTERNATIONAL CONFERENCE ON HIGH-ENERGY PHYSICS 


A. A. Tyapkin 


Translated from Atomnaya Energiya, Vol. 10, No. 1, 
pp. 80-83, January, 1961 


The Tenth International Conference on High-Energy Physics was held from August 25 to September 1, 1960 in 
Rochester (USA), Three hundred fifty physicists from 29 countries took part in the Conference. 


On the first days of the Conference, simultaneous sessions were conducted on four separate branches of high- 
energy physics, In these sectional sessions new experimental and theoretical results obtained in various laboratories 
during the preceding year were reported and discussed in detail, In plenary sessions, reports on specific branches of 
high-energy and elementary-particle physics were given, In these reports the results presented in the Conference 
were generalized and systematized, and on this basis, a general view was given of the present state of problems in 
the respective branches of high-energy and elementary-particle physics. 


In one of the sections of the Conference, experimental investigations of strong interactions of nucleons and 
pions were discussed, Material on this subject was summarized in the papers of J. Ashkin (Pittsburgh, USA) and O, 
Chamberlain (Berkeley, USA). 


In certain laboratories (Dubna, Berkely, CERN) work has been carried out concerning the precision with which 
the law of conservation of isotopic spin is satisfied, The most precise data have been obtained at Dubna in studying 
the forbidden reaction in which a neutral pion and an alpha particle are created in d-d collisions, The cross section 
of the reaction was found to be less than 2- 107 cm? and practically equal to the cross section of the electromagnet- 
ic procesd+d—> + y, 


The first experimental information was obtained on the very important problem of pion-pion interactions, The 
group working at Berkely with the large liquid-hydrogen chamber according to the method proposed by Chew and 
Low obtained data on the cross section of the ™1° interaction. For relative meson energies of three and five pion 
masses, the cross sections obtained were equal to ~150 and ~300 mb, respectively, 


Another group at Berkely studied the momentum spectrum of He® from the reactions 
He3-+- 19 
pt+d-j{ He®+n---n* 
70, 


Very interesting preliminary data were reported by Crowe, In the high-energy part of the spectrum a sharp maxi- 
mum with a width close to that of the apparatus was found, If the data on the maximum and its width are con- 
firmed by further measurements, this will indicate the existence of a new particle with a mass of about 2.3 pion 
masses, 


+ 1° at zero 


+ 
A value of (4-8) mb was obtained at Dubna for the cross section of the process 1 
energy of interacting m ~-mesons, 


The most precise data on the scattering of pions by nucleons was obtained this year at Saclay (France), In 
particular, for the scattering cross section of negative pions by hydrogen, the values of 45.8+1.8 and 58,0 # 1.8 mb 
were obtained as maxima at 600 Mev and 900 Mev, respectively. R, Glasser, et al. (Washington) carried out direct 
measurements of the lifetime of the neutral pion. The measurements were performed in photoemulsions by observ- 
ing Dalitz pairs from ™°-mesons created together with *-mesons in the decay of positive K-mesons. The mean 


life of mesons was found to be (2.3 #0.8) - 107° sec, 


The laboratories at Rochester, Dubna, Harwell, and Harvard presented new data on elastic scattering of nucle- 
ons by nucleons which were carried out with the purpose of establishing the scattering amplitude, 
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In Rochester, the depolarization factor D for angles of 30° and 60° in the center of mass system (CMS) andthe 
triple scattering parameters R and A for seven angles in the interval from 30° to 90° (CMS) were measured for p-p 
scattering with 210-Mev protons, For p-n scattering, measurements of polarization were performed for 12 angles in 
the interval from 30° to 90° (CMS). 


The correlation coefficient of polarization C, (90°) in p-p scattering was measured at Dubna for proton ener- 
PP 
gies of 315 and 630 Mev. 


The results of phase analysis of p-p scattering at energies of 95, 150 and 310 Mev were reported at the Confer- 
ence by the Moscow group. 


Several experiments on elastic scattering of nucleons by nucleons were performed at Liverpool, Berkely and 
Dubna in order to determine the constants of pion-nucleon interaction. 


The Conference gave its approval to the program of the Photomeson Laboratory of the P, N, Lebedev Institute 
of Physics, Academy of Sciences, USSR (FIAN). This program, which was presented in the paper of A. M, Baldin, 
involves the analysis of experimental data on the photoproduction of mesons in the near-threshold region, 


In the field of photoproduction of mesons, a great deal of research has been done on the second and third reso- 
nances of the meson-nucleon interaction, Measurement of proton polarization in the reaction y +p—~p+ n° was 
done in Frascati (Italy). The differential cross section for the photoproduction of positive pions in the resonance 
region was studied at Frascata and at the California Institute of Technology (USA), The data obtained agree well 
with the assumption that the second resonance is due to a dipole transition in the d3/z ~state and the third resonance, 


to an electric quadrupole transition in the f5/,-state. Both resonances correspond to the isotopic spin state T = %,, 


Some work on the strong interaction of nucleons and pions was performed at energies greater than 1 Bev. 
Measurements of proton polarization elastic p-p and p-Be collisions were made with the accelerator at Saclay for 
proton energies of 1,7 Bev. 


The double scattering of antiprotons with momenta 1,5-2,0 Bev/c was investigated by the group at Berkely 
with the 72-in liquid hydrogen chamber. The magnitude of antiproton polarization averaged over the angular range 
6° - 25° (CMS) was found to be 0,51 #0,09 under the assumption that the antiproton magnetic moment is equal in 
magnitude and opposite in sign to that of the proton, Study of the annihilation of antiprotons in the chamber showed 
that the fraction of annihilations leading to two pions or to two K-mesons does not exceed 1/400 and 1/1000, respec- 
tively, of the total probability of annihilation, 


The scattering of antiprotons by protons and neutrons at an energy of 1.2 Bev was investigated by Wenzel's 
group (Berkeley), In this work, the cross sections (elastic and inelastic) for these collisions were shown to be equal. 


The reports of the Joint Institute for Nuclear Research (OIYal)(V.Petrzhilka, V. I. Veksler) contained a large 
amount of material on the study of pion and nucleon interactions in the energy region from 3 to 9 Bev, From a com- 
parison of experimental data obtained on the synchrophasotron with theorétical investigations at OIYal(D. I. Blokhintsev 
et al.) and FIAN (D. S, Chernavskii) it was concluded that the interaction of nucleons and pions occurs mainly in the 
pion shells of nucleons, 


Using the large CERN accelerator, von Dardel's group measured the total cross section for the interaction of 
pions, protons, and antiprotons with protons in the energy range from 3 to 10 Bev, The pion-proton cross sections in 
this range have values of the order of 30 mb, The proton-proton cross sections decrease with increasing energy, haw 
ing at 10 Bev a value of about 42 mb, The cross section for p-p scattering at 10 Bev exceeds the cross section for p-p 
scattering by 20 mb, 


A separate plenary session was devoted to the structure of elementary particles, New experimental data onthe 
electromagnetic structure of nucleons were reported in the papers of R, Hofstadter (Stanford, USA) and R. Wilson 
(Cornell, USA), New measurements of form factors were made for q’ = 25 (in reciprocal fermis), The difference 
between magnetic and electric form factors for the proton was established. With increasing energy, the magnetic 
form factor decreases and the electric form factor approaches a constant value, In all earlier work the equality of 
these form factors had been assumed, The experimental data on electron scattering by the neutron continue to be 
compatible with the assumption that its electric radius is zero, 


A great deal of experimental material on strange particles was presented at the Conference, In the section 
dealing with this subject, 27 papers were given, 
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The paper of E, Turkot (Cornell University) was devoted to the photoproduction of strange particles, Study of 
the reaction y + p— A’+ K* at energies of 910-1040 Mev showed that the differential cross section for this reaction 
is isotropic in the CMS and amounts to about 1° 107" cm?/sr,_ In an investigation of the reaction y +p > 2°+ Kt 
it was found that the excitation curve is proportional to the cube of the momentum of the K' -meson, The photo- 
production of K* -mesons in deuterium for two values of the y-quantum was also studied, It was found that within 


the limits of error of measurement the cross section for production of K* -mesons with a neutron is equal to that with 
a proton, 


The group working with the 72-in liquid-hydrogen chamber at Berkeley studied the generation of strange par- 
ticles in ™” -p collisions in the near-threshold region, In this work, in order to determine the relative frequency of 
hyperons, the peculiarities in behavior of different terms in the expression for the differential cross section for the 
reaction t+ + p— A° + K° were investigated in their dependence on the threshold energy of the reactions: 


K*. 

M. I, Solov'ev reported on new data (obtained by the Dubna group with the propane chamber) on the generation 
of strange particles in pion-nucleon collisions at an energy of 8 Bev. I. V. Chuvilo (Dubna) reported on the generation 
of strange particlesby 3-Bev negative pions. Considerable interest was aroused by data obtained here on an asymme- 
try in the decay of A-hyperons, which indicated the possibility of nonconservation of parity in strong interactions, 


Investigations carried out at Dubna showed that the picture presented at the Kiev Conference of the generation 
of unstable particles by 7-Bev ™-mesons can be explained fully without assuming the existence of a new particle, 


The first results on the creation of strange particles by 16-Bev ™ -mesons were presented at the Conference 
(obtained at CERN with the 32-cm hydrogen bubble chamber), In this work, a value of ~1 mb was obtained for the 
cross section for generation of A-hyperons and ~3 mb for K’-mesons, In the CMS the A-hyperons created are 


strongly directed backward, This result indicated that in the 16-Bev region the same laws hold as at energies inves- 
tigated earlier at Dubna, 


In the propane chamber at Saclay, 81 cases of the production of =*-K* pairs by hydrogen were obtained in a 
beam of positive pions with a momentum of 1.15 Bew/c, The cross section for this reaction was found to be 0,26 
£0.05 mb, The resulting angular distribution of particles produced in this reaction does not agree with the distribu- 
tion required by isotopic invariance, However, the statistical accuracy does not permit a final conclusion. 


M. Good (University of Wisconsin, USA) reported on research on the interactions of 1,15-Bev/c K~-mesons in 
a hydrogen bubble chamber. In this work the analysis of momentum spectra of #*-mesons produced in the reactions 


K" + p-»A + nm suggests the existence of excited states of hyperons, 


The experimental work on weak interactions presented at the Conference was concerned mainly with research 
on the properties of 1 -mesons, 


V. A, Lyubimov (Institute of Theoretical and Experimental Physics, Academy of Sciences, USSR) reported on 
measurements of the helicity of 1-mesons, It was found that in n-» » decay the negative p-meson has a right-hand 
helicity (P. = +0,65 £0,22) and the positive one has left-hand helicity = -0,15 + 0,21), 


The first results were obtained in experiments to measure the difference in precession and cyclotron frequen- 
cies for -mesons, At Columbia University it was found that 


a 
—0.14 / 


According to preliminary data obtained at CERN, this quantity equals 
g—2 a 
(1.8 + 0.3). 
At Berkeley, experiments were conducted to determine the probability of a neutrinoless transition from  - 


meson to electron, It was found that the probability of this transition amounts to less than 4. 10°° of the probabil- 
ity of the neutrino transition, 


A, L, Mukhin reported on research carried out with the synchrocyclotron at Dubna on the intensity of radiation- 
less transitions in the mesic atoms of Pb, Bi, Th, U"™ and U . The experimental results of this work are satisfacto- 
rily explained by assuming the presence of a mechanism of energy radiation analogous to Raman scattering in optics, 


The asymmetry in angular distribution of neutrons emitted by nuclei of Mg and S in the capture of polarized y™- 
mesons was subjected to several investigations, In Chicago, measurements were made of the probability of 7 -meson 
capture in the isotopes ci" and Cl®® The relative probability of these captures was found to be 0,694*0,034, The 
ratio calculated by the formula of Primakov is equal to 0,82. At Columbia University, work was successfully complet- 
ed on the observation of muonium (yte™), It was found that the degree of polarization of muonium in argon depends 
strongly on admixtures of other gases, For the first time, scattering of 1-mesons was studied with an accelerator, Mu 
mesons with a momentum of 1,95 Bewe obtained with the Bevatron (Berkeley) were scattered by lead. Anomolous 
scattering was not observed, At Liverpool, measurements were made of the mean life of the 1-meson in connection 


with the discrepancy of results obtained earlier at Chicago and at CERN. The value found was T= 2,230 *0,006 


A large number of theoretical papers were presented at the Conference, A great deal of interest was elicited 
by the reports of R, Eden (Berkeley) and J, Polkinghorne (Cambridge University, Great Britain), which were devoted 
to proof of the double dispersion relations of Mandelstam, An important contribution to the study of the analytic 
structure of matrix elements was made in the report of A, A. Logunov (Dubna), The papers of A, P. Rudik (Institute 
of Theoretical and Experimental Physics) on singularities of scattering amplitude in perturbation theory should be mer 
tioned, as should his reports on the work of V. N. Gribov (LFT), V. B, Berestetskii and I, Ya, Pomeranchuk (ITEP) on 
the asymptotic values of cross section in the very high energy region. 


The following papers were concerned with practical application of the dispersion relations; G, Chew (Berkeley) 
on pion-pion scattering, M, Cini (Italy) on new numerical results obtained at CERN on pion-nucleon scattering, and 
A, M. Baldin (FIAN, USSR) on a new treatment of data on photoproduction based on the relativistic dispersion rela - 
tions. The report of D, V. Shirkov (Dubna) on a new method of obtaining equations for partial waves contained a 
series of numerical data on m-m and m -N scattering at threshold, 


Problems of elementary-particle theory were discussed at a special session presided over by W. Heisenberg 
(Munich), Considerable interest at this session was aroused by the paper of Y. Nambu (Chicago), in which the math- 
ematical methods developed by the Soviet scientist N. N, Bogulyubov in the theory of superconductivity was success- 
fully applied to field theory. Bogolyubov reported on the interesting results of the work of A. N, Tavkhelidze (Dubna), 
in which the methods of superconductivity theory were applied to field theory. 


In the paper of Heisenberg, some general properties of equations obtained earlier by him were investigated, 


D. I. Blokhintsev (Dubna) reported on attempts, based on the mathematical work of Lapo- Danilevskii, to solve 
some problems of elementary-particle theory without the use of the expansion in the interaction constant, 


Several papers were concerned with semiphenomenological theories, S, Drell (Stanford, USA) showed that at 
high energies y-quanta give a more effective method of obtaining high-energy pions than do nucleon-nucleon colli- 
sions, which lead to a high multiplicity of the mesons produced, 


The papers of M. A, Markov (Dubna) and T, Lee were devoted to processes involving high-energy neutrinos, 
These reports discussed the question of the possible existence of intermediate mesons in weak interactions and also 
the question, raised earlier by B, Pontecorvo, of the possible nonidentity of neutrinos involved in 6 -processes and 


neutrinos formed in the decay of mesons, The papers of R, Feynman and M, Gell-Mann (Cal, Tech.) discussed prob- 
lems associated with the conservation of vector current, 


The Conference, which proceeded in a friendly and business-like atmosphere, certainly played a significant 
role in uniting the efforts of scientists of various nations who are occupied with the study of the most important field 


of contemporary physics, The Proceedings of the Conference will be published by the University of Rochester early 
in 1961, 


After the Conference the Soviet delegation visited and inspected many American laboratories which are doing 
reseach in the field of high-energy physics, The delegates became familiar with the 240-Mev synchrocyclotron at 
the University of Rochester; the 1,2 Bev, strong-focusing electron accelerator at Cornell University; the 380-Mev syn- 
chrocyclotron at Columbia University; the 3- and 30-Bev proton accelerators at the Brookhaven Laboratory; the 3-Bev, 
high-precision proton accelerator with weak focusing, under construction at Princeton; the 6-Bev, strong-focusing 
electron accelerator under construction in Boston; and also with the laboratories of the Massachusetts Institute of Tech- 
nology. The Soviet delegates taking part in the International Conference on Instrumentation for High-Energy Physics 


became familiar also with the 960-Mev linear accelerator at Stanford University and with the 730-Mev synchrocyclo- 
tron and 6,3 Bev Bevatron at Berkeley. 
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The delegates from the Soviet Union were given a warm reception everywhere, both from the officials who con- 
ducted this Conference and from the American physicists, It should be mentioned also that there in a great deal of 
interest in the USA in Soviet physics, At the request of American physicists, D, I. Blokhintsev, V. I. Veksler, V. P. 
Dzhelepov, S. Ya, Nikitin, and A, S, Belousov gave survey reports at Brookhaven and at Stanford on research which 


has been done in the Joint Institute for Nuclear Research, the Institute of Theoretical and Experimental Physics, and 
the P, N, Lebedev Institute of Physics, Academy of Sciences, USSR, 


At the meeting of the Committee on High-Energy Physics of the Union of Pure and Applied Physics, it was de- 
cided to increase the interval between conferences to two years in order to permit more complete investigations, 


The Eleventh International Conference on High-Energy Physics and the Conference on Instrumentation is 
planned for July, 1962 in Geneva, The next conference on accelerators is planned for September, 1961 at Brookhaven, 


INTERNATIONAL CONFERENCE ON INSTRUMENTATION 
FOR HIGH-ENERGY PHYSICS 


A. A, Tyapkin 


/ 
Translated from Atomnaya Energiya, Vol. 10, No. 1, 
pp. 83-85, January, 1961 


The International Conference on Instrumentation for High-Energy Physics took place in Berkeley (USA) from 
September 12 to 14, 1960, Two hundred forty persons from 13 countries took part in the Conference, 


Seventy-three papers on new developments of methodology were given. The Soviet delegation presented eight 
papers covering part of the important work which has been done in the field of high-energy physics, 


The first session was devoted to problems of separation and focusing of high-energy particles, J, Murray 
(Berkeley) reported important results in the use of glass electrodes in electrostatic separators, Tests showed that the 
use of glass plates made it possible to double and even triplethe regularly used electric field strength. At Berkeley, 
separators with glass plates have already been used with fields of 150 kv/cm instead of the usual 60 kv/cm., 


In the paper of M. Good (University of Wisconsin) the problem of an energy limit for electromagnetic separa- 
tors was discussed, In the speaker's estimation this limit for the separation of K-mesons lies at about 10 Bev/c, 


Several reports were devoted to the production of large magnetic field. H. Furth (Berkeley) described impulse 
coils for obtaining magnetic fields up to 160 kilogauss in small volumes, The report of R. Post (Berkeley), which was 
received with great interest, presented a detailed discussion and also experimental data on obtaining constant, high- 
intensity magnetic fields in large volumes using sodium or aluminum coils cooled with helium. The electrical resis- 
tivity of these metals over a wide range is independent of the magnetic field strength, As an example, the speaker 


presented calculations of a coil, cooled by helium to 10°K, for producing a constant field of 60 kilogauss in a volume 
of 1 m*, 


A large amount of material on electronic methods of recording particles was presented at the Conference, A 
central place in this subject was set aside for new detectors of charged-particle tracks spark, discharge, and lumi- 
nescent chambers, The survey paper of M. S, Kozodaev on spark chambers developed in the Soviet Union was met 
with great interest, A, A. Tyapkin reported on research atOIYaland MIFlonthe conditions under which the discharge 
develops along the particle track in the chamber, 


An interesting report was given by S, Fukui (Japan) on the results of tests of a high-frequency discharge cham- 
ber, The 150x 100x 20mm chamber was filled to atmospheric pressure and neon and a half-percent admixture of 
argon, A 3.7-kv/cm, 2.8X 10°-cycle field from a magnetron is applied for about 107" sec with a delay of 0.5 psec 
after passage of the charged particles through the chamber, Along the track of the particle, a luminosity of the gas 
arises which diminishes more slowly than in the case of a spark chamber. 


The papers of J, Cronin (Princeton) and B. Cork (Berkeley) attested to the successes achieved in the develop- 


ment of spark chambers, The Princeton chamber contains 19 aluminum electrode-plates about 6 mm thick, situated 
10mm one from another and having a sensitive area of 18x 18 cm®, The chamber is filled with neon. This chamber 
has already been used for measuring the proton polarization, Figure 1 shows a photograph of two projections of the 
track of a proton which has undergone scattering in one of the plates of the chamber, 


An analogous spark chamber was developed at Berkeley and was tested with both neon and argon, An efficiency 
of 99% was attained with pulses of 10 kv amplitude in the case of neon and 15 kv in the case of argon (the interelec- 
trode spacing in this chamber was about 6mm), The recovery time of the chamber after a discharge is about 10 msec, 
The sensitivity time of the chamber without a clearing field is 10 psec, and with a field of 160 v/cm it decreases to 
0.5 usec, Interesting results were obtained when the chamber was tested in a magnetic field of 13,000 gauss parallel 
to the plane of the electrodes, In conclusion, the speaker showed electrodes made of graphite for a new chamber to 
be used in measuring the polarization of protons, D, Meyer (Michigan, USA) reported on a discharge chamber whose 
electrodes were made of a thin foil. A chamber with such an electrode system is useful in the study of the angular 
distribution of particles liberated from a liquid hydrogen target. 


Fig. 1. Projections of the spark track of a proton 


Fig. 2, Track of a cosmic ray particle 


Considerable success was achieved in several laborator- 
ies in the USA during the past year in the development of 
luminescent chambers, which were proposed earlier by E, K. 
Zavoiskii (USSR), Extensive work on the creation of scintilla- 
tion chambers became possible thanks to the appearance in 
laboratories of reliably effective electron-optical converters, 
Each stage of an electron-optical converter insures a 30-fold 
amplification for a voltage of 15 kv, Three such stages are 
usually used in amplifying light from a scintillator, 


A paper presented by the Berkeley group reported on 
luminescent chambers intended for investigation of elastic 
and inelastic scattering of negative pions by hydrogen in the 
energy range 1,5-2,5 Bev, A crystal of sodium iodide is used 
as the scintillator, In studying inelastic scattering the use of 
two luminescent chambers is recommended, 


Lande (University of Pennsylvania) reported on a cham- 
ber whose scintillator is built up of fine fibers made of scintil- 
lating plastic, Total internal reflection in the fibers insures 

an effective collection of light at the photocathode of an 
electron-optical converter, This type of scintillator has been 
chosen by several laboratories in the USA, in particular at 
Princeton University and Massachusetts Institute of Technology. 
The construction of the chambers is practically identical, The 
diameter of the scintillating fibers is 0.5 and 0.7mm, The 
fibers are arranged in uniform rows in two mutually perpendi- 
cular directions and are fastened in a single unit with the aid 
of Epoxy resin, An independent light-amplification system is 
used for the fibers of each direction. The first stage of the 
electron-optical converter insures image reproduction on a 
fluorescent screen with an emission time of about 1 psec, The 
control system for this time includes the second stage of am- 
plification, For convenience of observation and photography, 
after the third stage of amplification, the image is focused on 
the photocathode of a television kinescope. Dr, Hill at MIT 
created a special scintillation chamber of dimensions of 12x 
12x 24cm which can be operated in a magnetic field of 50 
kilogauss for measuring the magnetic moment of the A° -par- 
ticle. A photograph of the track of a cosmic ray particle pass- 
ing through a fiber scintillator is presented in Fig. 2, 
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A great deal of attention at the Conference was paid to gas Cerenkov counters, The most important paper in 
this field was given by A, Roberts (Argonne National Laboratory, USA), who told of a new type of Cerenkov counter 
which permits precise measurement of the speed of a particle within a wide range of speeds and also determination 
of the direction of motion of the particle with high precision, In this detector, the Cerenkov light from the gas radi- 
ator is amplified using four photomultipliers, and the speed of the particle is determined from the diameter of the 
illuminated ring on a television screen, The direction of motion of the particle is determined by the position of the 
center of the ring. The relative precision of the velocity measurement is estimated to be ~10™%10"4, 


The paper of R, Schluter (Argonne National Laboratory) was devoted to the same problem of combining preci- 
sion of velocity measurement with wide measureable angle, In this detector, the Cerenkov light from the separate 
angular intervals is collected by fiber light-guides at the photocathodes of separate photomultipliers, 


I. V. Chuvilo gave a report on the gas Cerenkov counter used for the separation of K*-mesons at the synchro- 
phasotron of the Joint Institute for Nuclear Research, 


New ionized-particle detectors based on p-n junctions in semiconductors are of significance. The paper of L, 
Yuan (Brookhaven National Laboratory), which was received with much interest at the Conference, discussed ioniza- 
tion losses in the sensitive layer of a silicon detector. This work showed the possibility of clearly separating positive 
pions of momentum 700 Mewc from protons of the same momentum on the basis of amplitude analysis, The energy of 
liberation for pions and protons of this momentum amounts to 110 and 200 kev, respectively, 


C, Wiegand reported on a very complete multichannel system of scintillation counters developed at Berkeley, 
the impulses from whose photomultipliers are fed directly to the memory of a computer. This apparatus is especially 
constructed to obtain information on the m= interaction on the basis of a study of the processes 1 + p> 1” + m*4n 
and nt + p-~*t+4n*t+n, The determination of neutron energies up to 50 Mev is accomplished by time-of-flight 
measurements, In the future, similar systems of scintillation counters directly connected to computers will undoubt- 


edly replace the usual counting systems utilized in coincidence and anticoincidence circuits, 


At one of the sessions, proposed neutrino experiments were discussed in the presentations of C, Cowan (Wash- 
ington), K, Greisen (Cornell University), L, Lederman (Columbia University), and others, Theoretical problems re- 
quiring experimental tests with neutrinos were reported by C, Yang (Princeton), D, L Blokhintsev told about work in 
this field that has been carried out in the Soviet Union, Particle detectors planned for neutrino experiments have 
grandiose dimensions, For example, at Columbia University, a spark chamber of dimensions 2X 2X 1m is being de- 
veloped for the study of scattering of neutrinos by electrons, The electrodes of this chamber will be made of iron one 
cm thick, and the total weight of this system will exceed 15 tons, 


Specialists in bubble chamber technique discussed experiments in two separate sectional meetings, At these 
sessions, reports from the Radiation Laboratory in Berkeley were presented on a xenon bubble chamber, on experi- 
ments with a large hydrogen chamber (L, Alvarez), and on a combination of a bubble chamber with a hodoscopic 
counting system, A paper by Glaser (Berkeley) compared various types of bubble chambers, 


V. P. Dzhelepov told of a meter propane chamber in a magnetic field developed at Dubna, The report of I. V. 
Chuvilo was devoted to the identification of particles in a xenon chamber without a magnetic field, A report from 
the Brookhaven Laboratory was presented on the construction of a 20-inch liquid-hydrogen bubble chamber, A report 
from Princeton University was devoted to a rapid-cycling bubble chamber with propane and hydrogen, C, Ramm 
told of the propane chamber developed at CERN. 


The fruitful exchange of information on the newest technological developments that took place at this Confer 
ence will undoubtedly benefit the growth of the experimental basis of high-energy physics in the next few years, 
Also very important is the establishment of personal contacts between physicists and engineers working in various 
countries on the creation of the most complex apparatus for high-energy physics research, 


The Conference material will be published by the University of California at the beginning of 1961, 
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A DEVICE FOR THE MEASUREMENT AND AUTOMATIC CONTROL 
OF LIQUID DISCHARGE BY MEANS OF RADIOACTIVE RADIATION 


N. N. Shumilovskii and Yu. V. Gushchin 


4 
Translated from Atomnaya Energiya, Vol. 10, No, 1, 
pp. 93-94, January, 1961 


The solution of many problems in modern technology necessitates the development of new methods for measur 
ing and controlling various technological parameters, Among such problems are the development of methods for the 
automatic control of various industrial processes that are connected with the motion of highly aggressive and explo- 
sive liquids as well as liquids that are at high temperatures and under high pressures, 


From among the new methods, the contactless method of automatically controlling liquid discharge by using 
radioactive radiation deserves attention, One of the methods of measuring the discharge of liquids by using radioac- 
tive isotopes has been developed at the Institute of Automation and Remote Control, Academy of Sciences, USSR. 
This method is based on mechanical modulation of radioactive radiation, A multivane wheel, which is placed in the 
liquid flow, serves as the sensing element, One or several vanes of the wheel contains a radioactive isotope, which is 
pressed into the vane in such a manner that it cannot come into contact with the liquid, A lead collimator, where 
the radiation from the radioisotope is partially absorbed over a certain segment of its rotation path, is mounted onthe 
outisde of the pipe wall, The radioactive radiation receiver, which is connected to the measuring device by means 
of a cable, is located behind a screen, The Co™ isotope is used as the radioactive radiation source, Instead of this 
one, any other y-source whose y-quanta energy is sufficient to secure the necessary penetrability through the piping 
walls can be used, Of the radioisotopes that are most widely used at the present time, cs!87, Eyl, snM!8, etc, can 
be used, In the case of a different data transmitter design (with a lighter wall), B-sources can be used: se™, T?™ etc, 


As the vane wheel rotates in the liquid flow, the amount of radiation falling on the receiver per single vane 
wheel revolution changes from a given minimum to a certain maximum value due to the presence of the absorbing 
collimator in the device, The thus modulated radioactive beam falls on the receiver, for which a counter is used, 
The counter transforms the energy of )-quanta into electric energy in the shape of pulses, Through a cable, whose 
length can be equal to several hundreds of meters, electrical pulse packages arrive at the electrical measuring device 
input. The electrical measuring device serves for counting the number of pulse packages arriving from the receiver 
and for determining the total liquid discharge with respect to this count as well as for measuring the arrival frequency 
of packages, which makes it possible to determine the instantaneous liquid discharge, An electronic self-recorder 
can be connected to the measuring device, 


An experimental model of the radioactive liquid flowmeter (see figure), which was constructed at the Institute 
of Automation and Remote Control, Academy of Sciences, USSR, was investigated at the All-Union Scientific Re- 
search Institute of Glass in measuring black oil flow. In these experiments, a mechanical volume-type flowmeter 
with the counting mechanism removed, which was developed at the All-Union Scientific Research Institute of Glass, 
was used as the primary converter of black oil flow; instead of using packing glands, the device was made completely 
airtight, The tests were performed for discharges from 0,1 to 1.0 mYhr, The total discharge accuracy was main- 
tained within *0,5%, The maximum error in measuring the instantaneous discharge value was +1,5% of the scale 
upper limit, The new method made it possible to reduce the operational sensitivity threshold of mechanical flow- 
meters to a great extent, and to improve the measurement accuracy, 


At the present time, in collaboration with the Institute of Physics, Academy of Sciences, Latv, SSR, the All- 
Union Scientific Research Institute of Glass, and the Special Design Office(SKB) of the Avtoelectropribor Plant, the 
Institute of Automation and Remote Control, Academy of Sciences, USSk, is engaged in the development of a stand - 
arized design of the RZhr-1 device for the measurement and automatic control of the discharge of liquids flowing 
through a closed piping system as well as for the dosing of liquids, The device measuring section is based on stand- 
ardized URAP units, The device design will make it possible to measure the total and the instantaneous discharge, 
to dose a required amount of liquid with the possibility of adjustment to the assigned value, and to record the 
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readings by means of a self-recorder, The complete assembly 
of the standardized units, each of which performs the above- 
indicated functions, constitutes the measuring unit of the 
RZhR-1 device, In the first variant of the RZhR-1 device (the 
design of which was developed at the Special Design Office 
(SKR) of the Avtoelectropribor plant), which is designed for 
measuring black oil discharge, a specially constructed volume- 
type flowmeter using 6 -radiation is utilized as the discharge 
data transmitter, The data transmitter makes it possible to 
measure the discharge of black oil under pressures in excess 

of 6 kg/cm? in the discharge range from a few liters per hour 
to 4000 liters/hour, If a suitable data transmitter that is de- 
signed for measuring the discharge of different liquids is con- 
nected to the measuring device, no additional modifications 
of its measuring unit are necessary, 


The described method and the measuring equipment 
of the RZhR-1 device can also be used for contactless tacho- 
metry, i, e., for measuring the number of revolutions of 
rotating parts, shafts, etc, in locations that are not readily 
accessible, For this purpose, only a specially constructed 
data transmitter (primary transducer) can be used, 


Experimental model of the radioactive liquid 
flowmeter, 
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BRIEF COMMUNICATIONS 


Translated from Atomnaya Energiya, Vol, 10, No, 1, 
p. 95, January, 1961 


USSR, The Ninth Session of the Scientific Council of the Joint Institute of Nuclear Studies was held at the end 


of November in Dubno, Prominent scientists from 12 countries which are members of the Institute participated inthe 
work of this session, 


The Scientific Council heard reports concerning the most important results obtained in 1960 by Soviet and 
foreign scientists working in Dubno, 


The personnel of the High-Energy Laboratory succeeded in considerably improving the 10-Bev. synchrophasotron, 


USSR. On the initiative of UNESCO and the Academy of Sciences, USSR, an international symposium on pri- 
mary and initial processes occurring in living cells under the action of ionizing radiation was held in September in 


Moscow, Eminent biologists from USA, Great Britain, France, Belgium, Czechoslovakia, and the USSR took part in 
this symposium, 


The participants at the symposium brought forward a new point of view concerning the processes of radiation 
action on living cells while paying attention to injuries sustained by internal cell structures, In the reports by P. 
Alexander (Great Britain), Z, M, Bach (Belgium), and A, G, Pasynskii (USSR), it was emphasized that injuries to in- 
ternal cellular membrane structures greatly affect the metabolism process, The fonizing radiation action results not 


only in direct injuries to the cell nucleus, as was maintained for a long time by biologists, but also to other cellular 
organic components, 


Two basic methods for radiation protection were proposed recently, The first (prophylactic) method consists 
in introducing special protective substances into the organism before irradiation, These substances receive the main 
radiation shock and protect the organism. For instance, Prof. Z, M. Bach recommends a substance which reduces the 
action of radiation by 50%, The other method is the so-called substitutiontherapy, which consists of introducing into 
the organism unirradiated tissue which is most affected by the action of radiation, However, these two methods 
have not yet been sufficiently studied, and they do not always yield satisfactory results, 


USSR. A group of collaborators of the Central Scientific Research Laboratory of the Construction Ministry of 
the Latv, SSR succeeded in producing heavy concrete which is characterized by great density, 


This shielding material has been successfully used in the nuclear reactor under construction in Salaspils, 


USSR, Radioactive devices for measuring the thickness of pastes that are applied to various materials are now 
being utilized at the Kalinin Iskozh combine, 


These thickness meters provide the possibility of accurately controlling the paste layer thickness, which pro- 
vides a yearly saving of several millions of rubles, Until recently, paste layers of maximum thickness had to be ap- 
plied if the quality of leather substitutes was not to suffer. Equipment for such measurements was not available be- 
fore the development of radioactive thickness meters, 
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Books and Symposia 


Issledovanie kriticheskikh parametrov reaktornykh sistem [A Study of critical parameters of reactor systems]. 
Symposium. Moscow, Gosatomizdat, 1960, 118 pages, 36 kopeks ° 


The symposium contains original articles not previously appearing in printed form and dealing predomi- 
nantly with theoretical calculations and critical parameters in various reactor systems, 


Graphs and tables are given to establish the dependence of the critical parameters on the relative concentration 
and nature of fissionable material and moderator, as well as on fuel enrichment for a wide range of neutron energy 
spectra, 


The book is compiled for engineers and physicists working in the field of reactor system design, for experimen- 
tal physicists, and for specialists in the related fields dealing with processing of uranium- and plutonium-containing 
products, 


M, Gauzy, T, Kahan, Upravlenie yadernymi reaktorami [Nuclear reactor control], Translated from the French, 
Moscow, Atomizdat, 1960. 174 pages, 88 kopeks, 


The book contains the first two chapters of volume II of the French three-volume edition “Introduction a la 
Physique Nucleaire," 


The first chapter (chapter VII in the original) presents a survey of the basic characteristics of different particle 
detector types, The general physical operating principles of detectors are explained, and a concise description of 
accessories is given, 


The second chapter (chapter VIII in the original) takes up the problem of control of nuclear reactors, Tech- 
niques of reactor power measurement, control devices, steady state control, reactor start-up, and emergency scram 
systems are discussed, 


The translation will be found useful by engineering and technical workers in nuclear industry, 


Spravochnik po korrozii i iznosu v yadernykh reaktorakh s vodyanym okhlazhdeniem [Handbook of corrosion 
and wear in water-cooled reactors], Translated from the English. Moscow, Atomizdat, 1960, 402 pages, 1 ruble, 
89 kopeks, 


The handbook is devoted to corrosion and wear of structural materials where water is used as coolant in nuclear 
power reactors. 


The handbook consists of three parts, Part I presents general information on iron and steel corrosion effects, 
friction and wear, and takes up questions related to water technology. Part II describes techniques for testing resist~ 
ance to wear and corrosive attack on the part of some materials and combinations of materials, with tabulated data 
and discussion of various factors affecting corrosion and wear, Part III is devoted to special research on materials 
selection and development of concrete reactor designs, 


The handbook is compiled to service a broad readership of scientific workers, engineers, and technicians inter- 
ested in reactor design. 


Yu. M. Dymkov. Uranovaya Mineralizatsiya rudnykh gor [Uranium mineralization in the Kru’né Hory (Erzge- 
birge)]. Moscow, Atomizdat, 1960, 104 pages, 34 kopeks. 


Translator‘s note: all prices listed are given in terms of the new “heavy ruble,"” 
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The book is based on an analysis of factual material and theoretical concepts published in the Soviet and for- 
eign press, Information on vein formations and mineral associations of endogenic occurrences is presented for one of 
the world's oldest mining regions, the Krusné Hory [or Erzgebirge, between Saxony and Bohemia, SE of Dresden, NW 
of Prague — translator's note]. 


The book deals mainly with the genesis of uranium occurrences and first of all with the evolution of the hypo- 
genetic mineralization of the Krusne Hory. A good deal of space is given to the mineralogical characteristics and 
geochemical analysis of discrete stages of mineralization of various endogenic occurrences in the Krusne Hory. 


Literature data embracing some 160 titles are supplemented by findings from the author's investigation of ex- 
tensive mineral collections, The text is illustrated with original sketches and microphotographs, 


The book is written for geologists, mineralogists, and geochemists engaged in the study of uranium deposits, 
and for a wide audience of specialists interested in the theory of ore formation. 


Sbornik rabot po nekotorym voprosam dozimetrii i radiometriiioniziruyushchikh izluchenii | [Symposium 
of papers on some aspects of dosimetry and radiometry of ionizing radiations], Moscow, Gosatomizdat, 1960, 188 
pages, 1 ruble, 


The symposium contains 23 articles all by Soviet authors, Part I considers particular problems in the dosimetry 
of gamma and beta radiations, Part I] deals with neutron dosimetry, primarily with neutron dosimetric instruments, 
Part Ill discusses the general problems encountered in the radiometry of aerosols and gases, plus sume particular prob- 
lems of interest in the area of radiometric monitoring. 


A, V. Bibergal', V. U, Sinitsyn, N. I, Leshchinskii. Izotopnye gamma-ustanovki [Isotope gamma facilities]. 
Moscow, Atomizdat, 1960, 156 pages, 45 kopeks, 


The book deals with basic questions in irradiation techniques, It discusses sources of ionizing radiations for 
isotope gamma radiation facilities, presents examples of various types of gamma facilities, and specifies the require- 
ments applying to them, Structural components of gamma facilities, such as irradiators, systems for positioning ob- 
jects and irradiator units, radiation shielding, charging of irradiators, controls, etc, are described, Techniques used 
in calculating the dose field of sources of different geometry are explained, Design requirements affecting gamma 
facilities are listed, and examples of shielding calculations for gamma facilities are given. A list of reference litera- 
ture appears at the end of each chapter, 


The book is written for engineering and technical workers interested in problems of irradiation techniques, 


Pribory i metody analiza izluchenii, No. 2 [Instruments and techniques in radiation analysis]. Moscow, Atomiz- 
dat, 1960, 168 pages, 49 kopeks. (Moscow Engineering and Physics Institute. ) 


17 papers by the faculty and graduate students at the Moscow Engineering and Physics Institute, delivered at 
the 1957 conference of the Institute, are published in the symposium. Of these, four papers are devoted to questions 
directly related to the design of shielding devices, Four papers take up dosimetry problems regarding ionizing radi- 
ation, The last four articles describe a facility based on delayed coincidences designed to measure 107 to 107"sec 
time intervals; other topics are the determination of small phases in pion-nucleon scattering and problems involving 
attempts to increase the dispersion and resolving power of magnetic analyzers, 


Articles from the Periodical Literature 


I. Nuclear Power Physics 


Neutron and reactor physics, Physics of hot plasmas and controlled fusion. Physics of charged-particle 
acceleration 


Voprosy istorii estestvoznaniya i tekhn,, No, 9 (1960) 


——, p. 18-27, Reminiscences of Frederic Joliot-Curie, (Remarks at the reunion of scientists with participants 
of the session of the Bureau of the World Peace Council, devoted to the memory of Frederic Joliot-Curie: 
A. N. Nesmeyanov, D. V. Skobel'tsyn, John Bernal). 
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N. A, Figurovskii, 28-37, Reminiscences of Frederic Joliot-Curie (on the anniversary of his death), 


Doklady akad, nauk SSSR, 133, No. 3 (1960) 


A, G. Ershov et al., 554-57, Experimental investigations of constriction of an electron bunch in the 280 Mev 
synchrotron, 


Zhur, tekhn, fiz. XXX, No, 9 (1960) 


A. I, Anisimov et al,, 1009-18, A method for investigating the spatial distribution of electrons in a plasma, 
V. V. Yankov, 1019-23, On the stabilization of localized and flying plasmoids by electromagnetic waves, 
S. N, Breus, 1030-34, On the stability of a fluid current cylinder at finite conductivity. 

A, I. Gubanov, Yu, P, Lun'kin, 1046-52, Equations of magnetic plasma dynamics, 


A. I. Morozov, L, S. Solov'ev, 1104-08, On acceleration of a plasma in a coaxial tract, 


eksptl, i teoret. fiz, 39, No, 2 (1960) 


Li Ha-Yuan et al., 225-229, Study of the He® + H’ reaction. 


I. M. Samoilov, A. A, Sokolov, 257-259, Contribution to the problem of azimuthal instability of circulating 
currents, 


S. M, Osovets, 311-16, Dynamic stabilization of a plasmoid. 


I, Sh, Vashakidze et al,, 393-96. Investigation of the (n, p) reaction on the Li® nucleus, 
A. V. Timofeev, 397-99, Buildup of ion acoustic oscillations in an anisotropic plasma, 


V. P. Dokuchaev, 413-15, On the buildup of magnetohydrodynamic waves in a stream of plasma moving 
through an ionized gas, 


R, V. Polovin, 463-70, Contribution to the theory of simple magnetohydrodynamic waves, 
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V. N. Dmitriev, et al., 556-62. Energy distribution of threefold-fission fragments of U*™, 

R. V. Polovin, 657-61, On the escape of a plasma into a vacuum in the presence of a magnetic field. 

D, A, Frank-Kamenetskii, 669-79, On the natural oscillations of a bounded plasma , 

E. P. Sirotina, S, I, Syrovatskii, 746-53, Structure of low-intensity shock waves in magnetohydrodynamics, 
V. F. Aleksin, V. 1. Yashin, 822-26. Study of plasma stability by means of the generalized energy principle. 
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R. G. Denotkina et al., 1509-15. Determination of composition and dissociation constants of phosphate plutonyl 
(IV) complexes by the solubility method, 


A. E. Klygin, V. K. Pavlova, 1516-21. Spectrophotometric investigation of the complexing reaction of thorium 
with benzene -2-arsonic acid — (1-azo-2-)-1 ,8-dihydroxynaphthalene -3,6-disulfo acid ( arsenazo), 


Zhur. neorgan. khim. 5, No. 8 (1960) 


E. A. Kanevskii, G. R. Pavlovskaya, 1738-42. Polarographic investigation of sulfate solutions of hexavalent 
uranium, 


A. I. Alekperov, S. I. Zhdanov, 1743-47. Effect of anions on the reduction of uranyl ions on a dropping- 
mercury cathode, 


A. S. Solovkin et al., 1861-67, Extraction of uranyl nitrate using the di-isoamyl ether of methylphosphonic acid. 
V. B. Shevchenko et al., 1911-13, Extraction of plutonium from chloride solutions, with tri-n-octylamine. 


Zhur, prikladn. khim. 33, No. 7 (1960) 


G. V. Samsonov et al., 1661-64, Preparation of thorium carbides, 


Zhur, struckturnoi khim. 1, No. 2 (1960) 


L. V. Lipis et al., 135-44. Spectrophotometric study of the complexing processes of tetravalent plutonium in 
nitrate solutions, 


Zavod.lab. 26, No. 8 (1960) 


Yu. A. Chernikhov et al., 921-24. Successive chelatometric determination of thorium and sums of rare-earth 
metals, 


R. S. Volodarskaya, 925-27. Chelatometric method for determining thorium in magnesium alloys. 


R. S. Veitsman, 927-29. Photocalorimetric determination of zirconium in the form of a phosphoric-molybdenic- 
zirconium complex. 


L. N. Aleksandrov, 975-77. Radiometric determination of thorium in tungsten and molybdenum. 


Izvestiya akad, nauk SSSR, seriya geol. No. 7 (1960) 


A, A, Smyslov, 32-45. Significance of data on the radioactivity and thermal conductivity of minerals in 
metallogenetic investigations, 


Izvestiya vyssh. ucheb, zaved. Chernaya metallurgiya No, 7, (1960) 


V. G. Vlasov, A. G. Lebedev, 5-6, Dissociation of uranium trioxide. 


Inzhener. sbornik Inst. mekhaniki akad. Nauk SSSR 28 (1960) 


V. S. Lenskii, 97-133, Effect of radioactive radiations on the mechanical properties of solids (literature survey). 
Radiokhimiya II, No. 4 (1960) 


N. M. Adamskii et al., 400-410, Zirconium distribution in extraction of tri-n-butylphosphate, 


A. P. Ilozhev et al., 411-18. Distribution of butylphosphoric acids between aqueous solutions and a tri-n-butyl- 
phosphate phase, 


V. G. Timoshev et al, 419-25, Extractive power of neutral oxygen-containing organic compounds, 


V. I. Kuznetsov, T. G. Akimova, 426-30. Organic co-precipitants. XIV. Improved technique for co-precipitat- 
ing uranium from mineral waters, 


S. Z. Roginskii et al., 431-37, Glauconite ion exchange columns separate out radiocesium concentrates, 


S. Z. Roginskii et al., 438-45. Concentrates of radioactive cesium isotopes separated out on ferrocyanides of 
heavy solutions from solutions with high content of salt impurities, 


M. V. Vladimirova, Z. V. Ershova, 495-99. Effect of alpha radiation from polonium on concentrated sulfuric 
acid solutions, 
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A. A, Grinberg et al., 505-506, Contribution to the problem of instability constants of oxalate complexes of 
uranium, 


Trudy Inst. geol rud. mestorozhd., petrografii, mineral. i geokhim, akad. nauk SSSR No, 42 (1960) 


S. A. Brusilovskii, 58-99. Investigation of conditions governing precipitation of uranyl hydroxide from low- 
temperature aqueous solutions, 


Trudy Inst. metallurgii im, Baikova No. 4 (1960) 


M. E, Dritz et al., 74-83. Thorium-containing magnesium alloys. 


Trudy Inst, metallurgii im. Baikova No, 5 (1960) 


B. N. Melent'ev, N. S. Gertseva, 198-201. Polarography of large amounts of uranium. 


Trudy po khim. i khim. tekhnol. No, 2 (1959) 


N. N. Mironov, 257-60. Determination of thorium from lanthanum and cerium by the oxalate method. 


N. N. Mironov, 264-70, Methods for obtaining thorium-containing compounds and some physical- chemical 
techniques, 


Tsvetnye metally No. 8 (1960) 


I. L. Perlin, V. A. Fedorchenko, 88-93, Equipment of fabricating reactor fuel elements, and shielding devices 
(review of the foreign literature), 


Fizika metallov i metallovedenie 10, No. 2 (1960) 


Sh. Sh. Ibragimov, V. S. Lyashenko, 183-86, Fast neutron effects on properties of metals. 
F, P, Butra, 223-25, Neutron irradiation effects on the structure of molybdenum, 


Sh. Sh. Ibragimov et al., 316-17. Study of fast neutron effects on the characteristic temperature of iron and 
copper. 


Canad, J. Chem. 38, No. 4 (1960) 


A. Kertes, V. Kertes, 612-19. Extraction of mineral acids, Part III. Interaction between solute and solvent in 
the hydrobromic acid-water-tri-n-butylphosphate system. 
Canad, J. Chem. 38, No. 6 (1960) 


G. Caldow et al., 772-82, Infrared spectra of some uranyl compounds. 
Energia Nucl. 7, No. 9 (1960) is 

B. Brigoli, 590-96, Up-to-date distillation methods. I. 

A, Bassi. C. Corsetti, 638-48, Study of the structure of uranium metal. 


Industr. Chemist 36, No. 427 (1960) 


G. Miller, 447-53. Tantalum and niobium production. 


Industries Atomiques. IV, No. 5/6 (1960) 


—-—-, 97-104, Grenoble conference on plutonium metallurgy. 


J. Inorg. and Nucl. Chem. 14, No. 3/4 (1960) 


J. Korkisch et al., 247-50. Ion exchange in mixed solvents [in German].I. On the distribution of uranium be- 
tween alcoholic solutions of mineral acids and strongly basic Dowex-1 resin, 


P, Antal, et al,, 251-54, Ion exchange in mixed solvents [in German]. II. On the distribution of Th, Zr, and 
Ti between alcoholic mineral acid solvents and the strongly basic Dowex-1 anion exchange resin, 


J. Inorg. and Nucl. Chem, 15, No.1/2 (1960) 


R. Panzer, J. Suttle, 67-70. Uranium tetrachloride solvents with phosphorus oxychloride. 


G. Gibson et al., 110-14, Interaction between uranium trioxide and nitrogen pentoxide. 
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C. Banks, R. Singh, 125-32. Composition and stability of metal complexes with 5-sulfosalicylic acid. 
N. Isaac et al., 151-57, Extraction separation of transcurium elements from large amounts of curium. 
H. Irving, D.Edgington, 158-70, The synergism phenomenon in the extraction of actinides. 


M. Taube, 171-76. The effect of solvent polarity on the extraction separation of neptunium and uranium com- 
pounds, 


J. Korkisch, F, Tera, 177-81. Ion exchange in mixed solvents III. On the distribution of thorium chloride be- 
tween strongly basic anion exchange resins and alcoholic chloride solutions, 

G, Panson, C, Weill, 184-85, Production and crystalline structure of americium metal. 

C, Pistorius, 187-88. Infrared spectrum of Th(OH),. 

M. Lindner et al., 194-95. Half-life of Hf, 

R, Naumann, M. Michel, 195-96, Production of the long-lived isotope Ho'®, 

A, Aten et al., 198-99, Co-precipitation of astatine iodide with iodine. 


J. Nucl. Materials 2, No. 3 (1960) 


B. Blumenthal, 197-208, Structure of U-C alloys with low carbon content. 
D. Evans, G, Raynor, 209-15, Lattice periods in thorium-yttrium lattices, 
R. Blanchard et al., 216-24. Effect of hydrogen in tensile rupture in stainless steels, 


S. Gregg et al., 225-33. High-temperature oxidation of beryllium. Part II. Interaction between carbon dioxide 
and carbon monoxide. 


A. Young et al., 234-47, Plastic flow of alpha-uranium. 


L. Howe, W. Thomas, 248-60. Effect of neutron irradiation on the properties of Zircaloy-2, under tension. 

J. Lehmann, 261-68, A suggestion on the systematics of phases in uranium alloy. 

H. Newkirk, 269-73. Electron-microscopic studies of radiation damage in irradiated uranium dioxide. 
Energy 14, No. 148 (1960) 

——-—, 414-15. Fuel processing technology at the Eurochimique plant at Mol (Belgium). 

Nucleonics 18, No. 9(1960) 

D. Leeser, 68-73, Radiation effects on reactor metals. 


J. Carroll, R. Bolt, 78. Radiation effects on organic materials. 


J. DeMastry, R. Dickerson, 87-90, 132, 134, Niobium, promising high-temperature core materials, 
B. Levy, et al,, 128, 130, Obtaining hydrazine by irradiation of solid and liquid ammonia, 


Nukleonik 2 No. 5 (1960) 


R. Kiessling, 198-203. Swedish research on UO, from the viewpoint of fuel element applications. 
Nucleonika, V, No. 5 (1960) 


Z. Zagorski, 253-60, Determination of the G value of gamma-induced reactions by the method of constant- 
potential polarographic electrolysis. 


Nukleonika V, No. 6 (1960) 


L. Adamskii et al., 317-27. Simultaneous determinationof trace amounts of sodium and potassium by neutron 
activation techniques, 


W. Dtugosz, 329-39. Tracer studies of flotation processes. 


T. Urabnski, 341-53, Some problems in uranium extraction using dodecylphosphoric acid. 
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IVji Nuclear Radiation Shielding 


Radiobiology and radiation hygiene. Shielding theory and techniques. Instrumentation. 
Biofizika 5, No. 4 (1960) 


B. M. Isaev, 479-87, Physical foundations of radiobiological neutron flux experiments. 
Vestnik akad, nauk SSSR XXX, No. 9 (1960) 


A. M. Kuzin, 48-51. Hazard of increased C™ concentration in the atmosphere. 


Voenno-meditsinskii zhur. No. 7 (1960) 


L. A, Tiunov, 36-39, Prophylactics of radiation injury, using combinations of drugs and medicines (survey of 
the literature), 


I. V. Fedorov, G. V, Sazykin, 69-71, On feeding cases of acute radiation syndrome. 
Vrachebnoe delo No. 8 (1960) 


M. A. Khvoinitskaya, V. P. Pugachevskii, 93-94, The public health approach to work safety conditions where 
radioactive isotopes are in use in the metallurgical industry, 


Gigiena i sanitariya No. 7 (1960) 


Z. Ya. Gorodishcher, N. I. Mashneva, 56-60. Deactivation of drinking water containing radiophosphorus by the 
contact coagulation method, 


Izvestiya akad. nauk SSSR, seriya geograf. No. 4 (1960) 


V. G. Bogorov, B. A. Tareev, 3-10, Ocean depths and the problem of radioactive wastes disposal at sea, 
Med, radiologiya 5, No. 7 (1960) 


Yu, V. Novikov, I. P. Korenkov, 66-71. Reduction in the radioactivity level of atmospheric air as a result of 
the cessation of nuclear weapons testing. 


M. Ya. Maizelis, 92-95. Problems of medical radiology reviewed at the All-Union Conference on the assimila- 
tion of radioactive isotopes and nuclear radiations into the USSR national economy, Riga, April, 1960, 
Priborostroenie No. 7 (1960) 


M. L. Gol'din, 22-24. On safety techniques in the preparation and use of instruments based on the use of radio- 
active isotopes. 


Pribory i tekhn. eksper. No. 4 (1960) 


Yu. K. Akimov et al., 71-77, Separation of particles according to the degree of ionization in some scintillation 
counters, 


Yu. A, Egorov, D. S. Khukhorev, 136-37. On the use of light pipes in scintillation counters. 


L. 1. Zlobin, Yu. I. Lyubimov, 137-38, Transistorized preamplifier for scintillation counters, 
L. D. Soshin, 138, Transistorized multivibrator with millivolt operating threshold. 


G. N. Muskhelishvili, G. V. Zakomornyi, 139-41. Electronic voltage stabilizer for transistor circuit power 
supplies, 


Sovetskoe zdravookhranenie Kirgizii No. 2/3 (1960) 


I. P. Mashkov, 112-18, The Tashkent conference on the peaceful uses of nuclear energy (review of the medical 
panel presentation and discussion). 


Trudy Ural. otdel. Mosk o-va_ ispytatelei prirody, No. 2 (1959) 


N. V. Gorbatyuk, N. V. Timofeev-Resovskii, 181, 193-4, On critical tolerance norms for radioactive contami- 
nation of water and air, I, IL. 


L. K. Ponomareva et al., 201-05. Determination of cesium-137 in the water of open reservoirs. 
Khlopkovodstvo No. 8 (1960) 


4 
4 
102 


D. T. Kabulov et al., 48-51. Effect of radioactive isotopes on the growth, development, and crop yield of the 
cotton plant. 


Archives Environmental Health 1 No. 2 (1960) 


J, Felton, C. Rozas, 87-95. Deactivation of human skin contaminated by hot spills as an experiment, 


Atomkernergie 5, No. 7/8 (1960) 


K. Lindackers, 272-77. Design of biological shielding for reactors. 
H. Schneider, E. Schwerdtel, 278-81.Portable scintillation counter with differential discriminator. 
K. Godt, K, Sommermeyer, 282-85, RaD content in plants, determined by gamma-ray spectroscopy. 
K. Potzl, R. Reiter, 285-87. Quantitative determination of fission-product fallout in the Northern Alps. 
Atompraxis 6, No. 8 (1960) 

A. Berger, 301-308, Behavior of x-rayed potatoes. 

A. Micke, K. Woéhrmann, 308-16, Radiosensitivity of dry seeds, 


S. Krawezynski, 320-22. Effect of synthetic purifying agents and complexing agents on deactivation of radio- 
active waters (deactivation by chemical precipitation), 


Canad. J. Chem. 38, No. 6 (1960) 


W. Armstrong, D. Grant, 845-50, Aqueous benzoate system as a sensitive dosimeter of ionizing radiations. 


Canad. J. Phys. 38, No. 8 (1960) 


J. Wolfson, F. Terentiuk, 991-1010. Concentration of fission products in the air at sea level. in Canada during 
the IGY. 


Energia Nucleare 7, No. 9 (1960) 


E, Cerrai, 581-89, Problems arising in the disposal of liquid hot reactor wastes, 


E. Cerrai et al., 630-37. Clean-up of liquid wastes at the CISE radiochemical laboratory. 


Jaderna Energie VI, No. 9 (1960) 


V. Santholzer, 294-98, Increase in fallout radioactivity as a result of nuclear testing in the Sahara, 

F, Béhounek, V. Zelenkova, 299-302. Determination of beta activity of liquid wastes, 

J. Mackou, 303-11, Designs of automatic radioactivity measuring facilities, 

J. Koti, 312, Individual sonic dosimeter warns of gammas, x-rays, and neutrons, 
Energy 14, No. 148 (1960) 

A. Gould, 421-24, Planar film dosimeter and electroscope. 
Power 5, No. 54 (1960) 

R. Owen, 82-86, Scintillation counter review. 


Sci. and Engng. 8, No. 2 (1960) 


W. Pettus, 171. Age of U?%* fission neutrons in water. 

P. Khubchandani, 172, Thermal neutron scattering in graphite. 
Nucleonics 18, No. 9 (1960) 
G. Hine, J. Cardarelli, 92-94, 96, 100, Conical plastic scintillators show total gamma absorption, 
S. Sandler et al., 102-03, Dimethyl styrene yields more efficient scintillators, 
Nucleonics 18, No. 10 (1960) 


A. Jones, 86, 88, 90, 91. 


Uses of semiconductor detectors in health-physics monitoring. 


R. Tolmie, 92. Dual detectors kc; to radioisotope identifier. 
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Nukleonik 2, No. 5 (1960) 


W. Pohlit, H. Pohlit, 175-78, Paraffin moderator used for dosimetry in a mixed neutron field. 
W. Kern, 203-13. On fragment activity due to fallout. 
Nukleonika V. No. 5 (1960) 


I, Ztotowski, H. Wincel, 261-79. Mass-spectrometric studies of chemical processes occurring in a self-quench- 
ing Geiger-Miller counter loaded with heavy hydrocarbons, 


J. Liniecki et al., 301-13, Sr® content in bones of humans and animals and in milk in Poland during 1958. 
Nukleonika V, No. 6 (1960) 


S, Chwasczewski, 355-60. Pulse modulation of a high-frequency ion source. 
T. Radoszewski, 361-68, Scintillation counter characteristics. 


M. Miazek-Kula, 373-76, Deactivation of a surface contaminated by hot isotopes. 


C, Szwacka, 376-77, Deactivation of protective clothing worn in radiochemical plants. 
J. Flattau et al,, 377-78, Protective clothing for deactivation points, 


Radiation Res. 12, No. 1 (1960) 


N. Barr, M, Stark,1-4. Fluorescent compounds in chemical dosimetry. Gamma radiation effects on chinine 
fluorescence, 


M. Berger, D. Raso, 20-37. Monte Carlo computations of background gamma radiation. 
Radiation Res, 12, No. 3 (1960) 


W. Gordy, I. Miyagawa, 211-29, Electron spin resonance studies of mechanisms for chemical protection from 
high-energy ionizing radiations. 


Radiation Res, 12, No. 4 (1960) 


J. Hunt et al, 319-24, Effect of parametric gases on the electron spin resonance signal from activated charcoal. 


Radiation Res, 12, No. 6 (1960) 


S. Kronenberg, H. Murphy, 728-35, Energy spectrum of protons emitted by a hydrogen-containing substance. 
Radiation Res. 13, No. 1 (1960) 


W. Jesse, 1-17. Ionization of air by polonium alphas and mean energy of formation of ion pairs. 


S. Kondo, M. Randolph, 37-60, Effects of finite dimensions of an ion chamber on measurements of low-level 
photon sources, 


P. McClement, G. Failla, 61-91. Dose measurement in a small volume of tissue surrounding a radioisotope 
point source, 


Radiation Res. 13, No. 2 (1960) 


EF. Burke, R. Pettit, 271-85. Absorption analysis of x-ray spectra produced by beryllium-window tubes operated 
at 20-50 kv. 


Research XIII, No. 10 (1960) 


D. Stevenson, 383-89, Theoretical and practical deactivation problems. 


V. Radioactive and Stable Isotopes 


Labeled-atoms technique. Uses of radioactive radiations, Direct conversion of nuclear to electrical energy. 


Azerbaid, neft. khoz. No. 1 (1960) 


: 

104 


A. G. Khanlarova et al., 39-41, Tracer study of the dynamics of formation of a protective corrosion-resistant 
film on metals by AzNII-7 and TsLATIM-339 additives. 


Byull okeanograf. komis. pri Prezidiume akad, nauk SSSR No. 3 (1959) 


O. K. Leont'ev, V, N, Afanas'ev, 73-79. Experience in the use of radioisotopes for research on shifting of sand 
bars, 


Vestnik mashinostroeniya No. 8 (1960) 


V. I. Postnikov, 78-80. Introduction and assimilation of radioisotopes into machinery manufacture. Based on 
materials of the All-Union Conference at Riga, April, 1960, on the assimilation of radioactive isotopes 
and nuclear radiations into the USSR national economy, 


Vestnik sel’sko-khoz. nauki No. 7 (1960) 


V. N. Silin, 115-18. A powerful tool for the progress of science and industry (on the uses of radioative isotopes 
in agriculture), Contribution to the balance sheet of the Riga April, 1960 All-Union Conference on the 
assimilation of radioactive isotopes and nuclear radiations in the USSR national economy. 


Vinodelie i vinogradarstvo SSSR No. 5 (1960) 


P, I. Litvinov, 37-40. Tracer applications in the study of the regeneration of grapevine roots and in elucidating 
the relations between separate roots and vine shoots. 


Vodosnabzhenie i san. tekhnika No. 7 (1960) 


I, I. Mechitov, 1-5, Tracer applications in hydraulics research on the performance of sedimentation tanks. 


Geologiya i geofizika No. 3 (1960) 


E. M. Filippov, 103-08. Contribution to the problem of deepening research on minerals and ores via scattered 
gamma radiation techniques, 


Geol. nefti i gaza No. 7 (1960) 


N. A, Vasil'eva et al., 55-59. Experience in the use of tritium, the radioactive hydrogen isotope, to study the 
flow of injected water driven into an oil reservoir. 


Zavod. lab. 26, No. 7 (1960) 


T. M. Khrenkova, 833-35, Radiometric technique for measuring coating thicknesses, 


I. A, Berman, S, A. Khaletskaya, 836-37, Determination of the thickness of zinc coating by the tracer method. 


Izvestiya akad, nauk SSSR, seriya geograf. No. 4 (1960) 


Yu. N. Agashkin, 117-21. Tracer studies of snow thawing and runoff. 
Izvestiya vyssh, ucheb, zaved, Gornyi zhurn No, 6 (1960) 


B. M, Barbin, 139-42, Study of flotation feasibility of copper by the null method and by the tracer method. 
Izvestiya vyssh. ucheb. zaved., Lesnoi zhur. No. 2 (1960) 


A, S, Freidin, I. G, Krapivina, 117-20, Sterilization of lumber by ionizing radiations. 


Izvestiya vyssh. ucheb, zaved. Pishchevaya tekhnol. No. 2 (1960) 


V. F. Oreshko, K. A. Korotchenko, 17-21, Gamma-ray techniques in the study of starch degradation. 


Izvestiya vyssh, ucheb, zaved. Khimiya i khim. tekhnol. 3, No. 1 (1960) 


M. M. Grudinina, E, M. Aleksandrova, 176-78. Some problems related to emulsion polymerization of styrene, 
using the tracer method, 


Izvestiya vyssh. ucheb, zaved. Khimiya i khim., tekhnol. 3, No. 3 (1960) 


G. B. Fedorov et al., 395-401. Tracer study of the mobility, interatomic coupling and element distribution in 
zirconium and zirconium alloys, 
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Koks i khimiya No. 7 (1960) 


A. Z. Kulishenko, K. P. Medvedev, 5-10. Radiosulfur S® as tracer in studying desulfurization of coals. 


Neft. Khoz. No. 6 (1960) 


I. G. Zhuvagin, Yu. A. Akchas'yanov, 7-12, A new method using radioactive isotopes for checking the results 
of hydraulicking operations. 


Prikladnaya geofizika No. 26 (1960) 


A. A. Korzhev, 137-48. Contribution to estimating the porosity of permeable rocks, by the tracer method. 
Radiotekhnika i elektronika 5, No. 7 (1960) 


G. N. Shuppe, V. P. Vasil'ev, 1135-44, Tracer application to the study of oxide cathode processes and other 
problems in cathode electronics, 


Svarochnoe proizvodstvo No. 9 (1960) 


G. V. Balabina, 16-17. Thulium-170 tracer monitors weldments. 


Trudy vsesoyuz. nauchno-issled, inst, mekhanizatsii sel'sko-khoz. 28 (1960) 


A. Khamidov, 94-107. Use of gamma radiation in studying the effects of wheels on soil. 


Trudy vsesoyuz. neftegaz. nauchno-issled. inst. No. 29 (1960) 


F, A. Akelseev, F. Ts. Denisik, 32-43, Radioactive techniques for monitoring the development and exploration 
of oil fields. 


N. A. Vasil'eva et al,, 266-77, Results of tritium tracer investigations of injected water movement in an oil 
reservoir, 


Trudy naucho-isslef. inst. rezinov. prom., symp.#6 (1960) 


L. A. Oksent'evich, E. V. Zhuravskaya, 102-110. Nuclear radiation effects on crude and finished rubber. 
Khim. volokna No. 3 (1960) 


S. A. Nechaeva et al., 7-9. Study of the possibility of enhancing the heat stability of polyolefin fibers by ioni- 
zing radiation processing, 


Khimiya i khim, tekhnol. 3, No. 3 (1960) 


A. D. Sotskov et al., 452-56. Tracer studies of diffusion processes accompanying chemical and phase transfor- 
mations, 


Atompraxis 6, No. 8 (1960) 
N. Desrosier, 293-300, The present state of the art in radiation processing of foods. 


Industr. Chemist 36, No, 427 (1960) 


P. Connor, W. Hardwick, 427-33, Use of radioactive radiations for sizing particles. 


Industries Atomiques IV, No. 7-8 (1960) 


C. Kellershohn, P. Vernejoul, 52-61. Medical tracer applications. 
J. Juillard, 75-78, A plasma diode thermocouple. 
J. Brit. Nucl. Energy Conf. 5, No. 3 (1960) 


M. Crowley-Milling, 207-25, Radiation applications in industry. 
Nova Technika No. 6 (1960)[Czech]. 


Z. Spurny, 270-71, Opportunities for using radioactive wastes and high-level radiation sources in the national 
economy. 
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Nucl, Energy 14, No, 147 (1960) 
—--—, 369-70, New Amersham hot isotope laboratories, 


W. Huber, A. Klein, 371-72, The outlook for radiation processing of foods. 
Nucleonics 18, No. 8,(1960) 


——-—, 94, 96, 98, Review of papers heard on panels on hot isotopes and radiations at the 1960 ANS meeting 


F, Jackson et al., 102-105. Efficient tritium labeling with an electric discharge. 

R. Black, W. Kerwick, 106-108, Sealed millicurie floats detect pipe leaks, 
Nucleonics 8, No. 10 (1960) 

G. Eichholz, 116, 118-120, Industrial tracer techniques in Canada, 

B. Schumacher, 106, 109-110, 112, 114, Gaging gas density with fast charged particles, 
Nukleonika V. No. 5 (1960) 


J. Domanus, B. Osuchowski, 281-300. High-level Co™ sources used in gamma radiography of metal products. 
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SIGNIFICANCE OF ABBREVIATIONS MOST FREQUENTLY 


FIAN 

GDI 

GITI 

GITTL 

GONTI 
Gosenergoizdat 
Goskhimizdat 
GOST 

GTTI 

IL 

ISN (Izd. Sov. Nauk) 
Izd. AN SSSR 
Izd. MGU 
LEIZhT 


Stroiizdat 
TOE 
TsKTI 
TsNIEL 
TSNIEL~MES 
TsVTI 

UF 
VIESKh 
VNIIM 
VNIIZhDT 
VTI 

VZEI 


Note: Abbreviations not on this list and not explained in the translation have been transliterated, no further 
information about their significance being available to us. — Publisher. 


ENCOUNTERED IN SOVIET PERIODICALS 


Phys. Inst. Acad, Sci. USSR. 

Water Power Inst. 

State Sci.-Tech. Press 

State Tech, and Theor. Lit. Press 

State United Sci.-Tech. Press 

State Power Press 

State Chem, Press 

All- Union State Standard 

State Tech. and Theor, Lit. Press 

Foreign Lit. Press 

Soviet Science Press 

Acad. Sci. USSR Press 

Moscow State Univ. Press 

Leningrad Power Inst, of Railroad Engineering 
Leningrad Elec, Engr. School 

Leningrad Electrotechnical Inst. 

Leningrad Electrical Engineering Research Inst. of Railroad Engr. 
State Sci.-Tech. Press for Machine Construction Lit. 
Ministry of Electrical Industry ¢ 

Ministry of Electrical Power Plants 

Ministry of Electrical Power Plants and the Electrical Industry 
Moscow State Univ. 

Moscow Inst. Chem. Tech. 

Moscow Regional Pedagogical Inst. 
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A Milvstone in Seiotifie Communication... 


Advances in | vowne: 
| VOLUME 2 
Cryogenic | VOLUME 3 
VOLUME 4 
Engineering VOLUME 5 


Edited by K. D. Timmerhaus 


These 5 volumes, the proceedings of the Cryogenic Engineering 
Conferences held since 1954, form the only authoritative collection on the 
research and development to date in this new and exciting field. 


~ 


A working tool for industry! 
ROCKETRY + MISSILES + LIQUEFACTION OF GASES + FREE RADICALS 
SPECTROSCOPY ELECTRONICS MINIATURIZATION AIR SEP- 


ARATION + PURIFICATION + INSULATION -+- ATOMIC ENERGY 


are all vitally affected by the rapid research progress in low-temperature 
technology. 


In order to meet the growing demand for more information about 
cryogenics, Plenum Press has reprinted Volumes 1-4 in new hard-cover 
editions and presents 


(Special rate VOLUME 5 a ae 
of $35 for (Proceedings of the Fifth National Conference on s1840 


Volumes 1-4 


Cryogenic Engineering, 1959) (set) $40.00 


first four velumes) 


for the first time. 
cloth illustrated - ever 350 pages in each volume $13.50 per volume 


PLENUM PRES S 227 west 17TH STREET. NEW YORK 11, N Y¥ 


A subsidiary of Consultants Bureau Enterprises, Inc. 
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RESEARCH by Soviet EXPERTS 


Translated by Western Scientists 


SPECTRA AND ANALYSIS 
by A. A. Kharkevich 


The first handbook directed toward acousticians and others work- 
ing in those fields which require the analysis of oscillations— 
ultrasonics, electronics, shock and vibration engineering. This 
volume is devoted to the analysis of spectral concepts as they are 
applied to oscillations in acoustics and electronic engineering, and 
to a discussion of the methods of spectral analysis. Contents in- 
clude KOTEL’NIKOV’S theorem for bounded spectra, the spec- 
tra and analysis of random processes, and (in connection with the 
latter) the statistical compression of spectra. 


cloth 236 pages $8.75 


ULTRASONICS AND ITS 
INDUSTRIAL APPLICATIONS 
by O. |. Babikov 


This work is concerned with ultrasonic control methods which are 
applied in industry, and also with the action of high-intensity 
ultrasonic oscillations on various technological processes. Con- 
siderable attention is devoted to ultrasonic pulse methods of flaw 
detection and physicochemical research. It is an invaluable aid 
to scientific researchers, engineers, and technicians working in 
fields which make use of ultrasonic methods industrially, as well 
as being a convenient reference for a broad category of readers 
who might wish to become acquainted with the current state of 
ultrasonic instrumentation. 


cloth 265 pages $9.75 


Tables of contents upon request 
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